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Vibrational Frequencies of Semirigid Molecules: 
A General Method and Values for Ethylbenzene'’ 


By William J. Taylor * and Kenneth S. Pitzer * 


It is shown that in normal coordinate calculations of the vibrations of molecules it is 


sometimes advantageous to calculate the kinetic-energy matrix, rather than the reciprocal 


kinetic-energy matrix. 
matrix. 


Explicit formulas are given for the elements of the kinetic-energy 
Illustrative calculations are given for propane, toluene, and ethylbenzene. A 


semiempirical assignment of the vibration frequencies of ethylbenezene is made on the 


basis of these calculations and the spectroscopic data. 


I. Introduction 


For the purpose of calculating the thermo- 
dynamic functions (particularly the entropy) of 
molecules by statistical methods, it is important 
that the low frequencies of vibration be assigned 
When these frequen- 
cies have not been observed in the spectra, the 
assignment must be made either by analogy witb 


as accurately as possible. 


other molecules, or on the basis of a normal 
force constants 
Even when 
the frequencies have been observed, a preliminary 
assignment by analogy or calculation is often 
necessary for the identification of the observed 


coordinate calculation, using 


determined from other molecules. 


frequencies. However, a complete normal coor- 


dinate analysis for molecules with many atoms is 
so time-consuming as to be impracticable in most 


investigations. Therefore, an approximate meth- 


od of calculating the low skeletal frequencies of 
large molecules should be of considerable value in 
connection with thermodynamic calculations. 


his investigation was performed at the National Bureau of Standards 
it of the work of the American Petroleum Institute Research Project 
the “Collection, Analysis, Calculation, and Compilation of Data on 
Properties of Hydrocarbons” 
‘esearch Associate on the American Petroleum Institute Research 
t 44 at the National Bureau of Standards 
seciate Supervisor, American Petroleum Institute Research Project 
fessor of Chemistry, University of California, Berkeley, Calif 
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In this paper there is described an approximate 
method that is based on the replacement of the 
actual molecule by a semirigid model (of the same 
dimensions and masses) in which certain groups 
of atoms move as rigid units. In this way the 
number of internal coordinates of the molecule 
may be reduced sufficiently to make the caleula- 
tion of the low skeletal frequencies practicable. 
The semirigid model should be 
provided that all the frequencies of internal 
motion of the groups assumed to be rigid are 


satisfactory, 


either considerably higher than the frequencies 
to be calculated, or else do not interact seriously 
with them because of molecular symmetry or 
other factors. In order to reduce the error due 
to interactions, it is desirable that the necessary 
determined by 
(semirigid model) calculations on other molecules 


force constants be analogous 
in which the frequencies are known. 

The analysis given in this paper shows that 
the kinetic-energy matrix of a semirigid molecule 
is usually more easily calculated than the reciprocal 
kinetic-energy matrix. The reverse is true for a 
complete normal coordinate calculation for the 
molecule. Formulas are given that reduce the 
calculation of the kinetic-energy matrix to a 
consideration of the masses, moments of miass, 
and moments of inertia of the rigid groups. 








An illustrative calculation is given for ethyl- 
benzene by using force constants determined 
from propane and toluene. A semiempirical 
assignment of the vibrational frequencies of 
ethylbenzene is made on the basis of this calcula- 
tion and the spectroscopic data. These calcula- 
tions for ethylbenzene were completed in September 
1945, were and used in the tables of the thermody- 
namic properties of ethylbenzene issued November 
30, 1945, by the American Petroleum Institute 
Research Project 44. 


II. Mathematical Analysis 
1. Kinetic Energy of Semirigid Molecules 


The first step in the normal coordinate analysis 
of a vibrating system is the computation of the 
kinetic-energy matrix, or of its reciprocal. Certain 
general features of the calculation will first be 
considered. Let the system of n particles be 
characterized by 3n generalized coordinates, q, 
six of which refer to over-all translation and rota- 
tion of the system (external coordinates), and the 
remainder of which specify the configuration of 
the system (internal coordinates). The present 
discussion is concerned only with small distortions 
of the system from its equilibrium configuration, 
so that the internal coordinates may pe assumed 
to measure this distortion and to have the value 
zero at the equilibrium configuration. 

The kinetic energy is a quadratic form in the 
coordinate velocities, the q's, 


(1) 


In this equation g is a column matrix, ¢ is a row 
matrix (the transpose of g), and S is the square 
symmetric kinetic-energy matrix. In general, 
the elements of S are functions of the q’s (not of 
the g's), but because of the restriction to configura- 
tions near the equilibrium configuration, S is to 
be regarded as a constant matrix here. 

For reasons that will appear presently, it will be 
convenient to assume that the coordinates are 
of three types, and to partition the column 


2T=GS4. 


matrix g accordingly. 
qi 
q=| q (2) 
qs 
The kinetic-energy matrix, S may be partitioned 
into the corresponding submatrices, and it will 
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also be convenient to designate by G the recipro: al 
of S and to partition G, 


S; 1 Si Sis G, 1 G, 2 Gi; 
S = Sw Sx» Sx ’ G = S- : G2 Gy Ox } 
S; 3 Soy Sos G 13 6. G. 33 


The concept of a semirigid molecule, with 
which this paper is principally concerned, will 
now be introduced. A semirigid molecule, as 
defined here, is one in which some of the internal 
coordinates, say the group q;, remain 
throughout the molecular motion, so that q,;=0, 
and q;=0. Thus, if parts of the molecule move as 
rigid units in their equilibrium configurations, the 
internal coordinates of these parts will remain 
zero. 

The kinetic-energy matrix for a semirigid mole- 
cule reduces to 


zero 


[= Sx ; (4 
Si Sy 
as q,;=0. The equivalent expression for the 
kinetic energy is 

27 = Sgr + 2G: S292 + G2Snge- (5 


The form of the potential-energy expression 
will now be considered. 
nates do not influence the variable part of the po- 
tential energy (the assumption that q,=0 is 


equivalent to setting the force constants for these | 
In addition, there | 


coordinates equal to infinity). 
will be a group of coordinates, say gq, that do not 
occur in the expression for the potential energy, 
although g; occurs in the kinetic energy. This 


group of coordinates of zero frequency includes | 
the external coordinates and possibly (as an ap- 

proximation) some very low frequency internal | 
Final- 5 
ly, it will be assumed that the potential energy is 4 | 
remaining coordinates ) 


coordinates, in particular torsional motions. 


quadratic form in the 
(harmonic vibrations) so that, 


2V=GF mp, (6 


where Fy 
matrix. 

If the Lagrangian equations of a 
for which the kinetic and potential energies ar 


is the symmetric 


given by eq 5 and 6, respectively, are to be 


integrable, it is necessary to remove by a suitabl: 


As q;=0, these coordi- | 


potential-energy | 


system | 
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transformation the cross term between q, and q» 
in the kinetic energy. The existence of this cross 
term is a consequence of the fact that the initial 
coordinates, q have been defined in such a way 
that their conjugate momenta, p, =(O7'/0q,), are not 
independent of q@. As a result, the momenta p, 
are not conserved, as they should be, because the 
coordinates q, do not enter into the potential- 
energy expression, and are therefore not acted 
upon by any forces. 

As it is desired to find the proper kinetic-energy 
matrix for the (arbitrary) internal coordinates q, 
it follows that the transformation should leave 
these coordinates unchanged. The most general 
linear transformation of this type is 


a] T), | qi |. (7) 
q2 0 Ex 2 

where #4 is the unit matrix with the appropriate 
number of rows and columns. The matrices 7), 


and 7, are subject only to the condition that the 
transformed kinetic-energy matrix have the form 


Sy 0 

, Ss 

[ &: | 9) 

The equivalent expression for the kinetic energy is 


2T ‘sg T G2S2x4]2- (9) 


The transformed kinetic-energy matrix is obtained 
by multiplying the original kinetic-energy matrix 
from the right by the reciprocal of the transforma- 
tion matrix, and from the left by the transpose of 
the reciprocal (congruent transformation). There- 
fore, 

Si, Tr" 


0 —T,.T 3 
Su Sto Tx" = T,5"1 2 
- Ge! ie Ex 


Equation 10 is equivalent to the relations 


~ 


Ni'SuTr" 


Ti7'(S2—SuT,7'Ti2) 
(Si ~ Te 7, Su) Tir 


S»—T2Ti7'S2—S2Ti7'T2+ 
T eT SuTiz'Ti2 
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Equations 12 and 13 are equivalent, and yield 
q 1 7 
Tir'Te2=Sir'Sx (15) 


as the condition that there be no cross term be- 
tween g; and q in the kinetic energy. 

Equation 15 does not determine 7, or Ty 
separately, but only the product 7\;7'Ty, and 
therefore Sj, is not uniquely determined. How- 
ever, Sy» is unique, as the combination of eq 14 
and 15 yields 


~ 


Sx = Sx —SpSir'Si2 (16) 
The matrix S;, will be called the reduced kinetic- 
energy matrix for the internal coordinates (exclud- 
ing those that do not occur in the expression for 
the potential energy). A result equivalent to 
eq 16 has been obtained by Majantz [1]* by an 
argument less fundamental than that presented 
here. 

When applied to the external coordinates, the 
transformation of eq 7 is equivalent to the applica- 
tion of the conditions that the linear and angular 
momenta of the entire molecule are independent 
of the q’s. When applied to an internal tor- 
sional motion, the transformation is equivalent to 
the condition, previously stated by Sayvetz [2, 3], 
that the momentum conjugate to the torsion 
should be independent of the (remaining) q's. 

An important property of transformations of 
the type of eq 7 may now be stated. The proof 
is given in appendix 1. The coordinates g,; may 
be divided into two or more groups in various 
ways, and the transformation applied consecu- 
tively to each group. The over-all transforma- 
tions from q, to g, will not, in general, be identical. 
However, the theorem states that all possible 
transformations of this type lead to the same 
reduced matrix S.. The importance of this 
theorem is that in practice it is convenient to 
apply the transformation consecutively to trans- 
lation, rotation, and finally, in some cases, to 
certain internal coordinates. 

The Lagrangian equations for the system may 
now be written as 


A=Sugi=0 (17) 


SxoG2 + Fnq,=0 (18) 


‘ Figures in brackets indicate the literature references at the end of this 
paper. 








The first equation shows that the momenta p, are 
conserved. The second equation may also be 
written 
Gis + GsFrog (19) 
where 
(20) 


is the reduced reciprocal kinetic-energy matrix for 
the internal coordinates q. In eq 16, Sy» is de- 
fined in terms of the submatrices of S. It is 
proved in appendix 2 that G,, may be expressed 
in terms of the submatrices of G, eq 3, by the 
equation 


Gig = (Sa2) 7! = Gn — Gg Gs ‘Gay. (21) 


This result has been obtained previously by 
Wilson [4], although not explicity for the case in 
which coordinates of the types q, and q, are present 
simultaneously. 

As is well known, eq 18 and 19 lead to the 
following secular equations for the determination 
of the normal frequencies of vibration of the 
coordinates q.: 


F’..— XS,.| =0, 22) 


GF — Eel =0, (23) 
where \=42°*(v=normal frequency). In_ the 
present method, in which S. rather than G» is 
computed, it is convenient to convert eq 22 to the 
secular equation of a single matrix, in the form 


Si.Fez'—d—'Ex|=0. (24) 


In practice, Fy» is diagonal, or nearly so, and 
therefore easily reciprocated. The roots of eq 
24 are, of course, the reciprocals of the \‘s. 

It is possible now to state clearly the relative 
advantages of calculating S.. and Gy» (see eq 3, 16, 
and 21). 


In calculating S.o, the coordinates q, are 
ignored, but all the elements for g, and qg, must be 
Conversely, in calculating Gy», the 


computed. 
coordinates g, are ignored, but all the elements for 
gq: and g,; must be computed. It should also be 
noted that even with the simplifications deseribed 
later in this paper, the elements of S are more 
difficult to compute than those of @° (this state- 

Wilson [4] has described a now well-known method for the calculation of 
G A detailed exposition of the method has been given by Meister and 
Cleveland [5]. A method similar to Wilson's has been developed indepen- 
dently by the Russian workers M. Eliashevich, B. Stepanov, and M. Wolk- 


enstein. References to the papers of these workers are not given here, as we 


ire not directly concerned with the calculation of G 
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ment does not apply to the elements for le 
external coordinates, but these are just the «le 
ments of G that may be ignored). From those 
considerations it is obvious that the caleulat on 
of G.. is to be preferred when all or most of the 
internal coordinates are included in the calculation 
(so that the number of coordinates of type q, is 
zero or small). On the other hand, S.. is often to 
be preferred in calculations by the semirigid mole- 
cule approximation as the number of coordinates 
of the type q, is then relatively large, and the 
number of coordinates remaining in the calculation 
(type g.) is relatively small. 

If the molecule has the symmetry of one of the 
point groups, the internal coordinates fall into sets 
equivalent under the group. The coordinates in 
an equivalent set may be transformed into each 
other by the operations (rotations and reflections 
of the group. In this case, just as is the case when 
working with the reciprocal kinetic-energy mtarix 
[4, 5], it is desirable to construct new coordinates, 
symmetry coordinates, which are linear combina- 
tions of equivalent internal coordinates, and which 
reduce the group. That is, the symmetry co- 
ordinates are chosen in such a way that each trans- 
forms according to one of the irreducible repre- 
sentations of the group. The advantage of the 
use of symmetry coordinates is, that there are no 
interaction elements in the kinetic- or potential- 
energy matrices between coordinates which trans- 
form according to different irreducible representa- 
tions. The matrices therefore factor into sub- 
matrices (the number of identical submatrices for 
each representation is equal to the dimension, or 
degeneracy, of the representation), and the 
secular equation factors into the corresponding 
equations of lower degree [{6, 7]. 

It will that the symmetry 
coordinates may be constructed either before or 
after applying the transformation of eq 7, provided 
that none of the symmetry coordinates involv 
linear combinations of coordinates of types q, and 
q2. With this restriction the symmetry 
dinates may be represented by 


Uy 0 q1 
0 Ux qoy’ 


The matrices U’,, and [7 are usually, although not 
necessarily, orthogonal. The transformed kineti 
energy matrix is 


now be shown 


coor- 
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vhich is identical with the result obtained by 
ransforming the original matrix S, of eq 16 by 
i’. It will also be observed that the matrix 
S is unaffected by symmetry combinations 
ormed with the coordinates g, (such, for example, 
is a rotation of axes). 

When (’» is orthogonal the transformed kinetic- 
nergy matrix is U’,,Sj,0», but for nonorthogonal 
» the general form 0,;'S,,U,z' must be used. 
An illustration of the use of a nonorthogonal 
ransformation will occur in the calculations at 


he end of this paper. 


2. Calculation of the Kinetic-Energy Matrix 


(a) Equations for Generalized Coordinates 


Explicit formulas for the calculation of the ele- 
vents of the kinetic-energy matrix for generalized 
iternal coordinates are presented in this section. 
Jenote by m, the mass of the ith atom, and by 

the position vector of the ith atom with respect 
0 & point stationary in an inertial system. Then 
he kinetic energy of the system of particles is 


2T=DSm R,R,, 28) 
i 


here R, is the vector velocity of the ith atom. 

lt will be assumed here, as in the preceding 
ction, that the positions of the n particles are 
pecified by 3n generalized coordinates g, some of 
hich, as before, remain zero for a semirigid 


wlecule.6 We may then define a set of vectors, 


Smi= (OR,/Oq»). (29) 


is section specific coordinates will be indicated by letter subscripts, as 
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Thus §,,,is a vector pointing in the direction in 
which the ‘th atom moves when g,, increases from 
zero (for the equilibrium configuration) to a small 
positive value, all other q’s remaining zero, and the 
magnitude of s,,; is equal to the displacement of 
the ith atom per unit change in g,. Then, 


R= SOS ndms 


m 
and 
> ( 28m Gn )-( LS. ) 
1 


m n 


4 mdn( SOM Sime Sri 
! 


m,n 
But from eq 1, 
2T=) Swudnde, 


min 


where S,,, is the element in the mth row and nth 
column of the symmetric kinetic-energy matrix 
S. It follows from eq 31 and 32 that 


,, SOM S» "Bi, 


Sinn 


¥ =~ 
S, n= > \M Berd, 
1 


In order to proceed further, define an orthogonal 
system of x, y, 2,-axes in the molecule, and denote 
the unit vectors, directed positively along the 
axes, by x, y, and z Let the system be right- 
handed, so that xxy=z. The position of the 
origin of the system, and the directions of the 
unit vectors, must be uniquely determined with 
respect to the molecule for any arbitrary (small) 
displacements of the internal coordinates, but the 
choice of axes is otherwise arbitrary. In parti- 
cular, it is not necessary that the origin be at the 
center of mass of the molecule. 

The coordinates of the ith atom are xz,, 
24. For a given choice of axes these coordinates 
are completely determined when the values of the 
given. The position 


¥;, and 


internal coordinates are 
vector of the ith atom with respect to the origin is 


Fy, X2,+ YYit Z2;. (35) 


Three of the generalized coordinates will be 
taken to be the translations of the origin of the 
axes in the (instantaneous) directions of the x-, y-, 
and z-vectors. The r-vectors of eq 29 are identi- 
cal for each atom, for each translation, and are 
simply the unit vectors 








Three more coordinates will be taken to be the 
three rotations, by angles ¢, ¢, and ¢, about the 
r-, y-, and 2-axes. The S-vectors for these are easily 
shown to be the vector products, 


Sori X XP, Soy Y XTi, Soe BAXTi.- (37) 


For the remaining coordinates, the (3n—6) in- 
ternal coordinates, 

Sm, (OF ,/ Og m) (38) 
where, 


X(OF,/ Om) + ¥(OY:/ OG m) + Z(O24/ OG»). 
(39) 


( OF ;/ Of m) 


Explicit formulas for the elements of — the 
kinetic-energy matrix may now be obtained from 
eq 33 and 34, combined with eq 36 to 39. The 





The elements for rotation are the moments and 
products of inertia of the molecule about axes 
parallel to the s,y,2-axes, but passing through 
the center of mass of the molecule, 


my? + 2?)— M(yi- 
: 2m (2 +2?)—M (23 

2m (2? t+-y?)—M (xi, 

2m iyi — Mruyu 


Dom iy i2i— My 





- 2om,2.01— Mzytu 


The total mass of the molecule has been denoted 
by M, and the coordinates of the center of mass 
of the molecule, for the equilibrium configuration, 
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reductions of the scalar products of the s-vec ‘ors 
are made by the standard methods of vector 
analysis. 
then be applied explicitly to the translations. — [y 
this reduction, S,, in eq 4 and 16 represents th 
submatrix for translation, and is a 3X3 diag. 
onal matrix, with each diagonal 
equal to the total mass of the molecule, M 
Sx» represents the submatrix for rotations and 
coordinates, and Sj» the submatriy 
for the interaction of these coordinates with 
translation. The elements of Sy» are simply 
the z-, y-, and z-components of linear momentun 
of the rotations (about 
internal coordinates. 

The reduced matrix for the rotations and in- 
ternal coordinates, corresponding to Sj. in eq It 
which is obtained in this way may be written, 


element 


internal 


the x, y, 2-axes) and 


by ry, Y¥w, and 2. The first term in each equa- 
tion of 41 is the moment or product of inertia 


about the arbitrary origin; the second term arises ” 


from the reduction process of eq 16, and vanishes 
if the origin coincides with the (equilibrium 
position of the center of mass. 


The elements for the interaction of rotation and ~ 


internal coordinates may be written 
H,.,.=XHy, Hym=yHe, Hem 

where H,, is the vector 

HH, = Dome, X (Or,/ OY m) — ae X 2am s(OF:/ Oym) (43 
7 


The vector from the origin of the axes to the center 
of mass of the molecule has been represented by 


Iu, 


Py=M Som t= Xt +Yyu t+ 2eu. (44 
' 
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the origin of the z,y,2-axes coincides with the 

guilibrium) position of the center of mass of the 
molecule, the vector ry, and the second term in 

, 43, vanish. The second term arises from the 

duction process of eq 16. The first term in eq 
1} represents the vector angular momentum of the 
motion 9, about the origin, while H,, represents 
the vector angular momentum of q,, about the 
equilibrium) center of mass. The quantities 
Hom, Lym, And Hy», are therefore the components of 
angular momentum of q,, about axes through the 
center of mass parallel to the z, y, 2-axes. 

For internal coordinates, 


K nm=D>lm (Or ;/ OG m)?—M [2pm (or. Ogn) | (45) 


and for the interaction of two internal coordinates, 


K on = Km > lM (OF ;/ OF m)+ (OF ;/ Oy) 


M [2am (or, Ogm) |-{ 2ami(Or, oq, | (46) 
1G i 


The second term in eq 45 may also be written 
\/(Ory /Oq,,)?, and therefore vanishes if the position 
of the center of mass of the molecule is independent 
of dm. Similarly, the second term in eq 46 is 
M (Or y/OGm)*(OF/O,), and vanishes if the center of 
mass is independent of either g, or qg, (or both). 
However, the vanishing of the (equilibrium) 
vector fy is not, of itself, sufficient to cause these 
terms to vanish. The second terms in eq 45 and 
1} represent the kinetic energy of translation of 
the center of mass of the molecule, and arise in 
the reduction process of eq 16. The first terms in 
these equations represent the kinetic energy of 
the internal coordinates computed with respect 
to the arbitrary origin. Subtraction of the second 
term from the first yields A,,,, and K,,., which 
represent the kinetic energy of the internal co- 
dinates computed with respect to the instan- 
taneous position of the center of mass. 

The set of atoms associated with the internal 
‘ordinate gq, Will be defined to include those 
itoms which move with respect to the z,y,2-axes 
vhen the coordinate g, varies. This set is not to 
regarded as fundamental, because the atoms 
neluded in it depend, in general, on the definition 
f the axes. The vector (dr,/0q,,) vanishes by 
lefinition for all atoms not in the set associated 
vith qm. As a result, the summations in eq 43 

| 45 need extend only over the set of atoms asso- 
ited with g,. In eq 46 the two summations in 


requencies of Semirigid Molecules 


brackets extend over the sets of atoms associated 
with ¢g, and q,, respectively, whereas the first 
summation in eq 46 extends only over those atoms 
included in both of these sets. 

The final reduced kinetic-energy matrix is ob- 
tained from the matrix 40 by applying the reduc- 
tion process of eq 16 to the rotations, and to any 
internal coordinates that do not appear in the 
potential-energy expression. The elements of the 
matrix of eq 40 are not independent, in general, 
of the choice of the z,y,2-axes, but the reduced 
matrix is independent of the axes after the over- 
all rotations have been removed. The reduced 
matrix is most easily obtained numerically, by 
applying the reduction process of eq 10.to one 
row and column of the matrix at a time (so that 
in each step S,, in eq 16 is a matrix with a single 
element). This procedure makes it unnecessary 
to reciprocate directly the submatrix for the co- 
ordinates that are to be removed. 


(b) Equations for semirigid molecules 


The equations for the calculation of the kinetic- 
energy matrix for generalized internal coordinates 
have been given in the preceding section. These 
equations may now be applied to the two types 
of internal coordinates which will be found use- 
ful in specifying the configuration of a semirigid 
molecule. These are coordinates that represent 
either a translation or a rotation of their associated 
sets of atoms with respect to the z,7,2-axes.  In- 
ternal translational and rotational coordinates will 
be represented by q, and q,, respectively, where 
q, is the distance of translation, and q, is the angle 
of rotation, both measured from the equilibrium 
configuration. 

Let u be a unit vector in the direction of transla- 
tion for the coordinate q,, with components 
Mr, My, and uw,, and let = be a unit vector directed 
along the axis of rotation for the coordinate q,, 
with components £,, &,, and &,._ Let p, be the posi- 
tion vector of the ith atom associated with 4q,, 
with respect to an arbitrary point on the axis of 
rotation for g,. In terms of these vectors eq 38 
becomes, 

(Or,/Oq,) =H, (47) 

and 
(Or,/Oy,) =EX Py, (48) 
for translational and rotational coordinates, re- 
spectively. Equations 47 and 48 may now be 
substituted in eq 42 to 46, to obtain the elements 
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of the kinetic-energy matrix, eq 40, in vector form. 
These vector expressions may then be reduced by 
the usual methods to more useful scalar expressions 
involving the masses, moments of mass, and 
moments of inertia of the associated sets of atoms. 
The details of these reductions will be omitted and 
only the final formulas stated. 

Some additional definitions will be necessary. 
The total mass of the set of atoms associated with 
a coordinate g, (m=t or r) will be denoted by 
M,,, and the coordinates of the center of mass of 
the set of atoms by Sym, Yarm, ANd Zym. The 
corresponding quantities for the entire molecule 
are M, and ry, yy, and 2y. 

For each rotational coordinate, g,, define a set 
of orthogonal right-handed £, n, f-axes, with the 
origin at a point 7,,, Yor, Zr, on the axis of rotation, 
but otherwise aribtrary, and the £-axis coinciding 
with the axis of rotation. The direction cosines 
of the &, n, ¢-axes with respect to the 7, y, 2-axes, 
written in matrix form, are 


nN: 


The determinant of this matrix will be +1 if 
both the x, y, 2-, and &, », ¢-axes are right-handed. 

The following moments and products of inertia 
of the set of atoms associated with gy, will be 
required, 


¢?)) 


I, is the moment of inertia about the axis of rota- 
tion. The products of inertia, J;, and J;:, vanish 
if (1) the &-axis is a twofold, or higher, axis of 
symmetry for the associated set of atoms; (2) the 
n, ¢-plane is a plane of symmetry; or (3) if a three- 
fold, or higher, axis of symmetry lies in the », 
¢-plane, and intersects the §-axis. J;, (but not 
necessarily J;;) vanishes if (1) the ¢,&-plane is a 
plane of symmetry or (2) the n-axis is a twofold 
axis of symmetry. 

Let d, be the perpendicular distance from the 
center of mass of the set of atoms associated with 
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q, to the axis of rotation, or £-axis, and let \,, \, 
and », be the direction cosines, with respect to 
the z,y,2-axes, if the line along which this cen‘e 
of mass moves when q- increases from zero to 4 
small positive value. 

That is, 


and 


(nurs — Saernz)/d, 
(nserky—Saurny)/d, 


(narrk: “* Surn:)/d, 


where, 


nur M Som ns, ure M Sm ss. 
! i 


The sums in eq 50 and 53 extend over the atoms 
associated with g,, Atoms on the axis of rotation 
are not with q,, according to th 
definition given in the preceding section. It may 
be more convenient to calculate these quantities 
Define a,, a,, and a, to be equal to 


associated 


as follows. 
their respective cofactors in the determinant 


= (a,/d,), 4,=(a,/d,), d= (a,/d,). (56 

The equations for the elements of the kinetic- 
energy matrix, eq 40, may now be stated. 
translational coordinate, q,, H7,,, H.,, and H,, ar 
equal to their respective cofactors in the deter- 


minant 
es 
Mu: 


I, 
Mu, 


H,, 
Mu, 


(Iue—Tar) (Yare— Ya) (23¢1— Zar) 
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e translation arises from the stretching of a 
u,. My, and uw, are the direction cosines of 


yond. 


hor a rotational coordinate, q,, 


H,=H2+H 4 
H,,=H? fats 


H,=H®+H® 
vhere 1), H2, and H'? are defined by the equations 
A? =é-l; ge 
HY? =é,1; oT es 


v 


H ; S24 z4§ t: ; 


nel em 


nyl ey 


ind 12, HW, and IT? are equal to their respective 
~ofactors in the determinant 
H? H? i? 


M.d,2, 


(60) 


Ma», Man, 


' Iv ‘Yor Unt \Zor <u 
7..=0, and eq 59 reduces to 


EJ, WY =e. (61 


Bn many cases, J 


iH’, 


td:, HY 


Before stating the formulas for the remaining 
lements, A,,,, and K,,,, of eq 40, it will be neces- 


“Bary to differentiate between connected and un- 


onnected coordinates. If the sets of atoms asso- 
lated with two coordinates q,, and q, contain no 
toms in common, the two coordinates will be 
aid to be unconnected 
ected. Whether two coordinates are connected 
r unconnected will be dependent, in general, on 
The 
Bnly types of connected coordinates that will be 
Monsidered are those for which the set of atoms 


otherwise they are con- 
Bhe way in which the r,y,2-axes are defined. 


rr one coordinate, say ¢,, is completely contained 
ithin the set of atoms for the other coordinate, 
A special case is that in which the sets of 
Stoms for q,, and q, are identical. If q,, and q, 
Bre unconnected, this fact will be indicated ex- 
licitly by writing K,,,(—K,,) as K,.,.(=—K, »). 
iy, and q, are connected, and the set of atoms 
” 4, is completely contained within the set of 
Bioms for qg,, then K,,,(—K,,,) will be written 
ik Kinym). If the set of atoms for q,, and 
are identical, A,,..,— Ay... Subseripts ¢ and r 
ili again be used to indicate translational and 


requencies of Semirigid Molecules 





rotational coordinates, and a second coordinate 
of a given type will be distinguished by the use of 
primes. The equations for the several possible 
forms of K,,,, and K,,,, may now be written as 

(62) 


K,.=M,A—M"M,) 


K, MMM! (pew’ 4 MyM (63) 
Kiuey=M/(1 


K, MM ,M,A,, 


M ‘Ml, (un’, + (64) 
(65) 
(66) 


Kiin=M,(1—M"M)) An 


(67) 


kK, M Arne 
M"' Mid? 


MM, AL,A,, 


K,.=I, (68) 


K, M~OM_M’ dd’ ,.4.N +, ,+.';) (69) 
Kye) = (ET + &, + £,HO") + Mid) Ay 


M~'M,Mid.di(\,d,+- dA, + ALND) 


(70) 


The symbols A,,, Ax», and A,,), represent the 


following determinants: 


Air d,(u;r,z T My Ay T bzX,) 


In eq 70, 71)", 2", and H®’, are the quantities 
calculated from eq 59 (with r=r’). When eq 61 
applies, the first term in parentheses in eq 70 
reduces to, 


(74) 


If the center of mass of the set of atoms asso- 
ciated with a rotation, g,, lies on the axis of rota- 


tion, so that d,=0, g, will be called a balanecd 
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rotation. For a balanced rotation, 7® =H®= 
H®=0, in eq 58, and there are obvious simplifica- 
tions in eq 65 to 70, as A,,=0 (in addition to 
d,=0). A balanced rotation may be called acci- 
dentally balanced if the £-axis is not an axis of 
symmetry of the associated set of atoms. Such 
cases are very unlikely to occur. As a rule, the 
f-axis for a balanced rotation will be a twofold or 
higher axis of symmetry. In this case, /;,= /¢;=0, 
and eq 59 is replaced by eq 61. However, it is not 
necessary that a rotation be balanced for /;, and 
/¢; to vanish. 


Ill. Application to Propane, Toluene, and 
Ethylbenzene 


1. Kinetic-Energy Matrices 


The application of this method to certain vibra- 
tions in ethylbenzene will be described in this 
section. In order to obtain the necessary force 
constants for ethylbenzene, it also has been neces- 
sary to set up the kinetic-energy matrices for 
propane and toluene. In the semirigid model of 
ethylbenzene for which the calculation has been 
carried out, it is assumed that the benzene ring 
and the attached hydrogens move as a rigid body, 
and also that the methylene and methyl groups 
move as rigid bodies. That is, the force constants 
for all deformations of the benzene ring and 
attached hydrogens and for deformation of the 
H—C—H angles in the methylene and methyl 
groups have been set equal to infinity. The force 
constants for the stretching of the two C—C bonds 
in the ethyl group have also been assumed infinite.’ 
In order that the force constants obtained from the 
observed frequencies of propane and _ toluene 
should be on a comparable basis, corresponding 
assumptions have been made for these molecules. 
The coordinates remaining in the problem are 
therefore the rocking and twisting motions of 
these rigid groups. 

The total number of internal coordinators in 
ethylbenzene is 48; of these, 8 are retained in the 
present calculation, and 2 are assumed to have 
zero frequency. In calculating the kinetic-energy 
matrix by the present method it is necessary to con- 
sider only these 10 internal coordinates, plus the 


? As no bond stretchings are included in the present calculations, they will 
not illustrate the use of the formulas for translational internal coordinates 
given in section II, 2b. 
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3 over-all rotations, and then to remove the * ¢o. 
ordinates of zero frequency, as discussed jp 
section II,1. In order to obtain the reciprocal 
kinetic-energy matrix for the 8 internal coordi. 
nates, it would be necessary to set the matrix up 
for the 46 internal coordinates of nonzero fre- 
quency, and then to remove 38 coordinates, as 








discussed in section IT,1. 






Before enumerating the coordinates, the sym- 
metry properties of the molecules will be reviewed 
briefly. The propane molecule has the symmetry 
(,,, with the nondegenerate representations A 
A,, B,, and B,. Representations A and B ar 
symmetric and antisymmetric, respectively, to 
rotation by 180 degrees about the twofold axis 
Representations with subscripts 1 and 2 are syn- 
metric and antisymmetric, respectively, to reflec- 
tion in the symmetry plane. The toluene an 
ethylbenzene molecules have the symmetry (C), 
with nondegenerate representations A’ and A’ Pf 
which are symmetric and antisymmetric, respec- 
tively, to reflection in a plane perpendicular to th 
plane of the benzene ring.* 

















The notation adopted for the internal coordi-f 
nates is as follows. In-plane and out-of-plan 
(with respect to the symmetry plane) rocking 
motions are denoted by o; and x;, respectivel) 
Internal rotational motions are denoted by ¢ 
The subscript 7 indicates the group, as follows 
Propane, methyl (1 and 2), and methylene (3 
ethylbenzene, methyl (1), ethyl (2), and methylew 
(3); and toluene, methyl (1), and methyl grou 
(2).°. On this basis methylene twisting is repre 
sented by ¢. The C—C—C angle bending (v 
propane and ethylbenzene) is denoted by a. The 
coordinates represent in each case the angular 
deformation from the equilibrium configuration 
In addition the over-all rotations about th: 
r.y,z-axes will be denoted by @,, @y, and ¢, 
respectively. 





- gent 














kinetic: 






onnect 





repre 










* The equilibrium configuration assumed for the ethylbenzene molec 
is that in which the plane of the benzene ring is perpendicular to the plas 
of the C—C—C angle inthe ethyl group. The angleof rotation of the meth) 
groups in these molecules does not affect the quanvities entering into '¥ 
calculations. 

* The axes for methy! and methylene rocking pass through the carbon atom! 
of the CH; and CH: groups, respectively, so that these are motions of (® 
hydrogen atoms only. The axis for ethyl rocking passes through the phen! 
carbon atom, so that the entire C:H, group moves as a unit. The motion® 
toluene of the entire CH: group as a unit about an axis through the phen! 
carbon atom has been called “methyl group” rocking to distinguish it free 
methyl rocking. Methyl rocking may conveniently be denoted by Ci 
rocking, and methyl group rocking by —C Hs rocking. 
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the case of propane it is necessary to con- 
t the following symmetry coordinates: 


(1/y¥2)(o,+ o2) 
(1/¥2)(o- 
(1/y2)(m, 
(1/y¥2)(a,— 





’ (1/ ¥2) (di — 2) 


These transformations are orthogonal, with the 


natrix, 


(1/¥2) 


spec: ge rom eq 26 there are obtained the relations 


D th 


vps 


(at ,a*) = (1/2)[(o1,01) + (o2,62)| + (01,02) 


ordi: 


(o ,o }= (1/2)[(o,,0;) (o2,02)]— (0;,¢2) 


plat 

“king 
vely 
von (o*,q) 
lows 


+ a) = (1/2)[ (01,01) — (o2,02)|=90 


ge ore CT: 


(1/-y2)[(o1,q) + (o2,q)] 





(a~,q) = (1/-¥2)[(o1,q) — (o2,9)] 
vlet 


rroul 


ALR CIN PE Mio 


here (ot,ot) is the diagonal element in the 
sinetic-energy matrix for o*, (o*, ¢~) is the element 
io fa (equal to (o7,o*)), ete. 
Thy represents any coordinate other than 01, 9, a", 
ro. Similar relations hold for the x and @ co- 
ation The elements for the internal coordi- 
- thefmpates are first calculated from the equations given 
m section II, 2b, and the elements for the sym- 
etry coordinates are then calculated from eq 76. 


os 


‘epre- :, 
onnecting of and ¢ 


gular 
brdinates. 


d ¢; 


The notation used for the coordinates, their 
nokeofYMetries or representations, and their conven- 
he p“EBional descriptions are summarized in table 1. 
e meth ~ . ° ° ° 
into opmnteraction terms in the kinetic- and potential- 

nergy matrices between coordinates belonging to 
on atoms 7 . ° ° . 
ns ot ofailferent representations vanish identically. 

e pheny 
notion = 
e phet 
h it free 
by CE 
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Tare 1.—Notation and symmetry for coordinates 


Symmetry Coordinate ¢ Description 


PROPANE (Ca) 


C Hy, rocking 
cC—C—C bending 


Over-all rotation 
CH; rotation 


CH: twisting 
C Hs rocking 


Over-all rotation 


C Hy rocking 
CH, rocking 


Over-all rotation 
CH; rotation 


CH, rocking 
CH: rocking 


LBENZENE (Cy) 


Over-all rotation 


CH; rocking 
CoHs rocking 
C Hi: rocking 
C—C—C bending 


Over-all rotation 
Over-all rotation 
C Hy rotation 
‘»Hs rotation 


‘H2 twisting 
“Hs rocking 
‘2Hs rocking 
*H2 rocking 


TOLUENE (Cy) 


Over-all rotation 


C Hy rocking 
—CHs3; rocking * 


Over-all rotation 
Over-all rotation 
CH; rotation 


C Hy rocking 
—C Hs; rocking * 


* The coordinates above the dashed line in each representation were 
removed from the kinetic-energy matrices (see text) 








The dimensions used in the calculations are as 
follows. The benzene ring was assumed to be a 


plane hexagon, with the length of the C—C bonds 


1.39 angstrom units. The other bond lengths, in 
angstrom units, were 1.54 for the C—C bonds in 
propane and in the alkyl groups (including the 
bond joining the alkyl group to the ring [8]), 1.09 
for the C—H bonds in propane and the alkyl 
groups, and 1.08 for the C—H_ bonds on the 
benzene ring. Tetrahedral angles were assumed 
in propane and in the alkyl groups. In every case 
the numerical values of the elements of the 
kinetic-energy matrices (tables 2 and 3) have the 
units of mass in atomic-weight units multiplied by 
the square of a length in angstrom units. 


TABLE 2.—Kinetic-energy matrices for propane before 


removal of over-all and internal rotations 


1. 57 
0. 026 


3.16 
0 
0.045 


2. 358 





Tar.e 3.—Kinetic-energy matrices for propane, tolue 
ethylbenzene 


Propane (C >.) 


hit ir 
origu 
Th 
treat 
hisect 
1. 65 Oo. 084 
0. OS4 724 
kinet) 


and ¢ 


out 


The s 
the 

latter 
energ' 


nonor 





Equat 


Ethylbenzene (C 


where 

The 
of the 
matrix 


0. O89 
1. 66 


es Nr ST 
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. stated previously, the intermediate steps in 
calculation are dependent upon the choice of 
-axes, but the final kinetic-energy matrices 
er removal of over-all rotation) are not. The 
ulations for propane are described in some 
tail in order to illustrate the method, and the 
vice of axes for propane will therefore be stated. 
The z- and 2-axes were taken in the plane of the 
(°-C—C angle, with the origin at the vertex of 
the angle, and the z-axis bisecting the angle. The 
axis therefore remains a twofold axis of sym- 
metry as the C—C—C angle is deformed. The 
origin is not at the center of mass of the molecule. 
The coordinate @ in propane is most easily 
treated as follows. By definition the c-axis 
hisects the C—C—C angle. Denote the changes 
in these two half-angles by a and a. The 
kinetic-energy matrix is first set up for a and a, 
and then the following transformation is carried 
out 


a= (a, + Ay) 


0 (ay ay), 


The second equation is merely the condition that 
the z-axis bisect the C—C—C angle when the 
latter is deformed. The elements in the kinetic- 
energy matrix for (a@,—a,) are then ignored. The 
nonorthogonal matrix of this transformation is 


l l 
l —j 


Equation 26 leads to the relations 


(a,a) = (1/4) [ (ay ,a,) +2(ay,a2) + (a2,a2)} (78) 


(a2,9)], 


aq) = (1/2)[(a,q) 4 


where g is any coordinate other than a, a, or a. 
The following remarks apply to the calculation 
of the elements of the unreduced kinetic-energy 
If the origin of the 
£, », f-axes (the point 2,,, Yor, Zor) is located at the 
carbon atom on the axis (or at any point on the 
&axis for the balanced rotations ¢;, ¢, and @,), 
then J;, and J,:, vanish in each case, and 7‘), 
and H%? may be calculated from eq 61. 

H?, HS, and H vanish for the balanced rota- 


matrix, eq 40, for propane. 
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The diagonal terms K,, are 
For balanced rotations 
subscripts 


tions ¢), @:, and 3. 
obtained from eq 68. 
K,,=/;. Coordinates with different 
(1, 2, or 3) are unconnected, and K,, is calculated 
from eq 69. A,, vanishes if either rotation is 
balanced (¢-coordinates), or if the directions of 
motion of the centers of mass are perpendicular 
(thus r-coordinates are perpendicular to o- and 
Coordinates with the same sub- 
is obtained from 


a-coordinates). 
scripts are connected, and A, 
eq 70 (the first term being given by eq 74). 
From a consideration of the directions of the 
taxes, and the directions of motion of the centers 
of mass, for each coordinate, it is readily seen that 
the only terms that do not vanish are Koyo, 
(¢=1,2). 

The numerical values of the elements of the 
matrices of eq 40 (one matrix for each representa- 
tion) are given for propane in table 2. 

The next step is the removal of the over-all 
rotations by the reduction process of eq 16. In 
the subsequent calculations the frequencies of the 
¢ and @» internal coordinates (methyl rotations in 
propane and toluene, and methyl and ethyl rota- 
tions in ethylbenzene) have been assumed to be 
zero. This is justifiable because the potential 
barriers for these rotations’® correspond to suffi- 
ciently low frequencies for small oscillations from 
the equilibrium configuration. 
@, and @ have, therefore, also been removed by 
the reduction process of eq 16. 

The final reduced kinetic-energy matrices, one 
for each representation, for each molecule, are 
shown in table 3. 


The coordinates 


2. Calculation of Frequencies 


The force constants are now to be evaluated 
from the observed frequencies of propane and 
toluene. These force constants may then be sub- 
stituted in the secular equation for ethylbenzene 
to yield calculated frequencies for the latter 
molecule. 

The frequency assignment for propane has been 
taken from Pitzer [9], and the assignment for 
toluene from Pitzer and Scott [10], The fre- 
quencies (v) are summarized in table 4. The 


#” Propane, reference [9]; toluene and ethylbenzene, reference {11} 








normal coordinates are mixtures of the coordinates 
in each representation, and when this mixing is 
appreciable it is not possible to assign the fre- 
quencies in each representation to the coordinates 
without some ambiguity. However, the pairing 
indicated is fairly certain. Previous investigators 
have assigned the two 73, frequencies as ~ (748) 
and x; (1179), but the assignment of these fre- 
quencies in table 4 is supported by (1) the force 
constants obtained for @», ¢3, and 2, which should 
have approximately the proportions 2,1, and 1, 
respectively (these figures assume that interaction 
constants are negligible, and that the main con- 
stants arise from the constants for the four 
C—C—H angles), and (2) the approximate 
equality of the force constants obtained for 2 
and o. 


TaBLe 4.—Frequencies and calculated force constants for 
propane and toluene 


PROPANE 


Sym- 


) 
Description metry 


Coordinate 


C Hy rocking 
C-C-C bending 
C Hy, twisting 

C Hy rocking 

C Hy rocking 
CH, rocking 

C Hy rocking 

C Hy rocking 


TOLUENE 


C Hy rocking 
—C Hy; rocking * 
C Hy, rocking 
~C H, rocking * 


The next column of table 4 gives the value of 
d (eq 22 to 24) in the units 


A\=[v(em~')/1000)°. 


The last column shows the calculated values for 
the force constants, *, as calculated from the 
secular equations. They are the diagonal ele- 
ments in the potential energy, or F, matrix 
(off-diagonal elements have been set equal to 
zero as there are not sufficient data to determine 
them). The units of these force constants corres- 


pond to the units that have been used for ) ay¢ 
for the kinetic-energy matrices in table 3; in, th 
reverse calculation for ethylbenzene the frec wep. 
cies are again obtained in em“. 

Table 5 shows the force constants selected { 
the ethylbenzene molecule. Those derived froy 
toluene are o) from o, and m from m. Thos 
derived from propane are o; from (o*+o°)2 
m, from (r++ 2~)/2, and o3,73, ¢; and a@ from «, 
™;, 3, and a, respectively. The averaged fore 
constants, (o* + o7)/2 and (x*+ 2~)/2, are the mai 
force constants for methyl rocking in propane 
The observed separation of the frequencies fq 
o*,o, and for x*, x~, in propane, is due to methiy!- 
methyl interaction constants and the twofol 
symmetry. 

TasLe 5.—Force constants and calculated frequencies fo 
ethylbenzene 


ETHYLBENZENE 


Sym- 


Description metry 


Coordinate 


C Hy; rocking 1. 067 

“:Hy rocking d 0. 0187 

“Hy, rocking 4 : 1.721 1, 310 
*-—C—C bending 0. 1437 Rt) 
‘He twisting 1.615 1, 270 
‘Hs rocking 1.144 1, 070 
“sH, rocking 0. 1025 320 
*H> rocking A” 5 . 004 775 


Substitution of these force constants in the see- 
ular equations for ethylbenzene then yields th 
values of \, and the calculated frequencies, show: 
in the last two columns of table 5. The observed 
frequencies for ethylbenzene are discussed in thy 
next section. 


IV. Vibrational Frequencies of 
Ethylbenzene 


The present calculations were undertaken as a! 
aid to the assignment of the vibrational frequencies 
of ethylbenzene.* These frequencies were require 
for the calculation of the thermodynamic function: 
by statistical methods [11]. The frequencies hav 
been assigned on a semiempirical basis whic! 
involved a detailed consideration of only th 


*F. G. Brickwedde, M. Moskow, and R. B. Scott [24] have also present 
an assignment of the vibrational frequencies of ethylbenzene. 
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st frequencies of vibration of the benzene ring 
also of the frequencies associated with the 
ation of the ethyl group. The remaining 
iencies were taken from the assignment for 
toluene [10]. A complete frequency assignment 
for ethylbenzene at this time is impracticable. 
several workers have reported Raman spectra 
13, 14, 15, 16, 17, 18] and infrared spectra [19, 
2), 21, 22, 23) for ethylbenzene. The 
observed frequencies in the Raman spectrum of 


ethylbenzene, are, approximately, 154, 305 (weak) 


399 (weak), 452, 485, 538 (weak), 5: 
and 620 em. The frequencies 480, 567, 


lowest 


> 
) 
53 (weak), 


and 
617 em7! have been reported in the infrared spec- 
trum in this region. It will be helpful to list 
the six lowest frequencies of toluene, as assigned 
by Pitzer and Scott [10]: 11, 216 (R); 18b, 340 
R); l6a, 405 (R); 16b, 467 (R); 6a, 521 (R); and 
6b, 622 (R, IR). The first number designates 
the mode of vibration (see fig. 6, of Pitzer and 
Scott), the second is the frequency in em™', and 
R and IR refer to the presence of the line in the 
Raman or the infrared spectrum, respectively. 
The vibrations 11 and 1I8b correspond princi- 
pally to the recking of the alkyl group perpendic- 
ular to and parallel to the plane of the benzene 
ring, respectively. The frequencies of these modes 
in ethylbenzene were calculated as 137 cm™' for 
11 and 320 em~' for 18b by means of the approxi- 
mate normal coordinate analysis of section III 
vibrations denoted 
respectively, in that section). The only possibility 


these were by o. and m, 
for the assignment of the 154 em™! frequency is 
therefore mode 11. The mode 18b was assigned 
the calculated frequency of 320 em™'. The 
frequeney observed at 305 em™' may easily be 
the overtone of 154 em”, 

The mode I6a has a frequency of 400 em~! in 
benzene and 405 em~' in toluene. The observed 
frequeney at 399 em in ethylbenzene may 
represent I6a, or an overlap of l6a with the 
carbon skeletal bending in the ethyl group, which 
should fall at about 380 cm~'. The toluene value, 
405 em™',was retained for ethylbenzene. The 
mode 16b has the frequency 467 cm~' in toluene; 
the nearest observed frequency in the spectrum 


of ethylbenzene, 452 cm~', was assigned to this 
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mode. The mode 6a in toluene gives rise to a rel- 
atively intense line at 521 cm~'. The variation of 
the frequency of the mode 6a in other monosub- 
stituted benzenes, as a function of the mass of the 
attached group [12], indicates that 6a should be 
lowered somewhat in ethylbenzene, since the 
effective mass of the ethyl group in this mode is 
considerably greater than that of a methyl group. 
Therefore, the frequency 485 em™~', which is the 
strongest observed line in this region, was assigned 
to 6a. The mode 6b has a frequency of 606 em~! in 
benzene, 622 cm~' in toluene, and remains essen- 
tially unchanged in frequency in other monosub- 
stituted benzene !12]. The frequency observed in 
ethylbenzene at 620 em~! was therefore assigned 
to 6b, 

The mode 7a, which is principally the stretch- 
-C bond joining the alkyl group to 
‘in toluene. 
in 


ing of the C 
the ring, has a frequency of 1,210 em 
As a strong line appears at about 1,200 em 
the ethylbenzene Raman spectrum, the toluene 
frequency was retained for 7a. The other C—C 
stretching mode in the ethyl group was assigned 
the frequency of the strong line at 965 em™ in the 
Raman spectrum of ethylbenzene as this is close 
to the average of the two carbon-carbon stretch- 
ing frequencies of 868 and 1,053 em~' in propane 
{9}. The symmetry that splits these two com- 
ponents in propane is absent in ethylbenzene. 
The symmetrical bending in the 
methyl group was assigned the frequency 1,385 
em™~', and the doubly degenerate unsymmetrical 
bending the frequency 1,460 em™'. The symmet- 
rical hydrogen bending in the methylene group 
was assigned the value 1,460 em~'. The value 
2,960 cm~' was assigned the C—H_ stretchings in 


hydrogen 


the methyl and methylene groups. These assign- 
ments were made on the basis of the frequencies of 
these modes in propane and toluene, and to some 
extent on the basis of the spectrum of ethylben- 
zene. 

The remaining frequencies of vibration of the 
ethyl group were assigned values calculated by 
means of the approximate normal coordinate 
analysis of section III (see tables 1 and 5). 

The complete frequency assignment for ethyl- 
benzene is shown in table 6. 








TABLE 


6.— Vibration frequencies of the ethylbenzene molecule 


Sym- Fre- 


‘ype of Vit 
Type of Vibration metry quency 


Vibrations principally characteristic of the benzene ring or pheny!] radical 


n¢c 
ISb (CoH, rocking 


lfia 
lfb 
fa 
fib 
‘ 
lob 
12 
la 
7b 
7a 
l 
Isa 
15 
vb 
Ye 
7a (¢ 
3 
vb 
1a 
sb 
sa 
4 
13 
7b 


an 
Wb 


‘sHs-CHeC Hy, angle bending 


“Hy 
“Hy 
“Hy 
“Hy 
“Hy 
“Hy 
"Hy 
"Hy 


‘H 
“Hy 
“Hy 


2 [~ stretching 


-C Hy stretching 


»Hs rocking 


“sHy-CyHs stretehing) 


Cn eet ll cel ee ct ce eel ee el cel el eel cl eel el oe ee 


pm D> De me 


Vibrations principally characteristic of the ethyl group 


rocking 


rocking 
rocking 
twisting 
rocking 
symmetrical bending 
symmetrical bending 


unsymmetrical bending 


symmetrical stretching 
unsymmetrical stretching 
symmetrical stretching 


nn nn 


unsymmetrical stretching 
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Appendix 1 " 


Let the coordinates q, be divided into two groups, 
and q, and partition S,, into the corresponding subme- 
trices. Then the kinetic energy matrix, eq 4, becomes 

S 


aa 


Sab Sa2 
S. b S bb S 2 
Ss Be Sx» 


nl ° ° . ' 
The removal of the cross terms in the kinetic energy ths 
contain gq leads to the reduced matrix 


(So — 84.55) (Sip — Sa8.z'S.) 

(S,, 8.8. Sa) (Syo . a2 'S,2) 
and the removal of the remaining cross terms that conta 
g» leads to the final reduced matrix for the coordinates 


i See section II, 1 (following eq. 16) 
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ition 81 may be rewritten in the form 


(3.=°+85'3, 0.4%. 


ad 


6.885" 


a SON WE Sis 
CG. a5 


Gy , ., ‘ 4 : ' . { 83) 


On the other hand, the simultaneous removal of the cross 
terms that involve ga and q» leads to the reduced matrix 


3 $ . . 1 . 
[Sa2 Soe ] Saa Sar Sw | 


Sat Spo S, | 


(S84) 


rhe matrices S’» calculated from eq 82and 84 are identical, 
provided that 
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where E,, and Ey, are unit matrices. But the last rela- 
tion is easily verified by multiplying out the matrices on 
the left. This establishes the theorem when q, is divided 
into two groups; the truth of the theorem for any number 
of groups follows by induction. 


Appendix 2 ™ 


From the relation that SG=E, where E is the unit 
matrix, and eq 3, there are obtained the relations 


SuGy + SnGr 4 
SiGy+ SpGu 4 


SieGig + SoGos 4 


S11Gy3 + Sy2Go + SiG (89) 
Now multiply eq 89 from the left by S.Si', and from the 
right by Ga' G,,, and add the resulting equation to eq 
86; next multiply eq 87 from the left by SjwS,;' and eq 88 
from the right by Gss'G,,, and subtract both of the result 
ing equations from the previous sum. The resulting 
equation may be factored as follows: 
(Ses — S955 °S1s) (Gos — GesGea Ges) = Es. (90) 


22 1 


Equation 21 follows immediately from eq 16, 20, and 90. 
WASHINGTON, September 24, 1946. 


2 See section II, 1 (eq 20 and 21 
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Analysis of a Standard Sample of Natural Gas by 
Laboratories Cooperating with the American Society 
for Testing Materials 


By Martin Shepherd *” 


This is a report of the analysis of a standard sample of natural gas by 30 laboratories 
cooperating with Subcommittee VII of Committee D—3 of the American Society for Testing 
Materials. The data are presented in a series of frequency distribution plots that show at 
a glance how the analyses from all laboratories compare with respect to each component 
determined, as well as calculated heating value and specific gravity. The heating value 
and specific gravity determined by analysis are compared with values carefully measured 
at the National Bureau of Standards. The analyses were performed volumetrically by the 


absorption and combustion methods, and the plots form a clear picture of this type of gas 





analysis in this country. 


standardization is evident. 


I. Introduction 


Subcommittee VIL of Committee D-—3 of the 
American Society for Testing Materials has been 
issigned the task of standardizing the analysis of 
raseous fuels. This is the second of a series of 
eports concerned with work preliminary to actual 
tandardization. In the first of these reports? it 
vas stated that the subcommittee was approach- 
ing their assigned task with certain reservations 
Mderived from an acquaintance with gas analysis 
methods, gas analysis apparatus, and gas analysts. 
rhe reason for such reservations usually becomes 
Wapparent when a gas analyst is asked to sub- 
Blilute a new reagent, pipette, procedure, method 
a Whole new apparatus or system of analysis 
or one that he has been using satisfactorily. 
uch a request is unamiable and might very often 


Chairman of Subcommittee VIL (Analysis of Gaseous Fuels) of Com- 
tee D-3 (Gaseous Fuels), American Society for Testing Materials 
tin Shepherd, Analysis of a standard sample of the carburetted 
8 type by laboratories cooperating with the American Society for 
Materials, J. Research NBS 36, 313 (1946) RP1704 
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Although some very creditable work is reported, the need for 


be unreasonable. Indeed, if it can be demon- 
strated by actual measurement that an apparatus 
or method now in use will give satisfactory results, 
there is no real justification for insisting that it be 
discarded in favor of another that may be spon- 
sored as the official instrument of some group. 

As it seems reasonable and fair to approve any 
apparatus and method capable of giving satis- 
factory analytical results, it remains only to decide 
what are satisfactory analytical results and which 
methods and apparatus will yield these results. 
To reach these decisions, two steps are being 
taken. 

First, the various purposes to be served by the 
analytical data are listed, and the necessary 
accuracy with which each component of each type 
of fuel gas must be known in order to serve each 
specific purpose is then estimated. These esti- 
mates afford the first criterion by means of which 
the suitability of analytical methods and apparatus 
may be judged, but they are subject to revision 
when more is known about the limiting attainable 
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accuracies of the analytical methods and also 
when the significance of the term “necessary” in 
the phrase “necessary accuracy” has been evalu- 
ated with more realism than enthusiasm. 
Second, the accuracy and reproducibility of 
existing apparatus and methods are being meas- 
ured by the direct procedure of conducting a 
series of cooperative analyses of standard gas 
types. (Although these 
samples are being prepared at the National 
Bureau of Standards, they are not to be confused 
with the Bureau's series of Standard Samples 
These gas samples are 


samples of various 


available for purchase. 
being prepared for free distribution to laboratories 
cooperating with ASTM D-—3% in this project). 

The cooperative analytical results obtained so 
far represent the first clear picture of the actual 
state of gas analysis in this country. 


II. Cooperating Laboratories 


These analyses are being conducted on a con- 
siderable scale. Cooperating laboratories include 
those of State and Federal agencies and of colleges, 
but the greatest contribution has come from the 
laboratories of the gas, petroleum, steel, and 
chemical industries. The laboratories are located 
throughout the United States. They are equipped 
with many types of apparatus, almost entirely 
modern, and employ different analytical methods. 
The type of apparatus used in the analysis of this 
sample is noted with the presentation of the 
analytical data. Of the 30 laboratories which 
analyzed the natural-gas sample, all but one (a 
college) employed chemists specifically trained in 
gas analysis. Taken as a whole, the work rep- 
resents the best present-day American practice. 


III. Preparation of Standard Gas Sample 
ASTM D3 VII-2 


Standard sample ASTM D-3-VII-2, natural 
gas type, was prepared to conform with the range 
of composition accepted for this test sample by 
Subcommittee VII: 


Percent hy 
rolu me 


CQ, 0 to 
Oy 0 to 
Ny Oto 5 
CH, 32 to 
CH, Oto 5 
CyH, 0 to 
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1,080 to 1,140 Btu/ft * 
0.630 to 0.750 


Heating value 
Specific gravity 


The sample was made up from the fol ow. 
ing compressed gases obtained from comme: cia! 


sources: 


Approximate 
purity 


Chief impurity 
Percent 

CO, 99. 9S Air 

No 99. 9 Oy 

CH,.. OS Ny, and C,H,4 

C,H, 9S C,H, and C,H, 

CyHg Cy Hy 


As this sample was not prepared from gases of 
known high purity, its actual composition wit) 


respect to all components is not known. It was | 


definitely established that no significant amount o! 
oxygen (i. e., less than 0.02 percent by volum 
occurred in the sample. In addition, its heating 
value and specific gravity were carefully measured 
The sample, therefore, has four definite purposes 
as it serves to measure (1) reproducibility, (2 
accuracy of the computed specific gravity, (3 
relation between the measured and the compute 
heating value, (4) ability to determine correct}; 
absence of oxygen in gas mixtures of this type. |: 
will not serve as a final criterion of accuracy (se 
footnote 3), since its exact composition is no 
known.* 

The sample was prepared under a pressure 0 
300 Ib/in.? in the large storage and mixing tanh 
used to prepare the mixtures for the Bureau: 
study of specific gravity instruments.’ After very 7 
thorough mechanical mixing, the uniformity of the | 


Lest there be any confusion as to the meaning of the terms accuracy a 
reproducibility, it should be stated that the conventional concepts are her 
implied. Accuracy is measured in terms of agreement with a true valu 
reproducibility is measured in terms of the mutual agreement of a series 
measurements of the same property of a substance. As a matter of 
venience, reproducibility has long served the gas analyst in lieu of a defir 
measurement of accuracy. Eventually this situation will be corrected 

‘Such terms as “known mixture” or “synthetic sample” are rather | 
lightly used in the literature. Unless the composition of a mixture is know 
to one order of magnitude better than may be expected from an analytic 
method, analysis of the mixture cannot he considered a measurement of | 
accuracy of the analytical method. If a mixture is “known,” the pur 
(exact composition) of each component must have been measured to a grea 
significance than the method used for its analysis will yield, and the porti 
ing of each component to form the mixture must have been done in suc! 
way that the relative amounts of all present are known again to a grea! 
significance than they can be resolved by analysis. To make such a mivtw 
usually seems too difficult, and the “synthetic sample” usually invoked 
not be truly definitive 

'F. A. Smith, J. H. Eiseman, and E. C. Cretiz, Tests of instruments for"! 
determination, indication, or recording of the specific gravities of gases, \! 
Miscellaneous Publication M177. (1947.) 
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Uber 


ectly 


sample was definitely established by sampling 
various parts of the tank and intercomparing 
portions by means of a 1—m Zeiss gas inter- 
neter to detect any difference in composition. 
interferometer is sensitive to changes of 0.01 
percent Or less of carbon dioxide and methane 
which were part of the standard sample. The 
suimple was then transferred to small sample 


S-ylinders as follows: 


|. The sample cylinders, equipped with needle 
valves having vacuum-tight packing, were con- 
nected to a single welded manifold by means of 
high-pressure unions. The large mixing tank and 
4 vacuum pump were connected to this same mani- 
fold. The manifold and opened sampling cylinders 
were evacuated to 0.0001 mm Hg pressure over- 
night. Thereafter, with vacuum pump discon- 
ected, the increase of pressure within this system 
vas less than 0.0005 mm during the first hour, 
ly a few minutes were required to transfer the 
was sample from the mixing tank to the sample 
ylinders. No significant contamination from 
Jeakage was possible under these conditions. 

2. The small cylinders were flushed with the 
sample from the mixing tank by alternate filling 
o 2 atmospheres and evacuating to 1 mm Hg. 
After the third flushing, the sample cylinders were 
filled to 300 Ib/in®. With full pressure on the 
manifold, the sample cylinders were then closed. 
The manifold was then opened to a water seal, 
and no leakage from the cylinders was observed 
wernight. 

3. All of the sample cylinders issued for the co- 
yperative analysis of this sample were thus filled 
imultaneously from the manifold connected to 
the single uniform source of supply. The samples 
were, therefore, identical at the date of issue. 

To avoid contamination while transferring the 
sample to the burette, the following illustrated 
instructions for a definite procedure were furnished 
the laboratories. 


INSTRUCTIONS REGARDING SAMPLER AND 
SAMPLING 


|. The eylinder containing the sample is shipped in a 
metal ease made from standard pipe fittings. Save this 
for return of the cylinder, which must be returned immedi- 
ate ipon completion of your analyses. 

lhe eylinder is equipped with a very special vacuum- 
po Nor close this valve with anything 
po Nov force the valve 


needle valve. 
the thumb and one finger. 
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shut. If too much twist is applied to the valve handle, 
the seat will be ruined. 

3. Remove the dust cap from the valve and attach the 
sampling fitting, which terminates in a thin brass tapered 
This sleeve fits a standard glass interchangeable 
joint, male part, size $ 10/30. This glass joint should be 
sealed to the brass sleeve by a thermoplastic cement of the 
DeKhotinsky type, or other cement capable of making a 
gas-tight joint. We find the arrangement sketched below 
desirable, 


sleeve. 














——About 50cm long 


Case for 
sampier 





Mercury seal 











EE ee 
Do not unscrew bottom cap 

Figure 1.—Flow diagram for sampling assembly. 
With cock of burette closed, gently open needle valve 
until gas bubbles out through mercury seal. When this 
line is flushed, take gas into burette and discard until 
horizontal portion of sampling line is flushed. About 50 
ml in all is adequate for flushing capillary lines. Leave 
connected to burette for next sample to be taken. Use 
mereury 
reduced pressure. 


seal as manometer to avoid sampling under 
If direct glass seal or interchangeable 
joint is not used in connecting to burette, flush sampling 
line each time with about 20 ml of gas. For convenience 
and safety, we prefer interchangeable joints secured with 
DeKhotinsky or a heavy lubricant. 

In no case did any analyst who made any of 
these determinations have any other knowledge 
of the properties of this sample, nor was any in- 
formation other than that given above imparted 
to anyone until all of the laboratories represented 
in this account had reported. During 1941, the 
specific -ravity was measured by E. C. Creitz and 
the heating values by J. H. Eiseman, both of 
the National Bureau of Standards. These values 
were turned over to the author, who, as chairman 
of Subcommittee VII, did not engage in the 
analyses conducted at this Bureau, nor did he 
consult these values until the Bureau's analyses, 
performed by Shuford Schuhmann, 
reported to him. 

In July 1944, R.S. Jessup, of this Bureau, re- 
determined the heating value. Also, during 1944, 
there were three determinations of composition 
made with the Bureau’s 
(manufactured by the Consolidated Engineering 


had been 


Mass spectrometer 
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Co.) under the charge of A. Keith Brewer. As 
the cooperative analyses had not then been re- 
ported, the first two measurements by the mass 
spectrometer could not have been made with a 
knowledge of the previously obtained results. 
The third analysis was made by Vernon H. Dibeler 
after all the other results from all sources had 
been reviewed by him, and proper corrections 
had been made for hydrogen. (See discussion 
in section V, 12.) 


IV. Apparatus and Methods Used for 
Analysis of Standard Gas Sample 


1. Apparatus 


It has been stated that most of the analytical 
apparatus used in this cooperative investigation 
was of the modern type. The apparatus may be 
classified with respect to the system of volume 
measurement, the arrangement of pipettes or re- 
action tubes, and the confining fluid. For con- 
venience, a set of abbreviations may be used. 

V,— Volumes are made comparable by means 
of a pressure-temperature compensator 
with manometer interposed between the 
compensating tube and the burette. 

Pressure within the burette is balanced 
against existing barometric pressure, and 
gas volumes are corrected from the ob- 
served pressure and temperature to a 
common basis, including a correction for 
changes in the saturation pressure of 
water. 

Pressure in burette is balanced against 

atmospheric pressure but no correction is 
made for changes in pressure or tem- 
perature during analysis. 

Volumes are measured by observing the 
pressure exerted within a constant vol- 
ume. 

Pipettes are connected to the burette by a 
manifold. 

A single pipette, connected to the burette, 
serves for all reactions. 

Pipettes are temporarily connected in suc- 
cession to the burette as the different 
reactions are progressively conducted. 

R,—The burette itself serves as a reaction tube. 

Hg—Mercury is used as the confining fluid. 

H,O (with appropriate subscript)—-An aqueous 
solution serves as the confining fluid. 


Thus the abbreviation V,R,Hg designates ay 
apparatus equipped with  pressure-tempers tun 
compensator and manometer, manifold conne: ting 
pipettes to burette, with mereury serving as th 
confining fluid; 23 of the 30 laboratories used this 
type of apparatus. Three laboratories used , 
V.R,Hg type of apparatus, two used a V,R,H.0 
type, another a Hemple apparatus (V,R,Hg 
another a V;R.H,O type, and finally one used , 
Bone and Wheeler apparatus, widely employed j 
Great Britain (V,R,H¢g). 


2. Methods of Analysis 


In general, the method of analysis used involve 
the conventional four steps, as follows: 


Method A: 
1. Absorption of carbon dioxide in an aqueous 
solution of potassium hydroxide. 

2. Absorption of oxygen in alkaline pyrogallo 
or in a chromous solution. 
single combustion over hot paltinum in th 
presence of excess oxygen, Mmeasuremen! 
of the contrection and carbon dioxid 
produced, and computation of all hydro- 
carbons present as CH, and C,H, using 
the conventional formulas: 


CH,=1/3 (4T7C—5CO.) 


C,H,= 1/3 (4CO,—27C), 
where 7( is the total contraction on burn 
ing, and CO, is the carbon dioxide produce 
by the combustion. 

Nitrogen by difference. 

Although 23 of the laboratories used this 

method, there were a few notable exceptions 

Four of the laboratories employed the fo- 

lowing method: 


Method B: 


1 and 2 as in method A. 

(One of the four laboratories then took fres 
portions of the sample for combustie: 
correcting the combustion data for t! 
CO, and O, originally found. The otle 
three used the general procedure of com 
bustion of a portion of the sample tak 
for absorption. ) 
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UeOls 


gallo 


in th 
men! 
OXid 


vdro- 
using 


burn: 


luce: 


\ single combustion over hot platinum in 
the presence of excess oxygen, measure- 
ment of the contraction, carbon dioxide 
produced, oxygen consumed, and_ the 

unreacted residue, and computation of 

all hydrocarbons as CH, and C,H, using 
the formulas: 


CH,=1/3 [7(7C+CO,) —90,] 


C.H,=1/3(60,—4(7C+CO,)! 


4. Nitrogen reported as the measured unre- 
acted residue. 

Method B has the advantage of eliminating 
rom the combustion data the additive errors 
used by the deviation of CO, from ideality, the 
oss of CO, by solution in water formed or origi- 
willy present, and the loss of CO, in any rubber 
‘onnections. All of these errors result in too high 

contraction and too low a produced carbon 
lioxide.” It will be shown that the few labora- 
ories using this method reported more consistent 
Whether or 
jot this consistency would persist if the method 


esults than the group as a whole. 


Bere used by all of the laboratories is yet to be 


Bietermined. 


In addition to the two chemical methods just 
One labo- 
atory employed one of the procedures ordinarily 
ised for the analysis of a fuel of the carburetted 
This procedure, method II de- 
previous report (see footnote 2), 


yutlined, there were a few departures. 


ater-gas type. 
cribed in a 
wluded absorption of unsaturated hydrocarbons 
mn fuming sulfuric acid and a fractional combus- 
ion over copper oxide, both preceding the com- 
ustion over hot platinum in the presence of 
xeess OXYgeN. 

In addition to these chemical methods, two 
iboratories used procedures involving distillation 
s follows: 


has been repeatedly demonstrated. See Martin Shepherd, An im- 


pparatus and method for the analysis of gas mixtures by combustion 

tion, BS J. Research 6, 121-167 (1931) RP266; Martin Shepherd 

bh R. Branham, Critical study of the determination of ethane by 

n over platinum in the presence of excess oxygen, BS J. Research 

RP625; Joseph R. Branham and Martin Shepherd, Critical 

e determination of ethane by explosion with oxygen or air, J 

\ BS 18, 377 (1934) RP715; Joseph R. Branham, Martin Shepherd, 

i Schuhmann, Critical study of the determination of carbon 

by combustion over platinum in the presence of excess oxygen, 
NBS 26, 571 (1941) RP1396 
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Method X: 

1. CO, and O, by 
sample. 

2. A single combustion of the hydrocarbons, 
with measurement of contraction and 
carbon dioxide produced. 

3. A separation of the hydrocarbons by dis- 
tillation in one of the usual rectifying 
columns. 

4. Correction of the combustion data to take 
account of hydrocarbons heavier than 
C,H,. 

Also, one laboratory employed only distillation 

for the analysis of the hydrocarbon fraction: 


absorption in a separate 


Method Y: 
1. CO, and O, separately by absorption. 
2. Other components by distillation. 


Finally, this sample was analyzed by the Con- 
solidated mass spectrometer, under the direction 
of A. Keith Brewer at the National Bureau of 
Standards. 

It will be shown that two of the laboratories 
employing contributed 
the maximum and minimum calculated heating 
respectively about 50 and 40 Btu on 
either side of the measured values. 


these special methods 


values 


V. Results of Cooperative Analysis 
of Standard Gas Sample 


1. Manner of Presenting the Analytical Data 


All the analytical data submitted have been 
tabulated, together with the average values de- 
rived from each series of analyses, but contempla- 
tion of these tabulated data for some hours would 
not serve to reveal what may be seen at a glance 
when the same data are presented in a series of 
frequency-distribution plots. Accordingly, these 
plots, which amount to actual pictures of the 
analytical data, are presented first, and the dis- 
The tabulated 
data are given in an appendix, without discussion, 


cussion is centered around them. 


but with sufficient explanation concerning their 
arrangement to permit their detailed study by 
(See table 1, page 40.) 

Each circle appearing on the frequency-distribu- 


those interested. 


tion plots represents a value derived from a single 


analysis. The cireles are plotted equidistant on 
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the ordinate corresponding to their value. Thus 
the abscissas are values derived from the analyses, 
and the ordinates indicate the frequency with 
which these values occur. For example, in the 
plot showing the distribution of the results of the 
analyses for carbon dioxide (fig. 2) it will be seen 
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Ficure 2. Frequency-distribution plot for carbon dioride-0.1 
percent interval. 
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that 1 analysis gave 0.6 percent CO, 5 ans \yses 
gave 0.7 percent, 15 analyses gave 0.8 percent. and 
soon. This manner of plotting shows the distr. 
bution of all of the results at a glance, and {rms 
the basis of a probability curve, of which th 
maxima are the most probable values.  Thes, 
graphical pictures not only save hours of study of 
the tabulated data, but actually make the sig. 
nificance of the data so self-evident that lengthy 
and detailed discussion of each plot is hardly 
necessary. 


2. Carbon Dioxide 


All the analyses reported for carbon dioxide ar 
given in this plot (fig. 2). Analyses reported to 
hundredths of a percent were rounded off to tly 
nearest tenth percent. Figure 3 is a separate plo 
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Figure 3.—Frequency-distribution plot for carbon dioride 
0.01 percent interval. 


of those analyses which were reported to hur 
dredths of a percent. (The tendency to round of 
to the nearest 0.05 percent is apparent.) Th: 
upper section of this plot gives the results of fow 
laboratories that have generally beer in agreemen' 
during this series of cooperative analyses, av 
probably represents the best reproducibility to l 
hoped for between laboratories with this type 
sample. The total spread in this group is less thai 
that noted for the whole group and suggests tli 
desirability of taking the little extra trouble | 
adjust pressures and read volumes rather carefull) 
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. greatest frequency occurs at 1.0 percent. 
« 313 analyses reported for this constituent, 
are between 0.9 and 1.1 percent. These 
is are not unsatisfactory, but can be improved. 


3. Oxygen 


The analyses for oxygen are plotted in figure 4. 
These are particularly interesting as there was no 
ienificant amount of oxygen in the sample. One 
 two-tenths of a percent of oxygen is usually 
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Frequency-distribution plot for oxygen-0.1 per- 


cent interval. 
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reported in the analysis of a natural gas, and the 
extent to which oxygen really does occur has long 
been a question. Other constituents of the gas 
will dissolve in the reagents used to remove oxygen, 
and it is accordingly necessary to saturate these 
reagents with respect to the other gases and to 
maintain this equilibrium throughout the analysis. 
If this is correctly done, the error of solubility, 
leading to an erroneously high value for oxygen, 
should be minimized or eliminated. The present 
analyses show the extent to which this procedure 
has been successfully applied. Of the 294 analyses 
reported, 124 (only 42 percent) correctly indicated 
no oxygen; 58 (about 20 percent) indicated 0.1 
percent oxygen, which may be considered the 
average error of solubility that experience has 
shown may be expected; and 69 (about 23 percent) 
indicated 0.2 percent oxygen, which should have 
been the maximum error caused by solubility 
alone. However, 24 and 19 analyses indicated 
0.3 and 0.4 percent, respectively, and these values 
round off the probability curve. Beyond this, 
there are 20 values between 0.5 and 1.5 percent, 
which probably were the result of air contamina- 
tion during the transfer of the sample to the 
analytical apparatus. 

It will be shown later that the heating values 
computed from the analyses indicating no oxygen 
are in general higher than those reported by the 
whole group. 


4. Methane 


The data for methane are given in two plots, 
figures 5 and 6. In the first of these, methane is 
plotted to the nearest half percent; in the second 
plot, 36 values have been eliminated before plot- 
ting to the nearest tenth percent. 

Referring to the first plot (fig. 5), the greatest 
frequency appears at 77 percent, with an average 
of 76.5 percent if values lower than 73.5 and 
higher than 80.0 are discarded. The range 73.5 
to 80.0, which includes about 94 percent of all the 
analyses, is rather wider than could be desired. 
The analyses were generally reported to the 
nearest tenth percent, although some were noted 
to hundredths. When plotted to the nearest 
tenth (fig. 6, section A+B) the whole group loses 
the sharp peak near 77 percent, and the distribu- 
tion becomes widespread. This plot has been 
separated into two sections, and the letters A 
and 8, correspond to the two methods previously 
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Figure 5.—Frequency-distribution plot for methane—0.5 
percent interval. 
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Ficure 6.—Frequency-distribution plot for methane—0O.1 percent intereal. 


noted. The upper section gives the analyses these values are still widely spread, they * 
reported by the four laboratories that used more closely grouped than is true for all of " 
method B, which takes account of oxygen con- analyses. (Note that none is off scale.) Gro 
sumed during the combustion as well as the con- B indicates a somewhat lower average valv 
traction and carbon dioxide produced. Although 76.1 percent, for methane. 
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S. Ethane cent. About 89 percent of the analyses lie in 

the range 17.0 to 21.5 percent, with an average 

The analyses for ethane are given in two plots, value of 19 percent. The analyses are plotted to 
figures 7 and 8, which follow the same scheme the nearest tenth percent in section A+B of 
used for plotting methane. The first plot (fig. 7) figure 8. The distribution is again wide, even 
shows the greatest frequency to be at 18.5 per- more so than is shown, for 31 values have been 
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Frequency-distribution plot for ethane—0.5 
percent interval. 
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FicureE 8.—Frequency-distribution plot for ethane—0.1 percent interval. 
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dropped. The values noted in the upper section 
of this plot, yielded by method B, are more 
closely clustered and are higher than the whole 
group. The average is 19.9. 


6. Nitrogen 


Nitrogen is always of especial interest in a 
combustion analysis, because no matter what 
method is used for the combustion the nitrogen 
(or more properly, the inert fraction), should 
always be the same. It is a measure of complete- 
ness of combustion and may indicate leakage 
during sampling or the analysis. In the case of 
method B, nitrogen is measured as an unreacted 
residue, while it is computed by difference from 
100 percent in method A. The analyses for 
nitrogen are plotted to the nearest tenth percent 
in figure 9. This scale corresponds to the scale 
used in figures 4, 6, and 8. A condensed scale 
(0.5 percent) for nitrogen, to correspond to the 


seale of figures 5 and 7, has not been presented, 


7. Calculated and Measured Specific Gr« vity 


The accuracy with which the specific g: ayity 
was computed from these analyses has becn ¢. 
rectly determined, since the specific gravity of this 
standard sample was measured by E. C. Creit, 
of the National Bureau of Standards,’ Thy 
measured value was definitely known to one & 
two orders of magnitude more than could be ey. 
pected from the analytical data. The measure 
value (referred to dry CO,-free air=1) was 0.682) 
+0.00005, and this has been indicated on tly 
plot (fig. 10), which shows the values calculate 
from all of the analyses, plotted to the neares 
0.002. The calculated values were all compute 
from the following values given in the Interna. 
tional Critical Tables: 

. 5290 
. 1053 
. 97208 
. 5544 
. 0493 
. 5624 








6 vawes 
> OFF SCALE 
@ O06 
os,070e 


@eeeeeeeeee 
@eeeeeseoeoeeee 
eeeceeeoooooooe 
eeeeeeeee 





eeeeceeeoeeeoe 
eeeeeeees 








Lh] 


NI 
c. 


PERCENT 
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even though it was not possible to include all of 


the analyses in figure 9. Note that six values 


are off scale to the left. 

The plot for nitrogen shows a rather wide dis- 
tribution. Values reported for method B are 
more closely grouped than is found in the over-all 
picture, but the average, 3.0, differs only slightly 
from that of the whole group, which is 3.1 for the 
range 1.5 to 5.0 percent. 


Frequency-distribution plot for nitrogen 


NITROGEN 


0.1 percent interval. 


The greatest frequency is noted at 0.672, low 


than the measured value by 0.010 or 1.5 percet! 
The calculated specific gravities are plotted to th 
nearest 0.001 in figure 11, section A+B. 
though the greatest frequency is still seen! 
0.671 to 0.672, the distribution now appears rath 
wide. Section B shows the values given 
method B, and in this case the distribution is s@ 


* See footnote 5. 
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Figure 10.—Frequency-distribution plot for calculated specific gravity-0.002 interval. 
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Fiaure 11.—Frequency-distribution plot for calculated specific gravity-0.001 interval. 


ide, although the range is considerably improved. A comparison with the specific gravity measured 
{ 37 of the 272 values are eliminated, the by the various instruments ordinarily used for 
verage of the whole group from 0.657 to 0.685 this determination shows a generally better score 
pbout 88 percent of all values reported) is 0.672, for laboratory type instruments designed for this 
r 1.5 pereent too low. The average of the B measurement, but a somewhat spotty score for 
oup is 0.675, somewhat closer to the measured recording instruments. Three specific-gravity 
alue, although 1 percent too low. balances yield average values of 0.680, 0.682, and 
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0.681. An effusion apparatus gave the average 
value 0.689, 1 percent too high. One recorder 
gave the correct value, 0.682, but three others 
recorded the values 0.690, 0.687, and 0.696, too 
high, by 1.2, 0.7, and 2.0 percent, respectively. 


8. Calculated and Measured Heating Value 


Because the heating value of this standard 
sample was measured, it is possible to estimate the 
accuracy with which the heating value was calcu- 
lated from the analyses. Unfortunately, this 
estimation cannot be made as closely as desired 
because of the uncertainty in the measured heating 
value. Two series of measurements were made 
with the Junkers calorimeter, the first by J. H. 
EKiseman and the second by R. S. Jessup. The 
following values were obtained: 


Btu/ft®, ° C, 30 in. He 


1100. 
1097. 
1098. 


(1102. 4) 
(1099. 1) 
(1100. 2) 
Average 1098. § (1100. 6) 
1106. 

1108. ! 
1107. ¢ 


(1107. 8) 
(1109. 5) 
(1109. 1) 
Average 1107. § (1109. 1) 
Average of 
series 


(1105) 


* Since this report was first issued for study by the ASTM Committees 
involved, Mr. Jessup has recomputed the measured heating values according 
to the new procedures of ASTM Tentative Method of Test for Calorific 
Value of Gaseous Fuels by the Water-Flow Calorimeter. This tentative 
method has been approved by Committees ASTM D-Sand D-3-III, ASTM 
designatian D-900-46T. The values in parentheses are those obtained by 
the new computation, which makes a correction of the measured heat of 
combustion to the standard temperature of 60° F. The original values noted 
on the frequency distribution plots have not been altered, although official 
action of the ASTM will probably make this desirable eventually. 


In addition to the Junkers determinations, other 
independent measurements were made. The 
Thomas recording calorimeter gave the value 1,100 
Btu/ft® at the time the first Junkers series yielded 
the value 1,099. In addition, the values 1,099 and 
1,102 were obtained from a measurement of the 
heat of combustion of this sample and a measure- 
ment of the carbon dioxide produced on combus- 
tion. These values were derived as follows: 


30 


(1) The heat of combustion was measu!od by 
J. W. Knowlton according to the proc dun 
described for the determination of the heats oj 
combustion of methane and carbon monoxide 
Three calorimetric experiments were performed q 
a mean temperature of 25° C, with the water vapor 
but not the carbon dioxide, removed from tly 
sample prior to combustion. The heat of com. 
bustion so determined was 842.68+ 0.20 inter. 
national kilojoules per mole (44.010 g) of (C0, 
collected after combustion. 

(2) The stoichiometric relationship betwee 
sample and CQO, present after combustion wa 
determined volumetrically by burning the samp} 
(containing its original CO,) over hot platinum \: 


the presence of excess oxygen in an all-glasd 


apparatus.” Seven determinations were made }y 
Shuford Schuhmann. 

Volume of CO, present after combustion 9 
1.0000 volume of sample ASTM-2: 


. 1679 
. 1677 
. 1672 
. 1678 
. 1665 
. 1681 
. 1669 


Average. 1. 167440. 0005 


The tabulated values represent the measure 
decrease in volume when the carbon dioxide we 
absorbed from the mixture of oxygen, nitrogen, aw 
carbon dioxide after the combustion of 1.00) 
volume of gas, corrected for the measured carbo 
dioxide originally in the sample. The measur 
ments were made at approximately 25°C and 7 
mm Hg, corresponding to the conditions durin 
the measurement of the heat of combustion 
The volume fraction of the carbon dioxide in th 
mixture was about 0.6. 

(3) By combining measurements (1) and (2 
the heating value of the sample may be calculate 

Converting int. kj/mole of CO, to Btusmole 


842.68 . ; 
‘ ° = 7QaC 5 / » ¢ hf] 
CO,: 1054607 99-95 Btu/mole of CO, pre 


after combustion. 


* Frederick D. Rossini, The heat of formation of water, BS J. Researctt 
(1931) RP259; The heats of combustion of methane and carbon monoxit 
37 (1931) RP260; The heat of formation of water and the heats of combust 
of methane and cirbon monoxide, A correction, 7, 329 (1931) RP345 

* A modification of the analytical apparatus described in BS J. Resear 
121 (1931) RP266, with the manometer-compensator described in J. Kese" 
N BS 26, 351 (1941) R P1382, and all rubber connections eliminated 
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volun 
f CO, 
ther h 
eviati 


ressur 


d by 
dures 
its of 
xide.! 
ied at 
apor 
nm the 
com. 
inter. 
CO 


verting Btu/mole of CO, to Btu/ft® of CO, 
ideal gas state at 0° C and 760 mm of mer- 


rv «799.05 Btu/22.4140 liters@1009.46 Btu/ 


s.3162 liters or per ft* of CO). 

Converting this expression to the condition 30 in 
273.16 762—13.25 
288.71 760 
940.95 Btu/ft® of CO, (ideal gas) present after 
urning this sample, at 30 in Hg, 60° F, saturated. 


, 60° F, saturated: 1009.46 


It now remains to correct the carbon dioxide 
easured volumetrically after combustion for 
eviation from ideality. This correction intro- 
uces the greatest uncertainty in the present 
leulation: 
The partial pressure of CO, in the residual 
olume after combustion was approximately 0.6. 
nder the conditions of the experiment (nearly 
C and 1 atm), the deviation of CO, from ideal- 
y is 0.0051, which is to say that it occupies only 
9949 of the volume 1.0000, which represents the 
eal condition. However, the actual deviation 
CO, mixed with the oxygen and nitrogen in the 
roducts of combustion is not equally well known. 
me rough measurements made under essentially 
¢ same experimental conditions ® indicate a 
artial deviation of 0.0004; and if this is accepted, 
volume of sample is equivalent to 1.1679 volumes 
{ CO, in the products of combustion. On the 
ther hand, the conventional assumption, that the 
eviation of CO, is proportional to its partial 
ressure, would fix the partial deviation as 0.0030, 
id the corresponding relation would be altered, 
volume of sample now being equivalent to 1.1709 
olumes of CO, after combustion. The two 
eating values calculated from these two stoichio- 
etric relations are 


(940.95) (1.1679)=1,099 Btu/ft® (1,100) " 
940.95) (1.1709) =1,102 Btu/ft® (1,103) 


hus the uncertainty in the deviation of CO, from 
leality corresponds to 3 Btu, whereas the uncer- 
hinty in the stoichiometric determination itself 
as only 0.7 Btu/ft®, and the uncertainty in the 
‘termination of the heat of combustion was 


Joseph R. Branham and Martin Shepherd, Gasometric method and 

pparatus for the analysis of mixtures of ethylene oxide and carbon dioxide, 

earch NBS 22, 171 (1939) RP1175. (See p. 185.) Joseph R. Branham, 

| Shepherd, and Shuford Schuhmann, Critical study of the deter- 

of carbon monoxide by combustion over platinum in the presence 
nyeen, J. Research N BS 26, 571 (1941) RP1396, (See p. 586.) 

footnote “a,” p. 30. Parenthetical values given hereafter refer to 

footnote 
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reported by Rossini as the equivalent of 0.2 
Btu/ft*. 

An attempt was made to take advantage of the 
accuracy of the determination of the heat of 
combustion by a gravimetric measurement of 
the CO, present after the combustion of a known 
volume of the sample. A number of such deter- 
minations were made by E. C. Creitz, of the 
National Bureau of Standards. The heating 
values calculated from these determinations and 
the measured heat of combustion were 

Rtu/fe 


1101.7 
1100.3 
1101.7 
1101.6 
1100.7 
1099.4 
1104.2 
1104.1 
1102.3 
1103.0 


1101.9 + 1.2 (1102.8) 


Several determinations obviously out of line, 
and all much lower, have been omitted from this 
average. There was no direct evidence that all of 
the gas had always been burned. However, some 
support is given the value 1,102 (1,103) Btu/ft*, 
derived from the stoichiometric experiment. 
These values are within 1 Btu of the average 
determinations with the Junkers calorimeter.” 
For the present, the average heating value, 
1,103 +4 Btu (1,105), obtained with the Junkers 
instrument will be accepted for comparison with 
the analyses. This value is noted in the two 
frequency-distribution plots that give the data 
for heating values calculated from analyses. All 
the data are plotted to the nearest 2 Btu/ft® 
in figure 12. The greatest frequency is noted at 
1,104 Btu. The data are plotted to the nearest 
1 Btu in figure 13, with 25 values (8 percent of 
those reported) omitted. Section B of this plot 
gives the values calculated from analyses per- 
formed by method B. These values are more 
closely grouped than the whole and are higher. 
‘2A standard sample of this type for the cooperative determination of 
heating value, using the procedure given in the new tentative standard 
proposed by Subcommittee D-3-III of ASTM on Determination of Calorific 
Value of Gaseous Fuels, might well be distributed. It would then be possible 
to compare the reproducibility of the measured heating value with that cal- 
culated from analysis. Another profitable excursion would be the accurate 
determination of heating value, either by improving the gravimetric measure- 


ment of carbon dioxide present after combustion, or by means of a more 
accurate calorimeter. 
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Figure 12.—Frequency-distribution plot for calculated heating value—2.0 Btu interval. 
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Ficure 13.—Frequency-distribution plot for calculated heating value—1.0 Btu interva’. 


The average of the B group is 1,108 Btu/ft®. which is itself expressed within +0.4_ peree 
The average of all the computed values taken from The thing desired of the calculated values 
1,086 to 1,113 Btu (about 85 percent of all values —s greater frequency at and near the true heat! 
reported) is 1,101 Btu/ft®. These averages for value. The thing desired of the measured val’ 
the calculated values are within 0.5 and 0.2 per- _ is a closer approach to whatever the true va 
cent, respectively, of the average measured value, — may be. 
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9. Laboratory Averages 


<» far the picture has included all of the deter- 
A 
icture of each laboratory with relation to the 
Figure 14 is a 
lot showing the average values obtained by each 


jinations made by all of the laboratories. 
thers may now be presented, 


»boratory for each of the constituents determined. 
.veral of the laboratories reported more than one 
verage value, and these have been included as 
parate points on the plots. Points with a slight 
ii] to the right indicate average values derived 
rom method B. All others were derived from 


too wide, and a need for standardization is indi- 
cated. 

The average calculated specific gravities and 
heating values reported by each laboratory are 
plotted in figure 15. All values reported are 
given in the lower sections, marked A+B+C. 
As before, the sections marked B give values 
obtained by the laboratories using method B. 
The section marked C gives the results of the 
laboratories which reported no oxygen in the 
sample. Apparently there is no correlation be- 
tween “‘no oxygen” and calculated specific gravity. 
However, those who found no oxygen reported 
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ethod A, with the very few exceptions to be 
ted in the next section of this report, which will 
ve results obtained by the three special methods 
eviously noted, 

The plot for methane has not included four 
lues higher than 78.8 percent, and that for 
hane has omitted four values lower than 17.0 
ut. In general, the spread of the averages is 
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ME THANE 


Frequency-distribution plot for laboratory averages of all constituents. 


heating values which were more consistently 
grouped around the measured value. This sug- 
gests the removal of some hydrocarbons during 
the absorption analysis, and again the inclusion 
of some air in the samples transferred for analysis 

an effect seen in the tendency for these values to 
wander to the left of center, in comparison with 


the C group. 
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10. Analysis by Method II 


As previously stated, this method is ordinarily 
used for the analysis of a carburetted-water gas 
and involves an absorption of unsaturated hydro- 
carbons in fuming sulfuric acid, as well as a 
fractional combustion over copper oxide preceding 
the combustion over the platinum catalyst in the 
presence of excess oxygen. 

Laboratory 21 used this method and reported 
the following average results: 

Percent 
CO, 1.0 
Op 0 
N; 5.3 
C,H, . 6 
H, 
Cs 
CH, 
C,H, 
Specific gravity 
Btu/ft*®. 


0. ! 
74. ! 
17. 6 

0. 6754 
1. 065. 9 

The analysis is interesting for many reasons. 
The CO, and O, are correctly reported. The high 
nitrogen is not explained by air contamination, 
and probably resulted from incomplete combus- 
tion or a faulty exchange of atmospheric inerts in 
the reaction tubes. 

The amount of C,H, found was high, and indi- 
cated loss of other constituents in the fuming 
sulfurie acid. Actually, there was about 0.2 
percent of propylene in this sample. Laboratory 
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Frequency-distribution plot for laboratory averages of calculated specific gravity and heating value. 


1 determined unsaturated compounds by absor 
tion in 10 drops of fuming sulfuric acid on gl 
wool, using as part of the volumetric apparatus 
reaction tube described by Shepherd." 
age of six determinations was 0.19 +0.04. 1) 
was checked later by the mass spectrometer (s 
following section “Analysis by Mass Spectromete 
of this paper), which indicated the presence 


.0.18 percent of C,H, as the only unsaturate presev’ 


The H, and CO reported were actually » 
present in the sample. The absence of hydrog 
was definitely established by a separation at (| 
temperature of liquid hydrogen. The absence 
carbon monoxide was established by a sensiti 
chemical test." 

Both the methane and ethane found were | 
low. The calculated specific gravity was |o 
The calculated heating value, however, was | 
minimum noted for the whole group, 37 Btu lov 
than the average measured value. 

These results taken alone afford no recomm 
dation for the double combustion method. On‘ 
contrary, this refinement of the single combust) 
method has not yielded the desired results. 

Laboratory 24 did not follow method II } 
did make an absorption in fuming sulfuric « 
“18 Modification of apparatus for volumetric gas analysis, J. Re 
NBS 26, 351 (1941) RP1382. See page 3M, figure 2, tube number 2 

'* Martin Shepherd, A preliminary report on the N BS colorimetric iné 


ing gel for the rapid determination of small amounts of carben mon 
Anal. Chem. 19, 77 (1947). 
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ng the contraction as C,;H,. The averages 
for two series of 10 analyses each were 0.6 
(5 percent. The complete results follow: 
Series 1 Series 2 
Percent Percent 
1.0 1.0 
0.0 0.0 
7 
. 5 
. 6 
3.2 


0. 672 
1103 


pecific gravity 0. 675 
u/ft >... 1106 
‘orreet values for CO, and O, are reported. The 
',H, is apparently lowered at the expense of the 
‘C,Hy.” The specific gravities are a bit low, but 
he heating values are very good. 


ll. Analysis by Distillation Methods X and Y 


Laboratory 28 analyzed the sample by method 
i, but in addition separated the sample by frac- 
ional distillation and noted in a secondary report 
he result of the correction for propane found by 


he distillation method. The data are 


Average Values 
by Method A 
corrected for 


Average Values 
by method A 


a 
= 


CO, 
N; 
) 
CH, 
C Hg 
Cy He, 


—=wooS © 
ore Sr oe Gr Ge Gr 


Specific gravity 
Btu/ft * 


rhis is very satisfactory work. The slight cor- 
rections applied did not materially alter the cal- 
ulated specific gravity or heating value—and did 
ot improve the calculated specific gravity. 
Laboratory 29 employed method Y. The com- 
ustion data were thus not calculated as such but 
‘orrected for CsH, found by fractional distillation. 
The average values obtained by this laboratory 


Percent 


0.8 
. 24 


Analysis of Standard Sample of Natural Gas 


Percent 
CH, 76. 98 
CLHg¢- 14. 1 
C;Hg.. - 3.0 
No 4. 83 


Specific gravity —_ - 0. 6838 


Btu/ft® : .. 1091.6 


Carbon dioxide is a bit too low, oxygen is high, 
and nitrogen is higher than the average value 
by a large enough amount to suggest air contam- 
ination. The calculated heating value is low. 
The calculated specific gravity is within 0.002 of 
the measured value. The method should have 
improved the accuracy of the calculated heating 
value, and did not. 

Laboratory 30 employed method Y, 
depended upon distillation for all constituents 
except CO, and O,. These average results were 
reported: 


which 


Percent 

CO, : 1, 00 
O» 0. 00 
CH, 80. 6 
CyH, 14.8 
CsHg *3.6 
No 0.0 


Specifie gravity 0. 6736 
Btu/ft® 1155. 2 
* Actually, the distillation indicated 2.4 percent heavier than CgH,, but 
this was called C3H, in the computation of specific gravity and heating value. 


The distillation method has long been favored 
over the combustion method for the analysis of 
natural gases, and it is well understood that the 
latter method cannot yield the true composition. 
It is accordingly surprising to find the maximum 
calculated heating value of the present series of 
analyses has been derived from the distillation 
procedure. This heating value is about 50 Btu/ft* 
too high. The specific gravity is 0.008 too low. 
The absorption analysis has reported the correct 
values for carbon dioxide and oxygen. 


12. Analysis by Mass Spectrometer 


The first two analyses by the Consolidated mass 
spectrometer were made under the direction of 
A. Keith Brewer. Notable amounts of hydrogen 
were reported, and the results were not concordant. 
The calculated heating value and specific gravity 
were not in agreement with the measured values 
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for the first analysis, but in good agreement in 
case of the second analysis. 

Because of the apparent presence of hydrogen, 
which is not to be expected in a natural gas in 
anything like the amount found, a separation 
of this sample was made at the temperature 
of liquid hydrogen. This experiment definitely 
proved the absence of hydrogen (i. e., the amount 
present could not have exceeded 0.001 percent). 
Accordingly, the original fragmentation patterns 
of the hydrocarbons used in the analyses by mass 
spectrometer were redetermined, and the cor- 


rected patterns were then used in computing the 
mass spectrogram obtained in a third analysis. 
This analysis and computation were made by 


Vernon Dibeler. The three analyses are 


Third 
analysis 


First Second 
analysis analysis 


Percent Percent 


0. 74 
76.7 
3. 90 


Percent 
H,. 0. 24 
CH, 68. 02 
N2- 13. 93 
CH, 12.55 | 14.7 
CH, 0. 22 0. 17 
CO, 1. 22 4 
C;H,_-.. 2. 72 2. 90 
Other hydrocarbons 0 0 
0. 680 0. 682 
1,102 1, 102 


Specific gravity - — - . 719 
Btu/ft® 992 


Although the first analysis is out of line, the 
second and third are in good agreement (with H, 
dropped from the second), and both yield excellent 
values for specific gravity and heating value. The 
C,H, found checks the chemical determination of 
unsaturated hydrocarbons made by laboratory 1. 
The correct value for CO, was reported, and no 
was found. The nitrogen is somewhat 
higher than the value indicated by chemical 
analysis. The 2.82 percent of C,H, found com- 
pares reasonably with the values 3.0, 3.1, and 3.6, 
reported by three laboratories using distillation. 


oxygen 


For ready comparison, the values obtained by the 
mass spectrometer (third analysis) are plotted in 
figures 14 and 15 as open circles. Methane and 
ethane are not directly comparable to the amounts 
found by combustion; and the propane and propy- 
lene do not appear. 
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13. Errors in Calculated Specific Gravity ang 
Heating Value 


In view of the additional data furnished by ily 
special methods just discussed, it will be well » 
reexamine the specific gravity and heating valy« 
calculated from the analyses. The presence of 3) 
percent of C,H, and 0.2 percent CH, in this sam. 
ple seems well established. (The figures hay 
been rounded off to the needed significance.) |: 
these C; hydrocarbons were computed as CH, an 
C,H,, it is well to check the errors thus introduce 
in the calculated values. 

Propane is calculated thus in the analysis: 


C;H,—2C,H,—CH,. 


Accordingly, the calculated specific gravity wis 

too high by 0.002: 
1.5624-—+2(1.0493) — 0.5544 

(0.03) (1.6442 calculated — 1.5624 actual) =0.00% 


Similarly, the calculated heating value was too low 
by 0.2 Btu/ft*: 


2528—»2(1759) —997 


0.03 (2528 —2521) =0.21. 

The calculation for the propylene depends i 
turn on the manner in which methane and ethaw 
were calculated. 

(1) When CH, and C,H, are calculated from 
total contraction and carbon dioxide formed, thes 
stoichiometric relations are used: 


CH, 
C,H, 


1/3(4TC—5CO,) 
1/3(4CO,—2TC). 


With propylene injected, TC and CO, are e 
pressed thus: 


TC 2CH, T 2.5C.H,+4 2.5C,H, 
CO,: -CH, T 2¢ ‘Hy + 3C3H,. 


Accordingly, CH, and C,Hg, calculated as abow 
are in error: 
CH, -CH,—5/3C,H, 
C,H, C,H, + 7/3C 3Hg. 


The 0.2-pereent propylene thus makes the ca: 
culated CH, too low by (5/3) (0.2)=0.33, and th 
C,H, too high by (7/3) (0.2)=0.47 percent. 
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Similarly, when CH, and C,H, are caleu- 
lated from the total contraction plus carbon 


dioxide and the oxygen consumed, according to 


the relations 


CH, z} 3(7(7C€ T CO,) 90, | 
C,H, =1/3[60,—4(7C+CO,)], 


and the (T04+-CO,) and O, are expressed 


T0+-CO,)>3CH,-+4.5C,H,+5.5C3H, 
O, »2CH,-4 3.5C,.H,+ 4.5C3H,, 


the caleulated CH, and C,H, are in error thus: 


CH, -CH,—2/3C,H,=0.13 percent (in this case) 
C,H, C,H, +5/3C3;H,=90.33 percent (in this case), 


The errors can now be evaluated: 

(1) When CH, and C,H, are calculated from 
TC and CO, alone, the presence of 0.2 percent 
of C;H, causes the calculated specific gravity to 
be too high by 0.0002 (which is not significant): 
(0.0047) (1.049) —- (0.0033) (0.5544) —0.002 (1.45) 
0.0002. 

The calculated heating value is 0.3 Btu/ft® too 
high: (0.0047) (1759)— (0.0033) (997.4)— (0.002) 
(2339). 

(2) When CH, and C,H, are calculated from 
(TC+CO,) and O,, 0.2 percent of C;H, makes the 
calculated specific gravity too low by 0.001 (not 
significant): (0.002) (1.45)—[(0.0033) (1.049)— 
(0.0013) (0.5544)] and the calculated heating value 
too low by 0.2 Btu/ft*®: (0.002) (2339)—[(0.0033) 
(1759) — (0.0013) (997.4)}. 

The errors for propane and propylene may now 
be combined: 

1) The calculation from 7C and CO, makes the 
caleulated specific gravity too high by 0.002, and 
the calculated heating value too high by 0.1 
Btu/ft*. 

2) The ecaleulation from (7C+CO,) and O, 
makes the caleulated specific gravity too high 
by 0.002 and the calculated heating value too low 
by 0.4 Btu/ft®. 

The errors in calculated heating value are of no 
interest, once they are known. The error in 
Specific gravity does not improve the situation, 
since the correction would move the calculated 


value further from the observed one. 
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14. Comparison of Analysts in Same Laboratory 


The extent to which the different analysts using 
the same apparatus and method will achieve con- 
cordant results has often been discussed and with 
varying opinions. One school holds that analysts 
are born and not made; while a second and more 
optimistic group believe that any high-school 
child can be quickly converted to the desired robot. 
The analyses submitted in this cooperative series 
prove nothing conclusively, but offer interesting 
comparisons to show that analysts using the same 
apparatus and procedure can be together or apart, 
which was to be expected. 

Two frequency distribution plots have been 
prepared to give this picture. It was not worth- 
while to plot each component of the mixture. 
Nitrogen and the calculated heating value were 
selected to illustrate these agreements and dis- 
agreements; nitrogen because it is a measure of 
the over-all behavior of the chemical analysis, 
and is comparable no matter what procedure has 
been used; and heating value because it is another 
measure of the over-all success of the analysis, 
and of especial engineering interest. 

The plots showing the comparative results of 
analyst versus analyst are shown in figures 16 
and 17. Each laboratory that contributed to 
this competition is numbered and separated by 
heavy horizontal boundaries. The analysts who 
entered each laboratory competition are lettered, 
and their individual efforts are separated by 
lighter horizontal boundaries. 

By reading these plots we find, for example, 
that analysts A and D in laboratory 2 were on 
good terms. Figure 16 for nitrogen shows no 
evidence for great concern until laboratories 22 
and 23 are examined. It is obvious that in both 
laboratories an analyst, M, had some difficulty in 
deciding upon a suitable value. Figure 17 shows 
in general a somewhat greater disagreement. 
There are differences between analysts A and B 
of laboratory 5, A and B of laboratory 14, and M 
and N of laboratory 22, which might well be 
straightened out. Laboratories 2, 17, 18, 20, and 
24 appear to have arrived at agreement with 
themselves if not always with each other. Here 
then is a brief picture showing the expected 
differences in talent. 















































a0 
PERCENT NITROGEN 
Fiaure 16.—Frequency-distribution plot comparing different 
analysts in the same laboratory with percent nitrogen as 
the basis for comparison. 
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Fiaure 17.—Frequency-distribution plot comparing differen 
analysts in the same laboratory, with calculated heating 
value as the basis for comparison. 
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Comparison of Laboratories Within Same 
Organization 


|: might logically be expected that laboratories 
‘thin th esame organization would have com- 
ared results and composed their differences, but 
he story seems to be much the same as for the 
ndividual analysts. The data available offer no 
ther information than this, for their number is 
pot sufficient to tempt a statistical statement as 
o the proportion of agreement and disagreement 
ithin companies. 

These data, again for nitrogen and heating 
alue, have been plotted in figures 18 and 19. 
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roure 18. 
aboratories in the same organization, with percent nitrogen 


Frequency-distribution plot comparing different 
as the basis for comparison. 


here are three groups or organizations repre- 
nted, divided on the plots by heavy horizontal 
undaries; the laboratories within each group 
re separated by lighter horizontal boundaries. 

Perhaps there is less reason to expect agreement 
nong the four laboratories in group 1 than in the 
ther two cases, for this is a group of Government 
boratories. Laboratories 5, 6, and 7 are in the 
ame Bureau, and 5 and 6 are in the same build- 
ig; but each laboratory employs different appara- 
s and methods. Laboratory 1 is in another 
The agreement between laboratories 1 
(5 in the matter of percentage of nitrogen is 
extended to laboratories 6 and 7. The agree- 
ent between laboratories 1 and 7 in the calculated 
ating value is not extended to laboratories 5 


ureau. 


nalysis of Standard Sample of Natural Gas 
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Figure 19.-—Frequency-distribution plot comparing different 


laboratories in the same organization, with computed 


heating value as the basis for comparison. 


In the case of laboratories 18 and 19 of group 2, 
both in the same company, and both operating in 
connection with different phases of the same gen- 
eral project and in the same location, there is only 
fair agreement on the percentage of nitrogen, 
and no good agreement on heating value. Both 
laboratories employed like apparatus and methods. 
In the case of laboratories 20 and 21 of group 3, 
both laboratories were in the same company, but 
were widely separated geographically, and used 
different methods. There is no agreement between 
these laboratories. Indeed, it was not found pos- 
sible to plot the results of laboratory 21 on the 
seale selected to include the other laboratories. 

These cases pictured in the last four plots serve 
to emphasize the need for standardization, and 
the further need for equally expert analysts. 


VI. Conclusions 


The data presented here give the first clear 
picture so far available of the actual status of gas 
analysis of this particular type in this country. 
As such they are worthy of considerable reflection 
on the part of gas chemists and plant engineers. 
In view of the considerable effort expended by 
everyone who took part in the analyses, and 
because of the importance of the data, the com- 
mittee believes that the results should be released 
for study by all laboratories engaged in this type 


39 








of volumetric gas analysis, even though the data 
need to be supplemented by further investigation. 
Some very creditable work has been presented 
and some that should have been better. The 
data indicate a considerable degree of standardiza- 
tion in the use of modern equipment, but, never- 
theless, offer evidence to support the claim that 
further important standardization is certainly 
needed to clear up the analytical picture, especially 
with respect to operating technic and methods. 
In particular, the relative accuracy and repro- 
ducibility of methods A and B herein described 
should be determined more thoroughly. Further 


Tasie lt. 


Each laboratory has been assigned a number, and the reports are arranged 
in the order of these numbers. At the left hand, just under the laboratory 
number, the method used is designated by the proper key letter At the 
right hand, under the laboratory number, is the key designating the ap- 
paratus used. These keys are given below. 

Referring to the tabulated data, column | gives the analysis number, 
column 2 the analyst, columns 3, 4, 5, 6, 7 give in order the percentage by 
volume of CO), Os, CHy, Cols, and N». Column 8 is used to report other 
hydrocarbons, or other gases, which may have been determined by special 
methods. Column 9 gives the total for all the constituents determined 
Columns 10 and 11 give, respectively, the heating value (in Btu/ft’,¢7.4) and 
specific gravity calculated from the analyses. The measured heating value 
was 1,103 +4 Btu/ft® (20 in, Hg, 60° F, saturated), and the measured specific 
gravity was 0.6820 —0.00005, referred to dry C O--free air. 

The keys to methods and apparatus follow, 


Method A 


1, Absorption of carbon dioxide in an aqueous solution of potassium 
hydroxide. 

2. Absorption of oxygen in alkaline pyrogallol or in a chromous solution. 

3. A single combustion over hot platinum in the presence of excess oxygen, 
measurement of the contraction and carbon dioxide produced, and 
computation of all hydrocarbons present as C H, and C2He, using the 


conventional formulas 


CHy=1/344 TC—5C Op») 
CyHe= 1344C O,—2TC), 


where 7'C is the total contraction on burning, and CO, is the carbon diovide 
produced by the combustion. 
4. Nitrogen by difference. 


Method B 


1. Absorption of carbon dioxide in a solution of potassium hydroxide 
2. Absorption of oxygen in alkaline pyrogaliol. 

(One of the four laboratories then took fresh portions of the sample for 
combustion, correcting the combustion data for the CO, and O, 
originally found. The other three used the general procedure of com- 
bustion of a portion of the sample taken for absorption.) 

3. A single combustion over hot platinum in the presence of excess oxygen, 
measurement of the contraction, carbon dioxide produced, oxygen 


correlation of the results obtained by the abso 
tion-combustion and the distillation methods is y 
indicated need. An investigation of the accura 
and reproducibility of that important new analy, 
cal tool, the mass spectrometer, is another jini 
cated need. And quite aside from the analytic, 
one or two standard samples might well be dist, 
buted to cooperating laboratories for the det 
mination of heating value, using the recent) 
proposed ASTM Tentative Standard. In th 
connection, the need for a very accurate meth, 
to determine heating value is becoming mo 
apparent. 


Analytical data reported by the cooperating laboratories 


consumed, and the unreacted residue. Computation of all by} 
carbons as CH, and C)H, using the formulas: 


CHy,=1/3[7(TC+C Oy —90)) 
CH= 1/3[602:—4(TC+C 02). 


4. Nitrogen reported as the measured unreacted residue. 


Method X 


1. CO, and O, by absorprion in a separate sample. 
2. A single combustion of the hydrocarbons, with measurement 
traction and carbon dioxide produced. 

3. A separation of the hydrocarbons by distillation in one of the wa 
rectifying columns 

4. Correction of the combustion data to take account of hyrdocarhe 
heavier than C:Hs. 


Method Y 


1. CO») and Oy; separately by absorption. 
2. Other components by distillation. 


Apparatus 


Volumes are made comparable by means of a pressure-temperis 
compensator with manometer interposed between the compensa 
tube and the burette. 

Pressure within the burette is balanced against existing barome 
pressure, and gas volumes are corrected from the observed prew 
and temperature to a common basis, including a correcticr 
changes in the saturation pressure of water. 

Pressure in burette is balanced against atmospheric pressure bu 
correction is made for changes in pressure or temperature dur 
analysis. 

V.—Volumes are measured by observing the pressure exerted wit! 
constant volume. 

R, — Pipettes are connected to the burette by a manifold. 

Ry—A single pipette, connected to the burette, serves for all reaction 

Ry — Pipettes are temporarily connected in succession to the burette © 
different reactions are progressively conducted 

Ry —The burette itself serves as a reaction tube. 

Hg— Mercury is used as the confining fluid 

H,O—(with appropriate subscript) An aqueous solution serves « 

confining fluid 
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TasLe 1.—Analytical data reported by the cooperating laboratories—Continued 


LABORATORY NUMBER 1 
Apparatus V;R,Hg-Shepherd 





Percentage by volume Calculated 


Anal. No. Analyst = 
CH, Cie Ne (CoH. Total Btu/ft * sp er 


0. 22 
17 
19 
21 
20 
17 


100. 01 
100. 00 
100.00 
100. 00 
100.01 
100. 00 
99. 98 
100. 00 
100. 00 
100. 00 


100.00 


Apparatus V;R,Hg-Shepherd 


Percentage by volume Calculated 


CH, CyHy No % Total Btu/ft * sp er 


76.05 19.60 100. 01 1103. 28 0. 6761 
vy. US 1107. 38 6762 
100. 01 1105. 74 OTM 
100. O1 1107, 82 6750 
76. 35 19. 64 100. 00 1106, 98 6740 
76. 67 19. : 100. 01 1108. 49 6728 
77. 24 18. { 2. 82 100. 01 1103. 90 HOS 


76. 02 19.85 
76. 26 19. 62 


76, 24 19.75 


en 


77.01 19. 2! 100. 01 1107.00 6722 
76. 45 19.72 100. 01 1100.00 6736 
75. 80 19 3. 100. 05 1106. 00 6762 


76. 41 19. SS loool 1106. 16 0. 6741 


LABORATORY NUMBER 3 
Apparatus V,R)-Shephe 


100. 01 
100.00 
100. OL 
100. O1 
Oo. Ol 
100. 02 


100. 08 
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TABLE 1.— Analytical data reported by the cooperating laboratories—Continued 


-ABORATORY NUMBER 4 
Method R Apparatus ViRiHg shepher 


Percentage by volume Caleulat 


Anal. No. Analyst 


CH, Coils N: (CsHs) Total Btu/ft * 


19. < 100. 0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100 0 


Fee we wow ww 


106. 0 


-ABORATORY NUMBER: 
Method .1 Apparatus V,;R:Hg-Bureau of M 


100.0 om 
100.0 
100.0 
100. 0 
100.0 


100.0 

100.0 

100.0 

100.0 

100.0 

100. 0 " Meth 


LABORATORY NUMBER 
Method .1 Apparatus V«R:Hg-Bone and Whe 


100. 0 
100.0 
100.0 
100.0 
100. 
100 
100 
100 
100 
100 


“oO mm ww « 


= 


> 
Se Be BNW Ww ww 


on 


2 Tinea H 
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TaBLe 1.—Analytical data reported by the cooperatiing laboratories—Continued 


LABORATORY NUMBER 7 
Apparatus V,R,;He-Allen 


Percentage by volume Calculated 
Anal. No Analyst . es Cae 
CH, Coll No (CyHs) Total Btu/ft 4 sp er 


1115.1 0. 6762 
1107.9 6765 
1104.4 6720 
1110.3 6776 
1109. 2 6745 
1109.6 

1110.2 6700 
1107.8 6788 
1103. 7 6742 
1108.7 . 6788 
1000. 2 6742 


ww ww 


Sw ww hw 


1107.8 0. 6761 


LABORATORY NUMBER 8 
Method A Apparatus V,R,Hg-Burrell 


100. 00 . 105 0. 665 
100, 00 . tis 
100 
100 
100 


LABORATORY NUMBERY# 
Method A Apparatus V);>R;Hg-Hempel 


1, 006 0. 673 
1, 105 678 
1, 07 676 
1, 109 679 
1, 14 679 
1, 104 674 
1, 108 678 
1, 109 thie) 
1, 107 671 
1,104 . 163 


1 
9 
2 
1 
2 
” 
” 
1 
2 
1 


1, 105 0.674 


-ABORATORY NUMBER 10 
Apparatus V,R,iH¢—Burrell 


. 678 
679 
671 

. 61 
65 
a9 

. 671 
660 


. 82 


SET STON Se me a TR renege 


STOP RANE 6 


0. 470 
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TaBLe 1.— Analytical data reported by the cooperating laboratories —Continued 


LABORATORY NUMBER Il 
Method A Apparatus \ | Ri Hy 


Percentage by volume Caleulat«{ 


CHy Cols No (C3Hs) Total | Bru/ft® 


Anal. No. Analyst " - 


16. 70 80. 00 2.00 100. 00 1, 082 
17.20 78. 80 2.70 ee 100. 00 1, O88 
17.00 79. 30 2.30 _ | 100.00 | 1,090 
17. 50 78. 30 2.80 100.00 | 1,089 
16.90 | 79. 9 1.85 100.90 1,092 
17.00 | = 79.95 1.65 ; 100. 00 1, 096 
16. 80 80. 10 1.70 : 100. 00 1, 095 
17.10 79.80 | 1.80 : 100. 00 1, 096 
17.40 73.70 | 2.50 | : 100. 00 1,092 
35 16.9 | 79.9 | 1.80 ; 100.00 | 1,004 


0.35 17.0% 79. 48 2.11 10. 00 1, 092. 4 


eeuocwesewnw- 


3 


- 
<4 
n 


LABORATORY NUMBER 12 
Method A Apparatus V,R;— Fischer 


100.0 | 1,001 
100.0 | 1,108 
1. 100.0 
0. 
4 
06 | 
(by diff.) 


LABORATORY NUMBER 13 
Method A Apparatus V,R,Hg—Burrell—U. 8. Stee! 


! 


LABORATORY NUMBER Ii4 
Method A Apparatus V,R,Hg¢—Elliott—Burr 


100.0 1,004 
100.0 1, O89 
100.0 1, O86 
100 1, 093 
100. 1, O86 
100 1, 007 
100 1, 004 
100 1,004 
100 1, 007 
100 1, 102 


“a - 


wean ww Ww 
“cw © Ww 


ower e#ewnwe 
eo 


~ 


100. 0 1, 093 
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TABLE 1.— Analytical data reported by the cooperating laboratories—Continued 


LABORATORY NUMBER 15 
Apparatus V;R; R:H,0-NaSO«vH S80. 
Modified Williams 


Percentage by volume Calculated 


Anal, No. Analyst 
CHa Cole | Na (Css) Total Btu/ft * sp ger 


20 5 Not 
Detd. 100 0. 6735 
20. 2 do 100 ; 6718 
22. do 100 6809 
pi) f do 100 Hi 
wo do 100 i A728 
21.3 § do 100.0 6769 
22. 2. do 100.0 6910 
do 100.0 M. 6831 
do 100.0 30. § 6742 
..do 100.0 { 6670 


NWN N NH WH Ww 
a a a 
se dsiz 


| 
x 


6767 


fethod A Apparatus V:.RiHg—U. 8. Steel, Modified 


LABORATORY MUMBER 17 


Method A Apparatus V,R\He-U 


ll ll a a 
oe & hw 


“Oar wersiss tw ® 


www Ww Ww hb 


yotnotes at end of table. 
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Tasve 1.—Analytical data reported by the cooperating laboratorices—Continued 


LABORATORY NUMBER 18 
Method A Apparatus V,H,Hg-Bureau © Mine 


Percentage by volume Caleula! 
Anal. No Analyst . , . 
CH, Cos N (CyHs) Total Btu/ft * 


18.9 
19. 2 


te 


1105. ! 

1101 

1106, 

1103. ¢ 

1110.! 

1107.5 

1105. 

1103.7 

1105. 2 this 
1107.9 (66. 


w 


cfc ee we +1 +t oS & 


“NNN NW Wt 


= 


1106. 1 6% 


100 1105.8 
100 1105.0 
100 1108.8 
100 1105. 2 
100. 1109.8 
100.0 1107.8 
100.0 1107.3 
100.0 1105.8 
100. 1108.3 
100 1108.2 


wenweenweewn www 


100.0 1107.2 


LABORATORY NUMBER 19 
Method A Apparatus VR, He 


0 hi 
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TABLE 1.—Analytical data reported by the cooperating laboratories—Continued 


LABORATORY NUMBER 2 
Apparatus V;R,Hg-U. 8, Steel 


Percentage by volume Calculated 


Anal. No Analyst 
CH, Clg Ne (Cys Total Btu/ft sper 


78 
78 
SI 
79. 
78 
79. 
78 


Not detd 0. 6613 
do 6626 
do 6512 
do Hnoy 
do , oH . 623 
do 6504 
do 6639 


17.6 
18. 2 
15.7 
17.5 
17.8 
17.3 
18.0 
18.9 


2 \ wo 


Se 
no em me 


do 6705 


i ee 


tS ww wD 
: a 


78 do HHSS 


78 7 ) 0. 6620 


0. 6615 
6549 
6714 
6504 
HOSS 


worn 
wi24 
6609 


ww ww hb 


fne60 


0. 6628 


LABORATORY NUMBER 21 
Apparatus V,R\Ke-Burrell 


centage by volume Calculated 


Anal. No, 
CoHs Ny ‘3 ; : Btu ft 


1070. 6 0. 6684 
1061. 5 6805 
1058, 2 6765 
1057.7 6748 


1071.5 6766 


1065. 9 


alysis of Standard Sample of Natural Gas 
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TABLE 1.—Analytical data reported by the cooperating laboratories—Continued 


LABORATORY NUMBER 22 
Method A Apparatus V;R,Hg and V:R,H.0.Nos 


Percentage by volume | Calcul 


Anal. No. Analyst — ce ae Vs - 
Cis Na (CyHs) Total Btu/ft ! 


1, 108 
1, 099 
1, 109 
1, 121 
1,101 
1, 100 
1, 102 
1,112 
1,113 
1,113 


ee 


1, 108 


ABORATORY NUMBER 23 
Method A Apparatus V,;RiHe 


3 
3 
9 
2. 
3. 
3 
4 


exec 


> ‘ 
2 
2 3.1 
2 1 
2 2.1 
2 2.3 
2 3.4 
2 44 
2 44 
2 2.5 
2 ‘ 
> ‘ 
? 


eon Oe we 
— ww Ww 
<= 


“oO 
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TABLE 1.—Analytical data reported by the cooperating laboratories—-Continued 


ABORATORY NUMBER & 
Apparatus V,R,Hg-Burrell 


Percentage by volume Calculated 


Anal, No. Analyst 
CHy Cite Na Total Btu/ft sp er 


0. 6714 
6723 

. 6751 
6727 

. 6749 
6749 

. 6762 


6779 


6785 


674 


6751 


1105 
1101 
1101 
1104 
1102 
1102 
1102. 2 
1105 
1104. 5 
1106. ! 


we te 


1103. ! 





BER 25 
Apparatus V,R,Hg-Burrell 


HRH : 6774 
HRH 7 : 5 6671 
HRH 76. 1: i 5. 6745 
HRH | 3.3 6646 
HRH 7.4 3. 82 ‘ fin | 
HRH 76 3.6 ; 6702 
HRH 7. 6 : 6663 
HRH mM 5 { 64g 
HRH 7. 4: : 7 6675 
HRH 4 3. 6 : 6706 


Heyl 


LABORATORY NUMBER 2% 
Apparatus V,R,Hg-Burrell 


645 
fi45 
659 
i590 
643 
OM 
(41 
(41 
656 


0.90 
sl 
vo 

1.00 

1.00 

O83 
vo 
w 

1.12 

0.90 


one KN HN & & & Ww 


657 


0.93 7% 3.17 3 . 0. 649 
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TaBLe 1.—Analylical data reported by the cooperating laboratories—Continued 


LABORATORY NUMBER 2 
Method A Apparatus V,R 


Percent by volume Caleul 
Anal. N Analyst 
‘Hs rotal¢ Btu/ft * 


1002. SH 
1109. 27 


1107.84 


LABORATORY NUMBER 2& 
Apparatus VR, Hg-Burrell-Oberfell and Distillation 


Erroneo 
Percent by volume Cak jas Saln 
otal co 
lotal Bin ft Propane 
Total hy 
Obtaines 
100. 00 1110.2 ObLtaine 
100. 00 1102.5 ‘ Heating 
100. 00 1107.6 ws \nalyst 
100. 00 1106. 4 
100. 00 owe 1 ‘ yASHID 
100.00 1111.0 
100. 00 1102.7 
100. 00 nee 0 
100.00 1092.9 


100. 00 lov | 


100000 1103.0 


LABORATORY NUMBER 29 


Method X Apparatus V;R;—Goek 
Distillation Apparatus 


100.0 love2 
100.0 1OsS9 
100.0 1091 
100.0 1090. 


100.0 1092 


Ave 76.4 : 100.0 1091. 6 
Maximum : 77.1 5 1092.8 


Minimum 7 ri r 1089.8 


See footnote at end of table 
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TABLE 1.— Analytical data reported by the cooperating laborateries—Continued 


LABORATORY NUMBER 30 
Apparatus V,R, He and 
distillation apparatus 


Percentage by volume Calculated 


Anal. No Analyst . 
CH, “2H No , Potal Bru/ft sp gr 


100. 00 D 0. 6817 


100. 00 (651 


100. 00 1155.2 0. 6736 


Erroneous result; apparatus not completely flushed with nitrogen in preparing for analysis 

jas sample exhausted after this analysis 
ital combustibles; numbers |, 2, and 3—wet manifold; numbers 5, 6, 7, and 8—dry manifold 

Propane assumed absent 

lotal hydrocarbon volume determined by difference, nitrogen directly determined. 

Obtained from average of absorption analyses 

Ottained from fractionation analysis of small sample. 

Heating value computed, using given factors and analysis on dry basis 

(Analyst B did not have enough time to complete any more of the analyses. 


Average heating value and specific gravity computed from average composition, 
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urification, Purity, and Freezing Points of 8 Nonanes, 
11 Alkylcyclopentanes, 6 Alkylcyclohexanes, and 
4 Butylbenzenes of the API-Standard and API-NBS 


Series” 


By Anton J. Streiff,? Evelyn T. Murphy,’? Janice C. Cahill,? Helen F. Flanagan,’ 
Vincent A. Sedlak,? Charles B. Willingham, and Frederick D. Rossini 


This report describes the purification and determination of freezing points and purity 
of 29 hydrocarbons of the API-Standard and API-NBS series, including 8 nonanes, 11 


alkyleyclopentanes, 6 alkvleyclohexanes, and 4 butylbenzenes. 


I. Introduction hydrocarbons [4].** The material labeled API- 
NBS is made available in appropriate small lots 
A previous report described the purification, on loan to qualified investigators for the measure- 
nd determination of freezing points and purity ment of needed properties.* 
{ 37 hydrocarbon compounds of the API- 
standard and API-NBS series, which were pro- II. Materials 
uced as part of the cooperative program on 
tandard samples of hydrocarbons of the National 
Bureau of Standards and the American Petroleum 
nstitute. This report describes the purification 
nd determination of freezing points and purity 
an additional 29 hydrocarbon compounds under 
his cooperative program, including 8 nonanes, 6 3,3-Diethylpentane. 
lkyleyelohexanes, 11 alkyleyelopentanes, and 4 2,2,8,5-Tetramethyipentane. 
utvibennenaé ® Ethyleyclohexane (B). 
fe. 7 4 1, 1-Dimethyleyelohexane. 
The final lots of the material labeled API- cis-1,2-Dimethyleyclohexane. 
standard are sealed “in vacuum”’ in glass ampoules trans-1,2-Dimethyleyelohexane. 


nd made available as NBS standard samples of Ethyleyclopentane (one-half) (/). 
prem trans-1,3-Dimethyleyclopentane. 
sented before » Divis De . namnintew . - 
, ore ¢ ne ] ivision of Petroleum Chemistry of the American n-Propyleyelopentane. 
\| Society, Chicago, Ill., September 1946, wl 
lsopropyleyclopentane. 


n brackets indicate the literature references at the end of this 
1,1,3-Trimethyleyclopentane (three-fourths). 

estigation was performed at the National Bureau of Standards —_—__—————— 
he work of the American Petroleum Institute Research Project 6 ‘ The allocation of the API-NBS samples of hydrocarbons is handled by 
ysis, Purification, and Properties of Hydrocarbons the Advisory Committee for the API Research Project 44 on the “Collection, 
h Associate on the American Petroleum Institute Research Project Analysis, and Calculation of Data on the Properties of Hydrocarbons” 
tional Bureau of Standards. (W. E. Kuhn, chairman), 

perative program was carried on under the API Research Pro- 5 (B) following the name of a compound indicates that, for the API-NBS 

mittee on Hydrocarbons for Spectrometer Calibration (W. J series, it is a second (and usually slightly purer) sample of the given com- 

irman). Further details are given im references [1, 2, 3). pound, the first sample of which is labeled (A). (See reference [5)). 


The starting materials were supplied as follows :° 

By the API Research Project 45 on the “Synthe- 
sis and Properties of Hydrocarbons of Low Molecu- 
lar Weight” at the Ohio State University, Columbus, 
Ohio, under the supervision of C. E. Boord: 


Puri cation and Purity of Hydrocarbons 53 





By the Hydrocarbon Laboratory at the Penn- 
sylvania State College, State College, Pa., under 
the supervision of F. C. Whitmore: 

2,4,4-Trimethylhexane. 

Ethyleyclopentane (one-half) (B). 
1,1-Dimethylevclopentane. 
cis-1,2-Dimethyleyelopentane. 
trans-1,2-Dimethyleyclopentane. 
1,1,2-Trimethyleyclopentane. 

cis, cis, trans-1,2,4-Trimethyleyclopentane. 
cis, trans, cis-1,2,4-Trimethyleyclopentane. 

By the National Advisory Committee for 
Aeronautics, through its Aeronautical Engine 
Research Laboratory at Cleveland, Ohio, and the 
Automotive Section at the National Bureau of 
Standards: 

a-Butylbenzene 

Isobutylbenzene 
sec- Buty lbenzene 
tert-But ylbenzene 


(B) 
(B). 


: | 2,2,3,4-Tetramethylpen- 
| 
(B). | 


tane. 
2,2,4,4-Tetramethylpen- 
tane. 
2,3,3,4-Tetramethylpen- 
tane. 


(B). 


TABLE 1. 


Hydrocarbon 
e for 
distillation 


Starting material * 


Compound * provided by 
Veh | a ee 
ume Purity | Kind 


By the Standard Oil Development Co. Ein, 
beth, N. J., through W. J. Sweeney: 


cis-1,4-Dimethyleyclohexane. trans-1,4-Dimet} yley; 


hexane. 


By the Gulf Oil Co. Fellowship at the \elly 


Institute of Industrial Research, Pittsburgh, p, 
through W. A. Gruse: 


1,1,3-Trimethyleyclopentane (one-fourth) 


By the API Research Project 6 at the Nation, 
Bureau of Standards, under the supervision | 
F. D. Rossini: 

n-Nonane. 2,2,5-Trimethylhexane. 

Table | summarizes the amounts of the starti 
materials and gives some additional informati« 
as to the source and purity. 


Information on the purification of 29 API-Standard and API-NBS hydrocarbons 





Volume oi 


Distillation / selected samp: 


Amount 
of hydro- 
car 
int 


| Num- 
- - Reflux 
ratio / 
(ap- 
prox.) 


Dis- 
| tilling 
he 


Results 
plotted 
in 
figure 


ate 
of col- 
lection 
of dis- 
tillate 


Areotrope- 
forming 
»@ 
substance tropic 
distil- 
late « 





PARAFFINS 


n-Nonane APIRP6« 


« 
volume 


49 





2,2,5-Trimethylhexane 
2.4.4Trimethylhexane 


2,2,3,3-Tetrameth ylpentane 
2,2,3,4-Tetramethyl pentane 


2,2,4,4-Tetramethylpentane 
2,3,3,4-Tetramethylpentane 


3,3- Diethylpentane 


n- Butyl benzene (8) 
Isobut ylbenzene (B) 
sec-Butylbenzene (2 


fert-Butylbenzene (B 


Foot notes on page 56. 
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APIRP6* 


Penn State 


APIKP45 


NACA 


do 


do 


APIRP45 


75 

. 09 
73 
0.73 


Si 


61 


Me. Cell 


Coll 


Cell 


Cell 


ALKYLBENZENES 


Reg 
do . 
Azeo Me Carb 
do do 
Reg 


Areo Me Carb 


~_- = 


J i on on 





trans 


vievelo 


PaBLe 1. 


pound « 


sion thyleyclopentane (B 


|-Dimethyleyclopentane 
»1,.2 Dimethyleyelopen- 
tane 

ne-1,.2-Dimethy leyclopen- 
tane 
ns-1,3-Dimethyleyclopen- 


Starti 
rmati 


tane 


Propyleyclopentane 


»propy leyclopentane 


|,2-Trimethyleyelopen- 


tane 


3-Trimethyleyclopen- 


tane 


is, trans-1,2,4-Trimeth- 


icvelopentane 


cis-1,2,4-Trimeth- 


trans 


cyclopentane 


vieVciohe\ane 


methyleyclohexane 


thyleyclohexane 


methyleyelohex- 


hyleyclohexane 
ethyleyclohex- 


Starting material * 
provided by 


APIRP45 
Penn State 


Penn State 


do 
do 


APIRP4S5 


do 

do 

do 

do 
Penn State 


Gulf-Mellon 


APIRP45 


Penn State 


APIRP45 
do 
do 


do 


do 


Oil Dev 


do 


stad 


Hydrocarbon 
charged for 
distillation 


Vol- 


— 
ume | Purity 


Information on the purification of 29 API-Standard and API-NBS hydrocarbons 


Kind « 


Continued 


Volume of 


Distillation selected sample 


Amount 
of hydro- 
carbon 
in the 
azeo- 
tropic 
distil- 
late « 


Num- 


Dis- 
- ber of 
tilling theo- 


col- 

umn 
num- 
ber / 


Rate 

of col- 
lection 
of dis- 
tillate 


Results 
plotted 
in 


Reflux 


Azeotrope- ratio ‘ 


forming 


subs “et 
ubstanc Que 


ALKYLCYCLOPENTANES 


Mole ©, 
99.47 | 
96. 29 | 
09. 86 
99. O58 
9. 90 


99. fis 


73. 6 
90.0 
4.6 


95.1 
97. 06 
OS. 67 


O85 
oS. 9 
98.5 


99. 66 


Reg. 


Azeo 
do 
do 


do 


Reg 
Azeo 
do 


Reg 
do 


do 
Areo 


Reg 


do 


Azeo 


do 

do 
Reg 
Azeo 


do 
Reg 


Azeo 


% by 


volume mil/hr 


wi «614.0 


Ethanol : 3 155, 8.0 
do 165 4.5 
do BS ’ 165, 4.5 


do 


Methanol 
Ethanol 


lsoprop 


Me. Cell 
Isoprop 


do 


Isoprop 
do 


Ethanol 


ALKYLCYCLOHEXANES 


| 


99. 50 


Reg 


rification and Purity of Hydrocarbons 


Ethanol 


‘| lsoprop 
[Cell 
Ethanol 

do 




















* (B) following the name of a compound indicates that for the API-NBS 
series, it is a second (and usually slightly purer) sample of the given com- 
pound, the first sample of which is labeled (A). See reference [6]. 

* The abbreviations represent the following laboratories: APIRP45; 
American Petroleum Institute Research Project 45 (formerly the American 
Petroleum Institute Hydrocarbon Research Project) at the Ohio State 
University, Columbus, Ohio. Penn State; Hydrocarbon Laboratory at the 
Pennsylvania State College, State College, Pa. NACA; National Advisory 
Committee for Aeronautics, Aeronautical Engine Research Laboratory, 
Cleveland, Ohio, and the Automotive Section, National Bureau of Standards 
Washington, D.C. Std. Oil Dev.; Standard Oil Development Co., Eliza- 
beth, New Jersey. Gulf-Mellon; Gulf Oil Co. Fellowship at the Mellon 
Institute of Industrial Research, Pittsburgh, Pennsylvania. APIRP6; 
American Petroleum Institute Research Project 6 at the National Bureau of 
Standards, Washington, D. C. 

© The abbreviations are Azeo., azeotropic; Reg., regular. 

4 The abbreviations are Cell., Cellosolve (ethylene glycol monoethy! 
ether); Me. Cell., methyl! Cellosolve (ethylene glycol monomethyl ether); 
Me. Carb., methy! Carbitol (diethylene glycol monomethy! ether); Isoprop., 
isopropanol. 

« Approximate value obtained from the actual volume of hydrocarbon 
recovered by extracting the azeotrope-forming substance with water in 
separatory funnels. 

f See reference [4] for further details. 

* Obtained by distillation from an Oklahoma petre'eum, 


Ill. Purification 


The procedure followed in the process of puri- 
fication and determination of purity was the same 
as that described in the previous report [1]. 

In addition to the name of the laboratory sup- 
plying the starting materials, table 1 and its 
footnotes give complete information for each 
distillation for each of the compounds. 

Details of the distillation apparatus and opera- 
tions are described in reference [4]. 

Figures 1 to 53, inclusive, show graphically the 
results of the distillations listed in table 1. These 
figures give, as a function of volume of hydro- 
carbon distillate, the refractive index (np at 25° 
C, to +0.0001°), the boiling point of the dis- 
tillate (at the controlled pressure of 724.5 mm 
Hg, to +0.01° C), the freezing point of selected 
fractions of hydrocarbon distillate (in air at 1 atm, 
usually with a precision near +0.003° C), and 
the purity of the hydrocarbon distillate. The 
letters W, YX, Y, and Z indicate the disposition of 
the material, as follows: W, returned to the labor- 
atory supplying the material; XY, blended for re- 
distillation; Y, used for the API-Standard mate- 
rial; Z, used for the API-NBS material. 

As demonstrated in the previous report, the 
blending of fractions of distillate for the prepara- 


TaBLe 1.—Information on the purification of 29 API-Standard and API-NBS hydrocarbons—Continuec 










* Obtained by distillation from commercial “‘pentenes” alkylates (peny. 
+ isobutane). 

‘ This starting material consisted of a blend of 2,2,3,4-tetramet )iylpen, 
(approximately 25%), “S-4" reference fuel (“isooctane”) (approxiny, 
66%), and n-heptane (approximately 9%). 

# Calculated from the measured freezing points of two separate lots » 
were blended together for this charge 

* The total volume of the API-Standard sample was 1,305 m! 

‘ The total volume of the API-NBS sample was 350 ml. 

= This is a second lot of a-propylcyclopentane supplied by Al’ Rese 
Project 45. 

* Fractions 272 to 292 from the distillation of isopropyleyclopentay 
column 10 (see figure 31). 

* This second lot of 1,1,3-trimethylcyclopentane was supplied by th: yy 
Research Project 45. 

» Fractions 41 to 60 from the azeotropic distillation in column 12 (see ty 

« This is a second lot of ethylcyclohexane supplied by the API Ree 


Nonon 
D 5- Tritt 
.4-Trin 










3,3-Te 


Project 45. 34-7 
* This is a third lot of ethylcyclohexane supplied by the API Ree 4.4-Te 
Project 45. 3.4-Te 


* Both cis and trans 1,2-dimethylcyclohexane were obtained from » 
material (see fig. 48). 

* The distillation was begun with isopropanol as the azeotrope-forn 
substance, but because of the relatively small percentage of hydrocartx 
the azeotropic distillate, cellosolve was added to complete the distillatic 

* The number of theoretical plates for this column was not determine 


- Dieth: 






















Butyl bs 
pbutylb 
-Butyll 


tion of material of the highest purity can be dm 
-Butyl 


safely only on the basis of freezing points | 
selected fractions. An example of a case whe 
the purest material is at the very beginning o! 
distillation is shown in figure 27 on trans-|, 
dimethyleyclopentane, and an example of a c 
where the purest material is at the end of th 


hyleyeh 
-Dimet 
1,2-Dir 


ns-1,2-1 


distillation is shown in figure 17 for see-butylbeqfiie.s-1 
‘ Propyle 
zene. propyl 
,.2-Trim 
Trim 
IV. Freezing Points, Cryoscopic Constanipm *""" 
and Purity 
Table 2 gives the following information { 
hyleyek 


each of the 29 compounds, except as otherwise i 
dicated: The kind of time-temperature curv 
whether freezing or melting, used to determine (i 
freezing point [7]; the freezing point of the actu 
sample, in air at 1 atm [7], for both the Af’ 
Standard and API-NBS lots; the calculated val 
of the freezing point for zero impurity {7]; ' 
value of the cryoscopic constant, determi 
from the lowering of the freezing point on | 
addition of a known amount of an appropri 
impurity [7]; and the resulting calculated amou 
of impurity in the API-Standard and the A! 
NBS material. 


-Dimet! 
-Din 
ns-1,2-D 
,1.4-Din 
ne-1,4-D 


19 


lee 
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1S, If 1S 
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TABLE 2.—Freezing points and purity of 29 API-Standard and API-NBS hydrocarbons 
Kind of Freezing point of the actual Calculated amount of im- 
time-tem- selected sample, in air at — in the actual se- 
perature latm cted sample ¢ 
observa- ees Freezing point for zero 
Compound * tions used | impurity in air at 
to deter- l atm 
mine the API-Stand- os arma API-Stand- 7 . 
freezing ard API-NBS ard API-NBS 
point > 


Cryoscopic 
constant, A 


PARAFFINS 


°C °C °c deg | Mole % Mole “ 

: F — 53. 555 — 53. 550 — 53. 53540. 010 0. 0888 | 0.08 +0. 04 0.06 +0. 04 
5-Trimethylhexane * and M — 105. 856 — 105, 834 — 105. 7800. 015 4. 0265 D +0.04 14 +0.04 
4-Trimethylhexane M — 113. 434 —113. 430 — 113. 380-+-0. 020 0535 2 +0. 11 -27 20.11 

+ Tetramethylpentane F — 10. 06 — 10.05 —9.90 +0. 05 . 0040 . 064+0. 020 . 060-40. 020 
.3.4-Tetramethylpentane F — 121. 221 —121. 178 —121.09 +0. 05 0027 - 03540. 014 -024+0.014 

» 4.4-Tetramethylpentane *and M —ih. 500 — 6. 580 —6. 54 +0. 08 0273 - 16 +0. 08 -1L +0. 08 
3,4-‘Tetramethylpentane -and M — 102. 137 — 102. 135 — 102. 123.40. 010 0869 =. 0514-0. 087 - 4440. 087 
B. Diethylpentane M —33. 118 — 33.116 —33. 110-40. 005 . 0223 O1I8+0. O11 - 01340. 011 


ALKYLBENZENES 


Butylbenzene (B) —87. 093 — 87. 970-0. 020 e(0.03885) 0.1240. 0.09 +0 
»butylbenzene (A ) 52 —51. 520 - 48 +0.08 0306 1540. 112+0 
re (lor Butylbenzene (BR ) —7 —75. 493 5. 470-+0. 020 . 0808 1240. .07 +0 
inte t-Butylbenzene (R } — 57.876 57. 850-+0. 015 . 0218 06+0 06 +0 


> whet ALKYLCYCLOPENTANES 


g of 
ng 0! 


ns-\ hyleyclopent ane 44 — 188. 452 43540. 010 0.0303 | 0.06 +0. 08 0.05 +0. 08 
a ca Dimethyleyclopentane , S02 —i0. 802 73 +0.04 404 .8 +0. 02 8 +0. 02 
1,2: Dimethyleyclopentane ; — 53. 927 —53. 910 53.85 +0. 04 . 0042 . B10. 016 . 025+0. 016 

of th ns-1,2- Dimethyleyclopentane ’ 63 —117.61 57 +£0.08 #(.0320) .19 40.10 13 40.10 
18-1,3- Dimethyleyclopentans 33. 767 — 133. 758 - 680-40. 020 . 0455 39 +0. 09 35 +0.09 
‘ropyleyclepentane ) 7.379 —117. 378 7. 34040. 020 0511 -20 +0.10 -19 +0.10 
ypropyicyclopentane ) 433 —111. 431 . 3750. 020 . 0846 20 +0. 07 .19 +0. 07 
2-Trimethyleyclopentane 21. 69 —21. 68 4 +0. 08 . 0080 . 015-40. 009 . 012+0. 009 
rimethyleyclopentane ) ab — 142. 55 44 +0. 08 040 -48 +0.32 - 4440.32 

is, (rans-1,2,4-Trimethyleyclopentane 66 —132.4 55 +0. 06 0385 .42 +0. 3540. 23 
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, trans, cis-1,2,4-Trimethyleyclopentane ) 85 — 130. 84 
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wing the name of a compound indicates that, for the API-NBS using the values of the cryoscopic constants and freezing points for zero im- 


1 second (and usually slightly purer) sample of the given compound, purity given in the preceding columns 
umple of which is labeled (A). See reference [6] 4 This cryoscopic constant was determined by the procedure given on 
rm! F indicates freezing and M indicates melting. See reference [7] for ex- page 371 of reference [7] 


on | mental details and the definition of the cryoscopic constant. ¢ Not determined in this investigation. From the “2” tables of the Ameri- 
The lues in this column were calculated as described in reference [7], can Petroleum Institute Research Project 44 [8] 
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Ficure | —Results of the first distillation of n-nonane. 


Regular distillation at 725 mm Hg in still 14 (8/28/45 to 9/15/45) 
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Ficure 2.—Results of the second and final distillation of n-nonane. 


distillation with ethylene glyccl monoethyl ether at 725 mm Hg in still 13 (10/13/45 to 11/14/45 
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Fiaure 3.—Results of the first and only distillation of 2,2,5-trimethylhexane. 


Regular distillation at 725 mm Hg in still 3 (6/26/44 to 8/6/44). 
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Fiaure 4.—Results of the first distillation of 2,4,4-trimethylherane. 


Areotropic distillation with ethylene glycol monoethyl ether at 725 mm Hg in still 13 (11/16/44 to 12/6/44), 
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Figure 6.—Results of the first distillation of 2,2,3,3-letra- 


methyl pentane. 
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2,2,3,3-tetramethyl pentane. 
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Fieure 8.—Results of the first distillation of 2,2,3,4-tetramethylpentane, 


Regular distillation at 725 mm Hg in still 6 (7/27/45 to 9/11/45). 
The portion marked ‘‘N” was discarded, See footnote i of table 1, 
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Ficure 9.—Results of the second and final distillation of 2,2,3,4-tetramethyl pentane. 
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Azeotropic distillation with ethylene glycol monomethy! ether at 725 mm He in still 10 (9/27/45 to 11/19/45). 
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Fiaure 10.—Results of the first and only distillation of 
2.2,4,4-tetramethyl pentane. 


Azeotropic distillation with ethylene glycol monoethy! ether at 725 mm Hg 
in still 4 (9/11/44 to 9/29/44). 
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Fioure 11.—Results of the first distillation of 2,3,3,4-tetramethyl pentane 


Regular distillation at 725 mm Hg in still 8 (7/4/45 to 7/31/45). 
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Figure 12.—Results of the second and final distillation of 
2,3,3,4-letramethyl pentane. 











Azeotropic distillation with ethylene glycol monoethy] ether at 725 mm He 
in still 8 (9/18/45 to 10/16/45). 


Purification and Purity of Hydrocarbons 





°c 


FREEZING POINT 


In °C 


BP ros 


Np AT 25°C 


“33.10 
-33.20 
-33.30 
-33.40 
-33.50 
-33.60 


14440 
144 3% 
14420 


1.4200 
14190 
1.4180 


PERCENT BY VOLUME 
































° 10 20 30 40 50 60 70 80 90 100 
| . T | T T T T T | ~ 
8 
R T PURITY | 3 
q FREEZING POIN Y TT PE 0.00 § 
y ee iit “7 99.80 3 
L. ait 4 99.40 
. | 43350 3 
5 -—{ 98.80 
| z 
, : 
| 
q BOILING POINT | p 
x REFRACTIVE INDEX 3 y 
p — —{3 } 
F 
iv 
|) 
NUMBER OF FRACTIONS | 
Hh — qo - im 37-90 yn 
ae ny “| | 
——_____ —»<_—______—- X w- 
| 
eee eee ee pio cg ng Siok os Bae 5 
0 500 1000 1500 2000 2500 3000 3500 
VOLUME IN ML 


Figure 13.—Results of the first distillation of 3,3-diethylpentane, 


Regular distillation at 725 mm Hg in still 13 (8/13/45 to 9/4/45). 
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Fiaure 14.—Results of the second and final distillation of 
3,3-diethyl pentane. 


Azeotropic distillation with ethylene glycol monoethyl ether at 725 mm He 
in still 13 (11/17/45 to 12/4/45). 
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Ficure 15.—Results of the first and only distillation of n-butylbenzene. 


Regular distillation at 725 mm Hg in still 9 (7/17/45 to 8/18/45). 
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Ficure 16.—Results of the first and only distillation of isobutylbenzene. 





Regular distillation at 725 mm Hg in still 8 (8/27/45 to 9/16/45). 


urification and Purity of Hydrocarbons 





PERCENT BY VOLUME 
° 10 20 3x0 40 so 60 7o 80 390 100 
T T 







































































T T T T Tt ’ ! t 
-7550} 499.9 — 
~7552F T TY 
?-7sseb FREEZING POIN a ‘See - a pees 
z -7556 ne “ | 4oo75 
—-7558F ater | + 
~7560F | 4 pees 
-7562F - \ 
2-564} / 4995, 
N -75.66F / || 4o008 
2 -7568F / 1} s 
* -75 70-4 | 995 
- = 
» a a || 4ee2 2 
+ 4, PERCENT BY VOLUME 
o 2 40 60 80 _ 100 \ 
| ny a ee -e | < 
° wool BOILING POINT— | : 
. = 17 ? | ; 
£ 168 8 t z | hoo 0 & 
; = i| 5 “75485 FREEZING POINTGPURITY | fs 
3! 5 -75.50 — 4o992 
» -75.52F -- ee SE 
£ 1 4900f iJ g -7554° Se: *s 
e .<eer REFRACTIVE INDEX— h 1 wy 775.56 | 49972 
+ | 4880% 4 w 2 
| & \seror NUMBER OF FRACTIONS aD q = : P 
friar a ia} 
eS SS ee , 4 “ 
——<«£o P in.70 a 
VOLUME IN ML . TO 
, : et lati = i71.60- 7 
Figure 17.—Results of the first distillation of sec-butyl- © 
, ' . : 7 | POINT 4 
benzene. fe mt .68 oon. nS Y 
a 171.40R-— | 7 
Azeotropic distillation with diethylene glycol monomethy! ether at 725 mm a 71.30 ~~ 
Hg in still 8 (10/2/44 to 10/19/44). 3 
a 
PERCENT BY VOLUME \$ 
") 20 40 60 100 
a 
7 T oe oe T T T T if bd 14900 \Z = 
o ~~ = 
. © 14890F REFRACTIVE INDEX = |> 
| @ & 14880 ad _ . s 
3 © 1.4870} 4 138, 
° > + 38.48 
> 499.9 5 - 
= -75.82 | Pe NUMBER OF FRACTIONS 
78.56 a FREEZING POIN TY PURITY | tooe & nal an -138.50 
§ -_---s-~ | i 3-56 64 || 
75.56 _—_-- a _ -138.52 
g -75.98- ~~ = a Y z+! 
4 
5-75.60} 7 dose g wh 138.54 
w -75.62- z | 
: 4995 > ° 500 1000 1500 
BOILING POINT - - . ‘ 101.80 
. 169.00} - | 4 Ficure 19.—Results of the third and final distillation 0: 101.70 
. 168.90/- Iw " butylbenzene. 
= 
|3 Regular distillation at 725 mm Hg in still 12 (8/17/45 to 9/4/# 14200 
° 1.4900 F ie a 
o be 4220 
_ — REFRACTIVE INDEX FS " 
z 
= 1.4880 a= 4 1.4200 
e 14870} 7] 4 
I r 1.4180 
| 1.4160 
NUMBER OF FRACTIONS | 
p) 6-54 
*) 
eis | 
1 1 
° 500 1000 1500 2000 


VOLUME IN ML 
Ficure 18.—Results of the second distillation of sec-bulyl- 
benzene. 


Azeotropic distillation with diethylene glycol monomethy! ether at 725 mm 
Hg in still 13 (6/20/45 to 7/12/45). 
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Ficure 21.—Results of the first distillation of ethyleyclopentane. 


Regular distillation at 725 mm Hg in still 15 (9/20/44 to 10/10/44 
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Figure 22.—Results of the second and final distillation of ethylcyclopentane. 
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Figure 24.—Results of the first and only distillation of cis- 
1,2-dimethylcyclopentane. 











Azeotropic distillation with ethanol at 725mm Hg in still 9 (9/29/45 to 10/24/45). 
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Azeotropic distillation with ethanol at 725 mm Hg in still 4 (10/2/45 to 10/28/45)- 
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Figure 26.—Results of the first distillation of trans-1,3-dimethyleyclopentane. 


Regular distillation at 725 mm Hg in still 4 (11/6/44 to 12/1/44). 
Fractions |, 2, and 107 to 140 (marked “X’") were retained for further processing by the API Research Project 4. 
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Ficure 27.—Results of the second distillation of trans-1,3- 
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Figure 28.—Results of the third and final distillation 
trans-1,3-dimethylcyclopentane. 


Azeotropic distillation with ethanol at 725 mm Hg in still 4 (lo2+ 


Fractions 67 to 96 (marked “ 


x” 


11/20/45). 
) were retained for further processin 


the API Research Project 6 


Journal of Research 


rific: 


‘4 





verve ewe |-FlhlUcrMPMrlCUCrlUC( fC 


in MOLE PERCENT (20.09) 


PunITyY 


° 
: 

z 
: 
° 
z 
N 
w 
w 
= 
ve 


BP 75 iN °C 


Np AT 25°C 


FREEZING POINT IN °C 


BP 725 N°C 


Np AT 25°C 


47.37 
-17.38 


PERCENT BY VOLUME 
60 80 





= 


= 


AI7.39F- 


snaof 


7.41 
+I7.42 


129.20 
129.10 


129.00}- 


128.90 


1.4250 
1.4240 
1.4230 


= 
= 
= 


5 


20 40 
T 


. = 2 T | . T 


FREEZING POINT — c PURITY 
= ~ 


BOILING POINT) 





~ 


- REFRACTIVE INDEX~. 





iq 





ae a" 


=W 


| 


FINAL VOLUME 


x 


NUMBER OF rnacrons | 


| 
able 


ma Yaz bbw 


s $ 
~~ @ 
E PERCENT (£0.10) 


3 
oa 
PURITY IN MOL 








.e) 


Ficure_29. 


Regular distillation at 


20 


Results of the first distillation of n-propylcy- 


500 
VOLUME IN “ 


clopentane. 


40 PERCENT BY vee 70 


($00 


725 mm Hg in still 4 (6/12/44 to 6/24/44). 


80 





U 


_——— 





1 


” 


FREEZING denen camels 


—_—_— —_———. __ 


q 


POINT. 


a 


BOILING 


~ 
\ 





fel] $|sshislie te faelas| el ge) see 


REFRACTIVE INDEX 


™ 


NUMBER OF FRACTIONS 


41-57 pe 67-79 $3] 


FINAL VOLUME 





ie) 


ee VOLUME a” 





Ficure 30.—Results of the second and final distillation of n-propylcyclopentane. 


Regular distillation at 725 mm Hg in still 13 (12/7/44 to 12/22/44 
See table | for the composition of the charge for this distillation 
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Figure 31.—Results of the first and only distillation of isopropyleyclopentane. 


Regular distillation at 725 mm Hg in still 10 (7/3/44 to 8/20/44.) 
The portion marked “X"’ was used as part of a charge of n-propylcyclopentane (see fig. 30 and footnote n of table 1). 
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Figure 33.—Results of the second and final distillation of 
1,1,2-trimethylcyclopentane. 
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Figure 35.—Results of the second distillation of 1,1,8- 
trimethylcyclopentane. 


Regular distillation at 725 mm Hg in still 12 (3/13/45 to 4/5/45) 
See table | for the composition of the charge for this distillation. 
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Azeotropic distillation with isopropanol at 725 mm Hg in still 4 (1/3/46 to 
1/20/46) 
See table 1 for composition of the charge for this distillation, 
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Ficure 39.—Results of the first distillation of cis, cis, trans-1,2,4-trimethyleyclopentane. 
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Regular distillation at 725 mm Hg in still 4 (12/5/44 to 1/5/45). 
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Figure 42.—Results of the first distillation of cis, trans, 
cis-1,2,4-trimethylcyclopentane. 


Regular distillation at 725 mm Hg in still 2 (7/19/45 to 9/2/45). 
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Figure 43.—Results of the second and final distillation of 
trans, cis-1,2,4-trimethylcyclopentane. 


Azeotropic distillation with ethanol at 725 mm Hg in still 9 (10/284 
11/23/45). 
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Figure 44.—Results of the first distillation of ethyleycloherane. 


Regular distillation at 725 mm Hg in still 8 (7/20/44 to 8/15/44). 
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Figure 45. 


1500 
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Results of the second and final distillation of ethylcycloherane. 


Azeotropic distillation with ethylene glycol monoethyl etber at 725 mm Hg in still 15 (12/16/44 to 1/3/45). 


See table | for composition of the charge for this distillation. 
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Figure 47.— Results of the second and final distillation of 1,1- 
ses Raa 


Azeotropic distillation with ethanol at 725 mm Hg in still 9 (3/30/45 to 5/7/45). 
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Ficure 48.—Results of the first distillation of cis- and trans-1,2-dimethylcyclohexane. 


Regular distillation at 725 mm Hg in still 5 (1/15/45 to 2/16/45). 
Fractions 7 to 112 were redistilled to obtain /rans-1, 2-dimethylcyclohexane (see fig. 50). Fractions 131 to 185 were redistilled to obtain cis-1,2-dimeth) 
hexane (see fig. 49). 
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Ficure 49.—Results of the second and final distillation of cis-1,2-dimethyleyclohexane. 


Regular distillation at 725 mm He in still 9 (5/7/45 to 6/4/45). 
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Fieure 50.—Results of the second distillation of trans-1,2-dimethyleyclohexane. 


Regular distillation at 725 mm He in still 7 (3/21/45 to 4/23/45). 
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Results of the third and final distillation of 


trans-1,2-dimethylcyclohexane. 


stillation with isopropanol and ethylene glycol monoethyl 
ether at 725 mm Hg in still 13 (9/5/45 to 10/11/45 

ute preceding the point marked “*M”’ was distilled with isopropanol 

trope-forming substance, and the remainder of the hydrocarbon 

ed with ethylene glycol monoethy! ether (see footnote ¢ of table 1). 
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52.—Results of the first and only distillation of cis- 
,-dimeth scl rane. 


Azeotropic distillation with ethanol at 725 mm Hg in still 9 (1/8/46 to 3/7/46). 
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Results of the first and only distillation of trans- 


1,4-dimethylcycloherane. 


Azeotropic distillation with ethanol at 725 mm He in still 10 (1/5/46 to 2/28/46) 


WasuHIncrTon, June 15, 1946. 
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rared Absorption Spectra of Some Experimental 
Glasses Containing Rare Earth and Other Oxides 


By Ralph Stair and Conrad A. Faick 


This paper gives spectral transmission data on soda lime glass containing rare earth 


(Y, Pr, Nd, 8m, Gd, Er) or other coloring oxides (V, Cr, Mn, Co, Ni, U) in the spectral 


region of 0.7 to 4.5 microns. 


I. Introduction 


In connection with experiments upon glass 
wectral filters and upon the physical properties of 
rtain other glasses as a function of their com- 
sition, infrared transmissions were measured 
pon a number of soda-lime-silica glasses contain- 
g some rare earth and other coloring oxides. All 
asses were made with the same basic composition, 
pmely: SiOg, 71.4 percent; Na,O, 18.4 percent; 
ul CaO, 10.2 percent, as computed from the 
itch materials. Calcium and anhydrous sodium 
rbonate of reagent quality and powdered crystal- 
ie quartz were employed for the base glass. The 
her materials were added in the form of oxides. 
These glasses were melted in a 25-ml platinum 
ucible, in a platinum-resistance furnace. The 
belting procedure followed very closely that used 
) previous experimental work [1]* and all glasses 
ceived the same treatment. Data on the spe- 
le refraction and dispersion of some of these 
asses have been published elsewhere [12]. 

The samples obtained often contained a few 
eds and a small amount of striation which may 
ave caused small errors in absorption measure- 
ents. However, as the complete infrared ab- 
rption curve was obtained through the same 


brackets indicate the literature references at the end of this 
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124217477 


area of the sample, relative values between the 
different wavelengths should be adequate. 

The rare-earth oxides used in the production of 
these glasses were from the James collection that 
had been acquired by the Bureau.' According to 
tests made by the Spectroscopy Section of this 
Bureau, these rare-earth oxides contained varying 
amounts of impurities, especially other rare-earth 
materials, but their purity was judged to be 
adequate for the purposes of this report. 

The praseodymium oxide contained small 
amounts of lanthanium and yttrium in addition 
to traces of calcium, copper, iron, magnesium, and 
silicon. The neodymium oxide was of high purity, 
containing only traces of calcium, copper, iron, 
magnesium, and silicon. The samarium oxide 
was especially good, containing, in addition to 
small amounts of calcium and silicon, only traces 
of magnesium and lead, with praseodymium doubt- 
ful. The gadolinium oxide contained about 10 
percent of europium, about 1 percent each of 
bismuth and magnesium, in addition to a small 
amount of calcium, traces of erbium, silicon, and 
samarium, and possibly traces of potassium and 
lead. The erbium oxide was not subjected to as 
thorough spectrochemical study, but it contained 


' Collection of rare-earth materials prepared by the late Professor James of 
the University of New Hampshire and subsequently purchased from the 
University by the National Bureau of Standards 
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about 1 percent of yttrium, smaller amounts of 
holmium and thulium, and traces of scandium and 
ytterbium. The purity of the yttrium oxide is 
not known. 

An examination of the different absorption 
curves shows no significant bands for yttrium 
(see fig. 2) or for europium (see curve 6 for gado- 
linium, which contains Eu as an impurity) nor 
the presence of overlapping of absorption bands 
that might be expected in mixed materials, where 
the same impurities are present for a number of 
the glasses. Hence, in most cases, at least, these 
impurities were in very small amount or else of a 
type that is known to produce relatively weak 
absorption within the region of 0.7 to 4.5 u. 
Reagent quality oxides were used for the other 
coloring materials. 

All glasses having a silica base show a high 
absorption in the region of 2.8 to 4.5 yw, which is 
caused by some combination of silica with water 
and other base materials of which the glass is 
composed. Because this silica absorption may 
be made to vary between wide limits [9] in the 
region of 2.9 u simply by controlling the amount 
of water within the glass, it is impossible to at- 
tempt to make accurate calculations on compar- 
ative transmissions within this silica absorption 
region. Furthermore, analysis of the oxide ab- 
sorption in the region of 2.8 to 4.5 uw is greatly 
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complicated as the result of the high absorp 
in the base glass. Either because of this hig 
and variable absorption in the base glass or becay. 
the oxides studied have no characteristic hg) 
in this region, no absorption bands between » 
and 4.5 » were discovered that could positiy, 
be attributed to the oxides added. 

The transmission of the base glass is give, ; 
figure 1. This glass shows no absorption bay 
at wavelengths shorter than about 2.8 4. Bey: 
this wavelength the silica absorption incre, 
rapidly with increase in wavelength. 

The glasses vary somewhat in thickness. [i 
several of them two thicknesses were availa) 
This is often useful in infrared studies, because 
proper selection of thickness of the specimen pe: 
mits the structure of the absorption spectrum y 
be exhibited to best advantage. In this work th 
more suitable sample was chosen. 


II. Instruments and Methods 


For a detailed discussion of the experiments 
procedure, reference is made to previous publica 
tions [3, 4, 5, 6, 7, 8]. The spectroradiomet 
consisted of a mirror spectrometer, with mirron 
50 em in focal length. A small fluorite prism » 
a Wadsworth mounting was employed. 

The observations were made by the “point } 
point” method. A portable vacuum thermopi 
of bismuth-silver, connected with an iron-cls 
Thomson galvanometer, was used for measuring 
the radiation. The source of radiation was 
Nernst glower, focused upon the entrance slit 
the spectrometer by means of a large coneay 
mirror of I-m focal length. The galvanomet 
readings were read visually and recorded. 1) 
sample was placed in front of the entrance sli 
the spectrometer by means of a mechanical ca 
riage which placed it in the same position for eat 
wavelength. The effective spectral slit widt 
varied from about 0.17 at 2 to 0.07 at 4.5 u. 


Ill. Spectral-Transmission Data 


A brief description is given of the glass and 
special transmission characteristics. A gene’ 
summary covering the kinds and amounts 
oxides added is given in table 1. The glass number 
correspond to the same samples examined in (! 
visible spectrum and reported on elsewhere (2) 
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1.—Amount of oxide added to the base glass 


Index of re- 
Thick- fraction 


Weight of oxide added ness | wavelength, 
| §803.A 


Percent 

Base glass 
3.0 V205 
1.0 CrO, 
1.0 MnO, 
004 CoO 

.05 NiO 

# NIO 
10.0 Y20Os 
10.0 PreOn 
10.0 NdzOs 
10.0 SmeO; 
10.0 GdyO, 
10.0 EreOy, 
3.0 UgOs 


os 
Zz 


l. Yttrium Oxide Glass 


No appreciable absorption in the near infrared 
caused by yttrium oxide (Y,O;) when added in 
» amount of 10 percent (see fig. 2). Neither was 
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Ficure 2. 


y absorption observed within the visible spec- 
im (see reference [2)}). 

\ comparison of the absorption curve for 
‘rum oxide glass with those for the base glass 


Earth Oxide Glasses 


and the erbium and gadolinium oxide glasses (see 
figs. 1, 6, and 7) indicates the possibility of some 
general absorption in the region of 2.7 to 4.2 yu. 
However, differences in thickness and amounts of 
water within the glass may account for the greater 
absorption in this region. 


2. Praseodymium Oxide Glass 


Two strong absorption maxima occur at about 
1.45 and 1.90 uw as the result of the addition of 
praseodymium (Pr,QO,,) in the amount of 10 per- 
cent. Weak bands occur at about 0.76 and 1.0 u 
The general shape of the transmission 


(see fig. 3). 
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Ficure 3. 


curve indicates the possibility of additional bands 
in the region of 2.4 and 4.0 u. Absorption in the 
visible spectrum gives this glass a light yellowish 
green color. 


3. Neodymium Oxide Glass 


The addition of neodymium (Nd,O;) in the 
amount of 10 percent results in three strong 
absorption bands at about 0.77, 1.63, and 2.43 yu 
(see fig. 4). It is to be expected that this material 
would produce a number of absorption bands in 
the near infrared judging from the number and 
intensity of the bands in the visible spectrum [2], 
which give this glass its bright-purple hue. 
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Ficure 4. 


4. Samarium Oxide Glass 


The addition of samarium oxide (Sm,Q,) in an 
amount of 10 percent gives rise to a very strong 
absorption band at about 1.47 yw, and weaker 
bands at about 1.1, 1.92, 2.62, and 4.1 uw (see fig. 
5). The strong band at 1.47 » almost completely 
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masks the one in the neighborhood of 1.1 uy, 
presence of so many bands within this part of t)y 
infrared spectrum suggests the possibility of oth 
bands at longer wavelengths. Hence it would \y 
interesting to study this material incorporate; 
into a special infrared transmitting glass.  Absorp. 
tion in the visible spectrum gives this glass a |igly 
yellowish hue. 


5. Gadolinium Oxide Glass 


No appreciable absorption occurs in the joey 
infrared (or in the visible spectrum) because of 
the addition of gadolinium oxide (Gd,0,) s 
fig. 6). 
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Fiacure 6. 


6. Erbium Oxide Glass 


The addition of erbium oxide (Er,Q,) in 
amount of 10 percent gives rise to one strug 
absorption band at about 1.52 » and two weak 
bands at about 0.80 and 0.95 y, respectively 
fig.7). A comparison with the transmission cu 
for gadolinium (fig. 6), which probably produ 
but little or no absorption in the infrared betwee 
0.6 and 4.5 u, indicates no other important absor?- 
tion bands shorter than 4.5 ~. Absorption in \ 
visible spectrum gives this glass a light-pink cole" 
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Ficure 7. 
7. Other Oxide Glasses 


The addition of vanadium (V,0;), cobalt (CoO), 
ickel (NiO), or uranium (U,Ox,) oxide in small 
ymounts (see table 1) produce strong absorption 
ands in the spectral region between 1 and 4u 
(see figs. 8, 11, 12, and 13). It is to be noted that 
nuch smaller concentrations of these metals were 
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Figure 8. 
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required to produce absorption of magnitude 
similar to that resulting from the addition of 10 
percent of some of the rare-earth oxides. In the 
case of the addition of chromium (Cr,QO;) and 
manganese (MnQO,) oxides, little change in ab- 
sorption occurs within the near infrared region 
(see figs. 9 and 10). 
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Ficure 9. 
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Figure 10. 
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Figure 11. 
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Fiaure 12. 


According to the work of Kiess, Humphreys, 
and Laun upon its atomic spectrum, uranium is 
a rare-earth type element and is the homologue 
of neodymium [10]. It has the same valence 
electron structure and we would expect, therefore, 
its behavior would be similar to neodymium. A 
comparison of figures 4 and 13 is interesting. 

Each of the metal oxides also produces a large 
amount of absorption in the visible spectrum 
resulting in a distinct color for each glass. The 
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colors resulting in the glasses of this group are « 
follows for the different oxides: Cobalt, pale blw: 
nickel, light to dark smoke to dark purple 
vanadium, light yellow to yellow green (similar \ 
novial B); uranium, light yellow (near novia! ( 
chromium, green to yellowish green; and mang 
nese, amethyst. 
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Properties of Water-Repellent Fabrics 


By John W. Rowen and Domenick Gagliardi 


A review and an analysis of the theory of water repellency of textile fabrics have been 


made. 
treated fabrics is discussed. 


repellency of textile fabrics is presented. 


ties of 11 commercial raincoat and 4 military fabrics. 


The physicochemical basis underlying the wettability, or water repellency, of 
A survey of the laboratory test methods for evaluating water 


A study was made of the water-repellent proper- 


For this study two of the more 


recent test methods were examined, the drop-penetration and the contact-angle tests. 


Two other, and older, test methods were also studied, the spray-rating and the hydrostatic- 


pressure tests. 


Several exploratory observations were made in an attempt to determine the 


mechanism by which water-repellent fabrics lose their repellency when exposed to rain. 


I. Introduction 


The National Bureau of Standards has been 
alled upon from time to time to supply informa- 
ion dealing with “water-repellent”’ fabrics. These 
equests have usually come from Government 
gencies and the public at large. In_ recent 
months there has been an intensified interest in 
his subject, and the Bureau decided to review 
he field of “water repellency”’ and to carry out 
laboratory study in order to answer present and 
uture questions regarding this subject. 

An examination of the literature [1] ' revealed 
he existence of a large number of chemical com- 
bounds and commercial treatments about which 
xtraordinary claims were made. Further inves- 
igation showed that a large number of laboratory 
st methods had been devised for evaluating 
‘water repellency.’’ Moreover, considerable con- 
usion was encountered concerning the interpreta- 
ion of the test results and the meaning of such 
erms as ‘“waterproof”’, ‘“‘water-repellent”’, ‘“shower- 
proof’, and “‘water-resistant”’ fabries. It appeared 


of ‘‘water repellency” and to carry out some 
laboratory studies on currently available materials. 
This report, therefore, will deal mainly with the 
fundamental nature of “repellency’’, the labora- 
tory methods used to evaluate it in fabrics, and 
some laboratory observations by the authors on a 
variety of raincoat fabrics. 


II. Theory of ‘Water Repellency”’ 


We must at the very beginning distinguish 
between “waterproof” and “water-repellent” tex- 
tile surfaces. There is still a great tendency to 
mention both terms simultaneously and inter- 
changeably, in spite of the numerous papers [2, 3, 4] 
that have been written to call attention to the 
great difference between them. A_ waterproof 
fabric is one in which the pores, the open spaces 
between the warp and filling yarns and between 
the fibers, are filled with appropriate substances, 
resulting in a fabric having a continuous surface 
and a very low air permeability. A water-repel- 
lent fabric is one in which the fibers are usually 
coated with a “hydrophobic” type of compound, 
and the pores are not filled in the course of the 
treatment. The latter types of fabrics are quite 
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permeable to air and water vapor. The charac- 
teristics of the two types of fabric surfaces are 
summarized below: 


Waterproof 





— 

Water repellent 
| 

| 


| 


| 
| Pores Filled | Unfilled. 


Water-vapor perme- | Very small | Small or large. 
ability. | 
Air permeability Small Usually large. 


| Chief characteristic | Extremely 
to passage of water | 

even under a hydro- | 
static head. 


resistant | Resistant to wetting by 
rain drops and to the 
spreading of water over 
the textile surface, but 
permits the passage of 
water under a hydro- 
| j static head. 








From the foregoing it is seen that a “water- 
repellent” textile surface may never be truly 
waterproof and that a “waterproof” textile surface 
does not necessarily have to be water repellent. 
Unless the pores of a raincoat or shower-resistant 
jacket are large enough to permit transpiration, 
the wearer will be extremely uncomfortable due 
to excessive perspiration. The water-repellent 
fabrics are, therefore, more suited in this respect 
for making raincoats or shower-resistant jackets. 
The subject of waterproofing of textiles is very 
extensive and beyond the scope of this paper, 
therefore the discussion is restricted to the sub- 
ject of water-repellent textile surfaces. 

As the term implies, a water-repellent textile 
surface is one that appears to “repel” water. In 
such cases, the water does not wet or penetrate 
the surface of the textile fabric. The tendency of 
a solid to resist wetting is a function of the chemi- 
cal nature of the solid surface, the roughness of the 
surface, the porosity of the surface, and the pres- 
ence of other molecules on the surface. When a 
drop of water comes to rest on a solid, it may 
assume one of the forms shown in figure 1, or any 
intermediate form from complete wetting to com- 
plete nonwetting. Which form the drop will take 
is determined by the above-mentioned variables. 

Examination of the angle of contact (fig. 1)— 
the angle formed by the tangent to the drop at the 
point of contact with the surface (angle measured 
through the liquid)—reveals that the size of the 
angle is related to the repellency of the surface. 
When the angle of contact is small (under 90 
degrees) the surface is said to be wettable [5, 6], 
and when the angle is large (over 90 degrees) the 
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A, Contact angle between 90 and 180 degrees; B, contact angle of ® der 
C, contact angle between 0 and 90 degrees. 
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reveals that the equilibrium condition is due 
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e following surface energies: Ysa, Yu; Ys, Where 
., is the free surface energy per square centimeter 
numerically and dimensionally equivalent to the 
urface tension) of the solid in contact with air, 7. 
the free surface energy per square centimeter of 
e liquid when in contact with air, and y,, is the 
surface energy per square centimeter of the 
id when in contact with the liquid. At equilib- 
‘um, the following relationship holds: 


Ysa= Ya COS 6- Yst- (1) 


Jupre [7] first derived the equation relating the 
ree surface energies with the work of adhesion: 


W= Ysa Tt Yia— Yat- 2) 
rom eq 1 and 2 it follows that 


W=v7i% (1+cos 6).? (3) 


herefore, by measuring @ and the surface tension 
f the liquid, one may calculate the work of 
dhesion between the liquid and the surface of 
he solid. It is clear that the smaller the work of 
dhesion, the smaller will be the wettability and 
he greater will be the repellency. For poorly 


The work of adhesion ( W’) defined above is the amount of energy necessary 

lestroy | em? of liquid-solid interface and to form 1 cm? of liquid and a 
cm! of solid surface in contact with air. The energy of a unit area of solid 
rface in air may be somewhat different from the energy of a unit area of solid 
urface in vacuum. The difference between these two energies will depend 
pon the relative humidity, the hydrophobicity of the surface, the presence 
transient gaseous or solids in the air, ete. It is physically impossible to 
parate a solid and a liquid without leaving the solid covered with at least 
part ofa monolayer [6, p. 62; 8]; however, it is desirable to define a work for 
ch a process, which is shown to be given by 


W" =YsotVie— Yer (1’) 


here y,. is the free surface energy per square centimeter of solid in vacuum, 
1, is the free surface energy per square centimeter of liquid in equilibrium 
th its vapor, and y,; has the same meaning as before. Now, for the clean 
lid in vacuum, the following relationship (analogous to eq 1) may be 
stulated 

Yer™ Vie COS 0+ Yai, (2’) 
here y,- is the free surface energy per square centimeter of the solid in 
juilibrium with saturating vapor of the liquid. It follows from eq 1’ and 2’ that 


W" = Ye0— Veet viel l+cos 8). (3’) 


If we assume yie= yi (Which approximation is almost within experimental 
ror), we have 

W" = Yeo Yeut W. “@) 
he QUANTITY yee yee may be looked upon as the free energy of immersion, 
constant temperature, of a unit surface of clean solid immersed in a satu- 
ted vapor and is given by [8] 


P 
rere RT| Fr In P, (5’) 


is Gibbs’ “surface density” of the vapor on the solid surface and 
©¢ MZ—T/V, where ¢ is the quantity of water vapor adsorbed per 


ram Of absorbent, M is the molecular weight of the vapor, = is the area in 


Umeters per gram of solid, 7 is the thickness of the surface region, 
e molal volume of the vapor. 
The right-hand member of 5’ can be integrated if an equation for the ad- 
rption isotherm is known or, if not, a graphical integration of the area under 
s! curve can be carried out. 


@ isot he 


repellent, easily wettable surfaces, the drop of 
liquid will assume a shape similar to C in figure 1; 
and for highly repellent, nonwettable surfaces, the 
drop will assume shape A. 

As textile fabric surfaces are not smooth, con- 
tinuous surfaces, but rather porous, screen-like 
surfaces, one must examine the above considera- 
tions and see how they apply to textiles. Figure 
3 shows the cross section of an idealized fabric. 


O09 


A 
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Ficure 3.—Drops on surface of idealized yarn. 


In this figure the parallel yarns are shown as 
circles. If the angle of contact between water 
and the surface of the fabric is larger than 90 
degrees, the equilibrium position of the water 
level of a drop of water will be, as indicated by the 
line A.A, well outside the fabric. If, however, the 
angle of contact is much smaller than 90 degrees, 
the water will penetrate the pores and the level 
will fall to some position such as BB. Cassie and 
Baxter [9] have shown that in the case of porous 
surfaces the apparent angle of contact is related 
to the continuous surface-water angle of contact 
in the following way: 


cos 6,=f; cos @—f2, (4) 
where 


6,—apparent angle of contact 

é=angle of contact 

f,=the fraction of the plane geometrical area 
of unity parallel to the rough surface 
occupied by the solid-liquid interface. 

f,=the fraction of the plane geometrical area 
of unity parallel to the rough surface 
occupied by the liquid air interface. 


They also derived the relationships between the 
“f variables, the distance between the fibers in 
the yarns, the radius of the fibers, and the angle of 
contact: 


r 6 
f=~tq(1- igs) (9) 


Ps 
i= qs 6, 





where 


r=radius of the fibers 

d=one-half the distance between the fibers 

é=contact angle between the water and the 
fiber. 


From the above, it appears as though the most 
water-repellent fabric will be the one on which the 
drop assumes a form as shown in figure 3, A. 
This condition may be obtained by adjusting the 
variables previously mentioned as being responsi- 
ble for the form of the drop. It is advisable to 
consider each of the variables separately. 

Chemical nature of the solid surface.—When puri- 
fied, the natural fibers, cotton, wool, silk, etc., are 
hydrophilic in character and hence the drops 
assume shapes similar to C (fig. 3). It is common 
practice, therefore, to treat fabrics intended to be 
water repellent, with various hydrophobic com- 
pounds—compositions of waxes, petroleum-like 
molecules, and soaps of polyvalent metals, which 
deposit long-chain hydrocarbon molecules on the 
surface of the fabric. When properly treated with 
a water-repellent agent, the surface of the fabric 
will cause a water drop to assume a form very 
similar to A in figure 3. 

Roughness of the surface.—Wenzel [10] pointed 
out that roughness has a peculiar effect on the 
angle of contact. He employed Freundlich’s [11] 
concept of “adhesion tension,’ which is defined as 
follows: 

AT=7ie COs 9. (7) 


It is noted that the adhesion tension is the differ- 
ence between the work of adhesion per square cen- 
timeter (eq 3) and the surface tension. Wenzel 
recognized that eq 7 was true only in the case where 
the surface was a mathematical plane. As most 
surfaces are not of this type, there is associated 
with each surface a roughness factor R, which is 
the ratio of the actual surface to the geometric 


surface. Wenzel showed the validity of the fol- 
lowing equation: 
. cos 6 
AT =V ta R ’ (8) 


Porosity of the surface.—Cassie and Baxter [9, 12] 
pursued the idea of the roughness factor and 
showed that the apparent adhesion tension of a 
porous surface was given by 


A T, = cos Ba ta = Yulfi cos 6—f.), (9) 
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where the symbols have the same mea) ing qs 
eq 4. Using a hexagonal idealized yarn patty, 
as a model, the above workers showed that ¢. 4 
bulk density of the yarn, was related to the ray 
of the fibers in the following way: 
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where p is the density of the fiber. With the gi 
of eq 4, 5, 6, and 10, they were able to objgi 
plots of contact angle versus ¢. These plo 
showed that the angle of contact approaches \w 
degrees as ¢ approaches small values (about ( 
g/ml). Their theory led them to the position thy 
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fabric structure was extremely important in ¢ aa 

production of water-repellent fabrics. bas 
Presence of other molecules on the surface—\\ 

a well-known fact that the surface properties « 

solids may be greatly altered by covering the sw. 

face with impurities, greasy films, dust particls 

water vapor and gases, etc. The presence of sw. 

face impurities may or may not be related to t! 

fact that numerous workers [12, 13, 14, 15, 16, ge" 

have reported not one angle of contact with 

particular system, but two angles, an advancing a 

and receding angle. One angle is obtained wh ®t 

water has advanced over a surface and the othe 

when water has receded from the surface. Th 

difference in contact angle is often called (i 

hysteresis of the contact angle. The cause of |! bari 

hysteresis is still not clearly understood. Kj ?~ 

cently, Harkins [18, 19], working with hydrophohi ware. 

materials and an improved method of measur _ 

contact angles, found no hysteresis. Langu spray 

[17] suggests that the hysteresis of contact angi *" 

especially with water, is due to the presence o! Offci 

surface layer of molecules with one end hyd Spray 


philic and the other hydrophobic, which w& 
overturned by the receding water. 

From the above considerations of the varaiegy 
involved in water repellency, one would expe sr 
that a fabric having a high contact angle wy 
be a water-repellent fabric. In such cases “4 Bune 
water level will remain well ourside the suriie fm * 
as shown in figure 3A. At this point it looks # 
though we may be able to write a prescription / 
the optimum conditions for water repellency, ™ 
a textile fabric in which each fiber is complet Con 
coated with the material having the highest so! on 
water angle of contact. It is unfortunate ‘hf 
this prescription has never been tested. 1 Int 
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for this is that the present methods of 
ating fabries for water repellency do not 
arantee the complete coating of each individual 
or. These treatments leave many of the inter- 
ined fibers without any surface coating. The 
esence of hydrophilic surfaces in the treated 
brie will of course greatly diminish the over-all 
istance of the fabric to wetting. 


ASOn 


_ Survey of Methods Used for Measuring 
Water Repellency in Fabrics 


The number of test methods devised for meas- 


ing water repellency appear to be almost as 


ARLE | 


Reference 
and year 


Test tempera- 


Name of test ture 


CLASS 


Hydrostatic pressure 

Suter hydrostatic pressure 
Mullen (high-pressure range 
Box test for waterproofness 
British (DSIR 


many as the number of laboratories engaged in 
this field. Comparative data from and references 
to the various reports on laboratory test methods 
have been assembled in table 1. The object in 
most of these tests has been to offer a sensitive 
and accurate method of predicting the relative 
performance of fabrics in the rain. Very few, if 
any, of the published papers, however, have fur- 
nished any data showing the correlation between 
the particular test method and the performance 
of fabrics in the rain. Examination of the liter- 
ature reveals that there is a considerable variation 
in the methods of testing and the interpretation 
of the results. It was found convenient for the 


Methods and testing conditions used by various laboratories in measuring water-repellent properties of textile fabrics 


Sample size Result units Extent of use 


i TEST 


em Wide, ASTM, and AATCC. 
Federal Specification CCC-T 
Supplement 
Do 
Limited 
Obsolete 


em 19la and 


rue Ib/in. 
432x20" 
4.8 cm? 


ml/min 
cm 


CLASS B TEST 


AATCE spray 
Drop penetration 


Water penetration (drip 
Bundesmann 

Impact penetration 

Spray test 

Kern 

Franz and Henning 

Official German rain 
Spray penetration 


CLASS ¢ 


Immersion 


[36] 


ASTM water-absorption (spray) method [37] 


Dynamic absorption 
Bundes 


(3s) 80°+1° F 
(39) ? 

[40] Room temp 
(41) ? 


mann 
Becker submersion method 
Esslinger submersion method 


Very wide; AATCC; Federal Specifica- 
tion CCC-T-191a and Supplement. 
Federal Specification CCC-T 

Supplement 
ASTM 
British Specifications; Germany 


min 19la and 
xe” 
14 om 
6” 
6” x6" 

4cem? 

12X26 cm 
28X38 cm 
12” x12” U 


xR” 
Limited 
Do 
Germany 
Do 
8. Navy. 


of drops 


TEST 


sx’ 


g’?¢R”” 


Wide; AATCC, 
Federal Specification CCC-T 
Supplement 
Do 


19la and 


a’ xr” 

14 em? 
5x75 cm 
15x 


British Specifications; Germany 
Germany 


5em Do 


CLASS D TEST 


ingle by tensiometer 
test 


ethod 


[42] 
(43) 
[44] 


ition mark signifies temperature not indicated. 
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0.196 in. ? 


Experimental. 
Do. 
Do 


Single yarn 
ria 








purpose of this report to arrange the different 
test methods in four general classes: 


Class A. Methods by which the hydrostatic 
pressure required to force water 
through a fabric is measured. 

Class B. Methods by which surfece wetting and 
penetration under the influence of 
falling drops is measured. 

Class C. Methods by which the absorption of 
water by the fabric, when immersed 
or manipulated under water, is 
measured. 

Class D. Methods by which the wettability of 
the surface of the fabric is measured 
by means of the angle of contact or 
some function of the angle of contact. 


Some timely comments about the more popular 
test methods in each group are given below. 

Class A.—In this class of test methods, the 
fabric is subjected to the action of water under 
pressure by a variety of means. Either the 
amount of water penetrating in a specified time 
or the pressure required to force water through the 
fabric is measured. The most widely used test 
methods in this class are the AATCC [21] and the 
Suter [22] hydrostatic-pressure tests. Both meth- 
ods give reproducible test values. It has been 
reported [45] that a correlation exists between the 
two methods. The following relationship is sup- 
posed to hold at 27°+3° C: 

A=1.32S—4, (11) 
where A is the AATCC test value, S is the Suter 
test value. 

The hydrostatic-pressure test values are de- 
pendent mainly on the pore size and the angle of 
contact. The bulge that occurs when the pres- 
sure is applied on the fabric requires a correction 
|46], but in view of the reproducibility of the test 
and the smallness of the correction factor, the 
extra measurements and calculations are usually 
avoided. 

Class B.—The tests in this class are in some 
cases very different from each other. They are 
all, moreover, quite different from those in class 
A. All the tests in this class subject the test 
fabric to the action of water drops. The num- 
ber, size, frequency, and energy of the drops in 
the various tests vary considerably. The ad- 
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vantages and disadvantages of thes, 
tests are briefly discussed below: 

1. AATCC Spray Test. This test is very widy 
used in the textile industry for control work. 1), 
test will distinguish qualitatively between treaiy, 
and untreated fabrics. It is not able, howoey. 
to distinguish between obviously different wat. 
repellent finishes or to predict the performance of 
a fabric in the rain. 

2. Drop-PenetrRatTion Test.—This test yy 
developed for and used by the United Stay 
Quartermaster Corps during the war. The energy 
of the drops used in this apparatus is appmy. 
mately 10,000 to 15,000 ergs. This valu x 
much larger than that of the drops in a clow. 
burst (3,000 ergs). It appears that by this te 
one may be able to arrange fabrics in what seens 
to be a proper order of protection in the miy. 
The hardness of the fabric backing has a pp. 
found influence on the test value. This is prob. 
ably also true for the other tests in this grow 
It should be mentioned that the values obtained 
in the drop penetration test may differ as much ss 
20 percent among each other, depending on the 
particular fabric and its test value. 

3. BunpesMaANN Test.—The Bundesmann tet 
was developed in Germany. It can be used 
two ways: (1) to measure the amount of wate 
penetrating the sample, or (2) to measure the 
amount of water absorbed by the sample in a given 
time. Again the drops of water in the test have 
energy of from 15,000 to 30,000 ergs, which is 4 
much larger value than is found in a cloudburst 
In this procedure the sample of cloth is in constant 
motion and is continuously rubbed on the under 
side during the test. When the amount of water 
absorbed is measured, the mean deviation [33) o 
the test value is reported to be about 3 percen! 
The method seems to offer some advantages over 
the other tests in class B and has been recently 
adopted by the British Sub-Committee of th 
General Technical Committee as a standard te 
method. 

4. Impact Penerration Test.—This test & 
volves the spraying of 500 ml of water on the te! 
specimen. At the end of the spray period, 4 
blotter beneath the fabric is weighed. Th 
increase in weight of the blotter represents the 
amount of water that passed through the cloth 
No data are available on correlation with natu! 
rain. 


Various 
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Test.—No temperature control is 
ified in this test. It appears to be similar to 
impact penetration test. No data are avail- 
le for correlating the results of the test with 
ults obtained with fabrics in natural rain. 
; Kern Test.—This novel method has been 
It employs a single drop of water of 
An electric circuit is used to detect 
. penetration of the water to the under side of the 
t specimen. The average deviation of the test 
ue from the mean appears to be smaller for 
aller size drops. Although there appears to be 
remendous variation in the results, depending 
mn the fabric, the author states that seven tests 
. sufficient for obtaining “exact results.” 
(Class C.—All methods in this group measure the 
istance that the finish offers to wetting by 
ter. As droplet penetration is not used in these 
thods, the part played by the structure of the 
vic is less important than in other tests in 
termining the test value. The water absorbed 
the fabric is measured by weighing the test 
imen after some form of partial drying. In 
ne of the tests (dynamic absorption, Becker, 
the sample is in motion during the test 
riod of exposure to the water. In almost all 
es the deviation of any particular test value 
mm the mean appears to become smaller as the 
posure time gets larger. The dynamic types of 
ts are definitely to be preferred to the static 
pes in this category. Again, it should be pointed 
t that little, if any, work has been done to corre- 
e these test results with the performance of 
ries in the rain. 
Class D.—Workers who have used the following 
ree technics have done so in an attempt to 
derstand the mechanism of wetting. Various 
thods for measuring contact angles [47] have 
The three methods which have 
‘n applied to fabrics are discussed below. 
|. WenzeL Mernop.—Wenzel appears to have 
n the first to attempt to apply methods of 
asuring contact angles to fabrics [10]. He 
ulified the tilting-plate method [47] to make it 
plicable to fabrics. Wenzel measured the con- 
‘t angle of many materials, waxes, cellulose 


5. SPRAY 


n described. 


elate, many metal stearates, rosin, and others. 


sdata indicated that the method was applicab‘e 
textile surfaces and was very reproducible. At 
‘Bureau it was found possible to reproduce the 
lues that Wenzel obtained on paraffin wax on 
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glass. It was found, however, that this method 
lacked control over the following variables: (a) 
Temperature, (b) rate of flow of the water over 
the solid surface, (c) rate of rotation of the plate 
in the water, and (d) cleanliness of the surface 
of the liquid. 

2. Werrinc Test.—This test measures the 
time for the contact angle on a yarn or strip of 
fabric to decrease to 90 degrees. This time has 
been called the ‘‘wetting time.’’ Cassie and Baxter 
claim that the wetting test is more sensitive to 
proofing efficiency than either the Bundesmann 
or the hydrostatic-pressure tests. 

3. Contact ANGLE BY THE 'TENSIOMETER 
Mertuop.— Wakeham and Skau [42] found that, 
by modifying an ordinary interfacial tensiometer 
and forcing a circular piece of fabric through the 
surface of the water, they obtained a relationship 
between the pressure necessary to force the disk 
of cloth through the surface of the water and the 
angle of contact of the water to fabric. The 
method does not appear to be capable of measur- 
ing angles greater than 120 degrees. Unfortun- 
ately, the upper portion of the measurable range 
is the most useful from the point of view of water 
repellency. 

The great dependence of all test values on the 
temperature of the water has been mentioned in 
several reports [43, 45, 33, 48]. It is quite sur- 
prising, therefore, that many of the methods listed 
in table 1 appeared to be uncontrolled in tem- 
perature. 

The above review of the methods 
measure water repellency indicates that no one 
method completely measures the phenomenon. 
It appears necessary that a combination of tests 
be used to evaluate the water repellency of a 
fabric. It appears possible, also, to select one 
test from each class and report the results in the 
same units. For example, the hydrostatic pres- 
sure is increased at the rate of 1.0 cm a second in 
the AATCC test. The test value could be ex- 
pressed in time units. This could also be done 
with the drop penetration, the dynamic absorp- 
tion, the wetting and other tests. One could then 
obtain a summary test value which would be the 
results of n number of test values, all expressed in 
time units: 


used to 


Tf (e,6,¢,.... ®). (12) 


It seems reasonable to assume that the value de- 
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rived from n number of tests would be a more To what extent does dry cleaning gy) fabri 
reliable measure of the repellency of a fabric than laundering diminish the water repelleney of weve! 
taking the value of only one test. It would re- fabric? ay 
main to be proved, however, whether such a sum- What causes a fabric to lose its repelloney dy 
mary value gives an accurate and reliable indica- ing the course of exposure to the rain? | 
tion of the performance of a fabric in the rain. In order to have available a background of ini 
Considerable information obtained during the mation and data about water-repellent {aby p 
war indicates that many of the test methods and common test methods, 11 typical commer 80 
given in table 1 are not applicable for evaluating raincoat fabrics and 4 Army fabrics were selec 2 7 
the repellency of woolen fabrics [49]. All the for study. The characteristics of these 15 wats spe. 
test methods mentioned require further investi- repellent fabrics are given in table 2. The wath $ 40 
gation to see whether they can be adapted for repellent properties were examined by the folloy 30 
testing all types of water-repellent fabrics. It is ing test methods: AATCC hydrostatic-pressy . 
also necessary that some of these methods be test, AATCC spray test, drop-penetration ts P| 
studied for the purpose of setting up commercial and angle of contact test (tensiometer type 
and Government standard methods of testing. Figure 4 shows the results obtained by the {vy 
IV. Evaluation of Water Repellency by test mutate. = test was cons according j — 
Typical Test Methods the directions prescribed in Gov ernment specifi 5 100 
tions or published papers. It is noted that iM 3 s0 
In planning the laboratory work, it was neces- most widely used test method, the spray-ratigf ¥ %° 
sary to take into account the questions that are test, did not show any differences among 12 of i © be 
asked by Government agencies, industry, and the 15 fabrics. Great differences among the vari bs 
publie at large. The following general questions fabrics are revealed by the remaining three tes f «0 
are frequently asked: As each test measures different properties of (mg ® °° 
To what degree is a particular fabric water fabrics (see section on survey of methods, p. \i7 bio 
repellent? it is not surprising that the 3 tests do not rank (| ol 


TaBLe 2.— Properties of 15 water-repellent fabrics 


WwaTER 
— 


Number of plies Thread count Air permeabil- 





' ity at a pres- 8 50 
Sample No Type of weave sure of 0.5 in. Weight ; 
‘ . 7 . of water across 
Warp Filling Warp Filling fabric «3 
; 
WATER-REPELLENT COMMERCIAL FABRICS 8 30 
[*) 20 
(ft */min)/ft? Oz/yd? 5 
I Twill-tackle 1 ; 125 65 4 a6 "4 io} 
2 do I 3 125 67 17.9 &.9 ; 
3 Plain-poplin 2 l 107 51 5.0 5.8 z Py i 
4 do 2 1 109 52 6.2 5.9 | 
5 Twill 2 2 135 Ss 6.4 6.4 
" Plain-poplin 2 I 107 55 a0 §.7 r 
7 do 2 1 107 M4 5.0 6.0 130 
8 T will-eabardime 2 1 127 M 16.2 6.5 
9 Plain-poplin 2 I 106 MM 4.6 6.3 8 0 
10 Twill i 1 104 60 23.6 7.3 g 
& 0 
i Plain-poplin 2 1 107 A) 5.3 6.2 ' | 
_— 
= —L= a tod j 00 
ARMY WATER-REPELLENT FABRICS zo 
v 9 
ethiatiaitittieninn aaiel Fee ES hnintall 9 
~_ 
z 
12 Sateen 2 2 115 70 3.9 9.1 8 80 
15 Oxford i 1 134 53 3.2 62 | 
| 
a Sateen 2 2 115 67 3.8 8.9 | 70 
» Poplin " 1 53 52 4.1 62 | 
| Figure 
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fabrics in exactly the same order. It is noted, 
wever, that all fabrics can be approximately 
ted into 3 general groups: good, fair, and poor. 
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Ficure 4.—Comparison of four water-repellency tests. 
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Certain samples always appear in a specific group: 
e. g. No. 5, poor; No. 3, good; No. 6, fair, ete. 
The fallacy of using only one test method to 
evaluate water repellency can be clearly seen from 
the data presented. It would be quite impossible 
to classify a particular water-repellent fabric on 
the basis of only one test. 

Both the hydrostatic-pressure test and the 
measurements-on contact angles reveal certain 
differences among the various samples not re- 
vealed by the spray test. The drop-penetration 
test, however, appears to be able to show differ- 
ences in the fabrics that were not revealed by any 
of the other tests. The limitations of a par- 
ticular laboratory test for evaluating water re- 
pellence are generally recognized. The most 
frequent argument against a certain test is that it 
does not test the fabric under conditions of actual 
use. Tests that simulate natural rain, therefore, 
are considered much more useful for practical 
testing. The drop-penetration test appears to 
be very useful in predicting the relative perform- 
ance of fabrics in the rain. It seemed desirable, 
therefore, to study the present group of fabrics in 
greater detail on the basis of this test. In view 
of the arbitrariness of the test value—the time 
to collect 10 ml of water passing through the test 
specimen—it was decided to study the rate of 
penetration of water through water-repellent 
fabrics for long periods of exposure. It was found 
that the “rate of penetration (ml/min) versus 
time curve” was very characteristic of the nature 
of the particular fabric. Figures 5 and 6 show 
plots of the rate of penetration of water against 
time of exposure, for 11 of the fabrics studied. 
The arrows in figures 5 and 6 indicate the time 
at which 10 ml of water had passed through the 
fabric. Each rate curve contains the characteris- 
tic portions: an induction period, an S-shaped 
portion of a curve, and a limiting rate region. 
It is believed that the penetration of water through 
a fabric in natural rain follows a similar pattern 
and that the curves obtained for the drop-penetra- 
tion apparatus are not a characteristic of the 
apparatus alone. These curves reveal that during 
the induction period little or no water comes 
through the cloth. When the resistance to pene- 
tration is finally broken, water passes through the 
fabric at a very rapid rate (S-shaped region). 
Finally, the amount of water that penetrates the 
fabric per unit time reaches a steady state. As 
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ost of the water-repellent fabrics tested are at 
st only somewhat resistant to penetration, it is 
oar that, after prolonged exposure of a water- 
pellent garment to the rain, the wearer would 
1t keep dry. The curves also indicate that the 
nger the induction period, the lower will be the 
yilibrium rate of penetration of water. In 
signing new water-repellent fabrics, therefore, 
e object would be to fix the various variables 
volved so that the resulting product would have 
long induction period and a low rate of pene- 
ation. Such a fabric has been designed. The 
irley type of cloths that the Army Quarter- 
aster Corps developed during the war have a 
ry long induction period (several hours) and a 
ry low equilibrium rate of penetration (approxi- 
bately 0.1 ml/min). These types of fabrics are 
ivhly effective for prolonged use in tbe rain. 

The above studies of water repellency lead to 
he conclusion that obtaining the rate of penetra- 
m of water with the drop-penetration apparatus 
ves a more complete picture of the repellency of 
fabric than the 10-ml test value by itself. This 
especially significant when one is considering 
mmercial fabrics (compare fig. 5 with 6). It 
also concluded that better evaluation of water 
pellency may be made when a group of tests is 
mployed. A study of the four test methods 
ven above reveals that, as the rank correlation 
nong various tests is very low, it is necessary to 
relate all tests with the behavior of fabrics in 
piural or artificial rain. The Philadelphia Quar- 
prmaster Depot has developed and constructed 


[50] a rain room for fabric and garment testing. 
By means of specially designed nozzles, it is 
possible to obtain rainfalls at the rate of 0.1, 1.0, 
and 3.0 inches an hour, having natural and re- 
producible drop-size distributions. 


V. Durability of Water-Repellent Finishes 


The durability of water-repellent finishes is 
generally taken to be the resistance that the 
finish shows to the action of dry cleaning or 
laundering. Durable water-repellent fabrics are 
those whose water repellency is only slightly 
diminished by dry cleaning or laundering. Those 
fabrics whose water repellency is totally or greatly 
impaired by dry cleaning or laundering are classi- 
fied as being nondurable. It should be remembered 
that the two terms are only relative and that the 
difference between them is a matter of degree of 
resistance and not necessarily of type of compound. 

The effect of dry cleaning and laundering on 
water repellency was studied, using only the 11 
commercial raincoat fabrics given in table 2. 
Standard methods of dry cleaning and laundering 
were used [51]. The results of these studies are 
given in table 3. By the use of three different 
test methods, it was found possible to distinguish 
between durable and nondurable repellent finishes. 
Samples 4, 7, 9, and 11 show a smaller decrease 
in all test values than do the remaining samples 
which reveal fairly large changes in all test values. 
It was found that samples whose water repellency 
had been totally lost by dry cleaning or laundering, 


TABLE 3.—Effect of dry cleaning and laundering on the water repellency of 11 commercial fabrics 


Spray rating 


Sample No 
Dry 
cleaned 
three times 


Laundered 


Control three times 


100 


9 





| 100 
| 


Control 


Hydrostatic pressure (centimeters | Drop-penetration time (time, in 


of water) minutes, to collect 10 ml of water) 


_ 
ars Laundered 
three times three times 


Laundered 


» 1 
three times Contro! 


Dry 
cleaned 
three times 





* Droy \etration time could not be measured as the swatches became completely saturated with water at the start of the test. 
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gave no test value in the drop-penetration test. 
These samples became completely saturated with 
water as soon as they were placed in the test 
apparatus. 

The lowering of the water repellency of fabrics 
by dry cleaning or laundering is commonly 
attributed to leaching of the repellent compound. 
The presence of residual soap on the fabric and 
distortions of the fabric structure are also con- 
sidered responsible for the loss in repellency. It 
has been observed, however, that simple wetting 
by water will also decrease the repellency of a 
fabric. Tests done at the Philadelphia Quarter- 
master Depot [52] have shown that water-repel- 
lent jackets were more easily wetted if they had 
been previously exposed to the rain. It was found 
that drying the jackets after the test exposure, 
even by heating in an oven, failed to restore to the 
garments their original water repellency. In view 
of the interest in the durability of water repel- 
lency, it seemed advisable to attempt to establish 
the rate at which some fabrics lose their repellency 
when exposed to the rain. For this study, four 
swatches each of five different water-repellent 
fabrics were tested five different days in the drop- 
penetration test. At the end of each test period, 
the swatches were allowed to dry at room tem- 
perature and then were conditioned overnight at 














NUMBER OF TIMES TESTED 


Ficure 7.— Lowering of the drop-penetration time owing 
to repeated testing. 
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65-percent relative humidity and 70° F. 4, 
results of this experiment are shown in figur ; 
It was found that the various fabrics lose thy 
repellency at different rates. The drop-pene 
tion test values get progressively lower cach tin. 
the swatches are tested. This loss in repelley 
appears to be an irreversible process. Attemp: 
to restore the original repellency by heating } 
vacuum at elevated temperatures were Unsuccos 
ful. An analysis of the factors that could ) 
causing the loss in repellency leads to the { 
lowing: Loss in the repellent agent by leaching 
creation of new uncoated surfaces due to swellizg 
of fibers, and changes in the surface geometry; 
the fabric due to the pounding of the high eney 
drops. Any or all of these changes could \y 
taking place, and they all would be irreversitk 
changes. 

In order to throw light on the loss in repellency 
two simple exploratory experiments were mai 
The changes in the contact angle of two of t 
fabrics were measured after various periods 
wetting in the drop-penetration apparatus. Aly 
the change in contact angle and the increas: s 
“standard moisture regain” * of one sample wey 
noted after prolonged periods of spraying. Tl 
spraying was done with the nozzle of the AAT( 
spray test and was continuous for the partici 
period of time. The results of these observation 
are shown in figures 8 and 9. The contact ang 
of the dry fabrics was found to be lower al 
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Ficure 8.— Change in contact angle of fabrics expos 
the drop-penetration apparatus. 


Journal of Researiggy ater 


delleney 
e made 


articulay 


rvation 
ct angi 


‘er altel 


ch period of wetting in the drop-penetration 

The contact angle of one of the fabrics was 
.» lowered after spraying the sample for various 
es. Swatches of this fabric, when sprayed 
nd then tested for standard moisture regain, 


showed an increase in rate of moisture absorption. 
It is emphasized that the increase is an increase in 
rate (per 24-hour period) as the “repellent” 
fabrics do not always reach equilibrium in this 
length of time. 
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EXPOSURE TO WATER SPRAY—HOURS 


cure 9.— Changes in the contact angle and the standard moisture regain of a 9-ounce water-repellent fabric after prolonged 
spraying with water. 


VI. Summary and Conclusions 


The theory of water repellency of textile fabrics 
s been reviewed with special references to the 
hore recent theories of the wetting of fabrics by 
ater. A survey has also been made on the vari- 
s testing methods that have been devised 
br measuring water repellency. The results 
the present investigation lead to the following 
nelusions regarding the status of water repel- 
ucy. There is definite need for a comprehensive 
udy of the role that the structure of a fabric 


lays in the phenomena of water repellency. In 


e past, the emphasis has been on developing 
ore efficient compounds. Data available indi- 
te that a better understanding of fabric construc- 
m as it applies to repellency, coupled with the 
w available water-repellent agents, will lead to 
me more nearly idealized type of water-repellent 
hrment. In regard to testing methods, it is re- 
lired that correlation be established between the 
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results of laboratory test methods and performance 
of fabrics in the rain. 

As already stated, the contact angle is influenced 
by the following factors: The chemical nature of 
the solid surface, the porosity of the surface, and 
the presence of other molecules on the surface. 
Again, any one or all of these factors could dimin- 
ish the contact angle during wetting of the fabric. 
The change in moisture regain of sample 24 
(fig. 9) shows that the rate at which a fabric 
absorbs moisture increases in proportion to the 
number of times the sample has been wetted. It 
is of interest to examine the surface factors that 
might be responsible for the increase in moisture 
absorption. The loss in repellent agent, the 
change in position of the fibers in the yarns, and 
the creation of new surfaces could all affect the 
rate of water absorption. The swelling of par- 
tially coated or uncoated fibers would also result 
in making available more hydrophilic surfaces 
(OH groups) [53]. 
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dsorption of Water Vapor by Untanned Hide and 
Various Leathers at 100° F 


By Joseph R. Kanagy 


Adsorption of water vapor by untanned hide and various tannages of leather was 
The Brunauer, 


. 


determined at relative humidities ranging from 0 to 96 percent at 100° F. 
Emmett, and Teller equations for adsorption of gases in multimolecular layers were applied 
Adsorption for all of the samples below a relative humidity of 75 percent is a 
Above 75-percent relative humidity, other factors exert 


to the data. 
function of the available surface. 
an influence, the most important of which is probably size of the capillaries. Untanned 
hide and chrome-tanned leather adsorb the most water vapor except at relative humidities 
above 75 percent, at which sole leather adsorbs more. The behavior of the sole leather in 
this region is ascribed to the presence of deliquescent materials used in the finishing process. 
The other vegetable-tanned leathers adsorb less moisture than either hide or chrome 
Calculated heats of adsorption indicate that the affinity of groups for water is not 
Methods by 


which the data may be applied to practical processes, such as fat-liquoring and tanning, are 


leather. 
appreciably altered by tanning and that the adsorption is physical in nature. 


pointed out. 


I. Introduction ' 


The amount of moisture that leather adsorbs or 
determined by temperature, relative 
umidity, degree of porosity, and the size of the 
Moisture has great practical significance 
because its amount affects the durability of leather, 
nd in articles such as shoes, gloves, and other 
arments, the comfort of the wearer. High mois- 
ure content accelerates deterioration and pro- 
ites mildew action. On the other hand, a 
unimum amount of moisture is required to keep 
ather properly lubricated and thus prevent 
racking. The adsorption of moisture affects 
omfort by its influence on weight, area, and 
isulation. Insulation is a function of the amount 
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of water adsorbed, water vapor transmission, rate 
of condensation, and evaporation. 

From data obtained in the study of the moisture 
adsorption by leather and untanned hide at various 
relative humidities, theoretical concepts as to the 
nature and magnitude of the surfaces may be 
obtained. Such concepts may be of value in 
determining the nature of tanning and in studying 
other processes, such as fatliquoring and dyeing. 

The effect of relative humidity on strength and 
stretch of leather has long been known, and it is 
now customary to make all such tests under con- 
trolled conditions of temperature and humidity. 
Veitch, Frey, and Leinbach [1]? studied the 
strength and stretch of vegetable-tanned calf 
leather at relative humidities of 35, 55, and 75 
percent. Wilson and Kern [2, 3] determined the 


? Figures in brackets indicate the literature references at the end of this 
paper. 
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variations in strength, stretch and area for both 
vegetable and chrome-tanned leathers. Evans 
and Critchfield [4] made an investigation of the 
amounts of moisture adsorbed by chrome and 
vegetable-tanned calf leathers at different per- 
centages of relative humidity at a temperature of 
21° C. Recently Green [5] determined the 
moisture content of chrome leather under dif- 
ferent percentages of relative humidity at 25° C. 
These studies have all been somewhat limited in 
scope and, in general, the results were obtained 
where static methods for obtaining the relative 
humidity conditions were used. 

Work on the determination of the internal 
surface area of hide and formaldehyde-tanned hide 
has recently been done by Zettlemoyer, Schweitzer, 
and Walker [6]. These determinations were made 
by studying the adsorption of nitrogen. Bull [7] 
studied the adsorption of water vapor by a 
number of different proteins, one of which was 
collagen. 

This paper gives results of an investigation on 
the moisture relations of collagen and of leather. 
This is part of a study that is intended to cover a 
wide range of temperatures and relative humidi- 
ties. The results given were obtained with col- 
lagen (commercial hide powder) and various tan- 
nages of leather at 100° F and at relative humidi- 
ties ranging from 0 to 96 percent. 


II. Description of Samples Used 


The collagen used was a sample of commercially 
prepared hide powder. A more highly purified 
collagen has been prepared since this work was 
started and will be compared with the commercial 
hide powder in future work. 

Six different samples of leather were used. As 
there was no precedent to determine the exact 
nature of the data that would be obtained, the 
selection of samples was made somewhat arbi- 
trarily, but included leathers which differed from 
each other in properties and represented some of 
those most important commercially. Sole leather 
was used because of its importance and the fact 
that the addition of salts and glucose in the 
finishing exerts an influence on the moisture 
adsorbed. This was compared with a similar 
tannage of leather from which the grease had been 
removed, and glucose and salts were removed by 
washing. Water-soluble material in the leather is 
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mostly vegetable tannin, which it is almos; in 
possible to remove completely by washiss 
Chrome hydraulic-packing leather was use ; 
represent a mineral tannage. Chrome-tetan 
leather represents the combination of the miyy. 
and vegetable tannage. Calfskin and beliy 
leather were selected to represent two types, 
finished vegetable-tanned leathers. The analy 

of the leathers are given in table 1. 


TABLE 1.—Chemical analysis of the leathers used 
water-vapor-adsorpltion measurements, based on moix 
free weight 


2 De- 
Water Hide 
aa — greeof Acid P 


Leather _ 
bles | stance ii.) 


Per- Per- 
cent cent 
Sole 5 36.5 33.0 
Sole, degreased 
and washed 10.8 
Chrome, hydraulic- 
packing 
Chrome-retanned 
upper, degreased 
Calfskin, vege- 
table-tanned 
Belting 


III. Methods of Test and Procedure 


The conditioning equipment was developed | 

Carson and Worthington for use in studying 
water-vapor permeabilities of sheet materials 
This equipment is shown in figures 1 and 2, ands 
detailed description is given in reference {8} 
It consists essentially of a sealed temperatur- 
controlled chamber in which a constant relativ 
humidity is maintained by means of the proper 
saturated salt solution. Air is circulated over thy 
salt solution and throughout the chamber cov 
tinually. It is equipped with a mechanism for 
weizhing the samples in the chamber. This s 
done by means of the rotating disk (36), figure? 
from which the samples (1) are suspended. Each 
sample is weighed by rotating it to the positi’ 
of the weighing arm (40), and transferring it from 
the disk to the extended arm of the balance \y 
means of a device on top of the chamber. 

The water-adsorption measurements were mat 
with leather samples ground to pass a 4-mm sie'' 
Weighed portions of the samples were placed 0 
small cylindrical baskets (approximately 2 in. 
diameter by 2% in. high, shown in fig. 3) made 
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The same sample of each particular material was 
J to obtain all the data in the adsorption and 
sorption procedure. The samples were allowed 


s come to equilibrium at a constant relative 
midity at 100° F. The salt solution was then 
banged, and the same samples were brought to 
builibrium at a slightly higher relative humidity. 
his procedure was repeated at 10 different points 
m 0 to 96 percent. The desorption determina- 
ons were made by reversing this procedure by 
arting at 96-percent relative humidity and end- 
vat about 0 percent. The salts used for main- 
pining the constant relative humidities and their 
alues as determined by the wet-dry bulb, and by 
irect vapor-pressure determination [9] are shown 














» table 2. 
To correct for the moisture adsorbed by the 
mtainer, an empty basket was used as a blank 
nd weighed each time weighings were made with 
Figure 1.—Egquipment used in obtaining data on adsorption 


of water vapor by untanned hide and various leathers at 
100° F. 
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various leathers at 100° F. 
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Fiaure 3. 


TABLE 2. 


Wire basket used to contain the leather samples 
used in the study of the amount of water vapor adsorbed at 
100° F. 


Relative humidities used in obtaining the results 
as determined by the wet-bulb method and the relative 


vapor pressure determined at the Institute of Paper 


Chemistry (9) 
Method 


Relative vapor- 


Wet-bulb pressure 


MgCl0, (Anhydroné 0 (approx.) 
Licl.. 12. 
KCyHyO, » 
MeCh.6H,O $1 

Ke Oy 43 

NaBr 53. 5 
NaNO; 61. i 

NaCl 74 

KeCrO, M.: 

K»S0, 04. 


the samples. The increase in the weight 0; 
blank was subtracted from the increase jp 4 
weight of the basket and the sample. 4, 
amount of moisture adsorbed by the basket , 
quite low and amounted to a total of only abo: 
0.35 percent for a 5-g sample of leather fo; ‘ 
relative humidities up to 96 percent. At » 
percent relative humidity, an additional! cor, 
tion of about 0.45 percent was necessary for ij) 
same weight of leather. | 

No difficulty was experienced in weighing \ 
samples at any condition of relative humidiy 
except at 96 percent. At this humidity, moisty, 
condensed on the rod leading from the balay 
arm to the sample. It was, however, possibj 
to make a correction for this so that the error) 
the determinations was minimized. 


IV. Results of Measurements 
1. “Dry’’ Weight of the Samples 


The base or “dry” weights of the samples use 
in the experiments were obtained from the equil- 
brium values in the presence of magenesium per. 
chlorate Anhydrone. A_ relative humidity 
about 10 (dry-wet bulb reading) was obtaine 
with this material by first permitting the attai- 


ment of equilibrium with calcium chloride. Ths 
brought the humidity to about 2 percent. The 
with the addition of magnesium perchlorat 
was possible to obtain a velue of approximate) 
for the humidity for 3 days. 

The customary method for determining moist: 
or obtaining a dry sample of leather is to heat |! 


material at 100° C for 18 hours. However, as ths 
treatment is known to alter the structure of leather 
and collagen somewhat, these materials migh! 
have been rendered unfit for use at lower temper 
atures. On completion of the tests, the sample 
were placed in an oven at 100° C for 18 hours 
determine the difference between the values 0 
tained after reaching equilibrium in the presen 
of Anhydrone and after heating. These resus 
are given in table 3. 

All the samples lost additional weight on bei: 
heated. The additional losses were all less tha’ 
1 percent and in four of the seven samples, the) 


were under 0.5 percent. 
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Untanne 
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Chrome | 
Chrome-! 

degress 
Sole lea 
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umidit: 
to 7-« 
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The 


COrTe 


JOSS! Die 


TTOr ip 


Lained 
ttain 


Difference in weights of samples in equilibrium 
ce of Anhydrone and after heating 18 hours at 
» dry-air stream 


Weights of samples 


Per- 
cent 


Differ- 
ence 


After 
heatin 
at 100° C 
in dry air 
stream 


After 
equilibrium) 
in presence 
of Anhy- 
drone 


Sample 


v v 
Untanned hide 6. 0226 5. 9928 
Sole leather : wi 4. 3611 4.3199 | 412 
Chrome leather | 3.27% . 0249 
Chrome-retanned leather, 
4. 7497 . 0231 | 


0. 66 
"4 
. 76 


. 
| 0.0398 


degreased am . 
Sole leather, degreased, | 
0061 
. 0096 
. 0135 


5. 0937 
6. 2011 
5. 4971 


washed ° » 
Vegetable-tanned calfskin 
Belting leather 


6. 2107 
5. 5106 


2. Adsorption Isotherms 


The time required to reach equilibrium at a 
finite relative humidity is shown in figure 4, 
here the rate at which equilibrium was reached 
11- and 54-percent relative humidity is given. 
quilibrium is reached in 1 day at a relative 
umidity of 11 percent, whereas at 54 percent a 
to 7-day period was required. 
Adsorption isotherms at 100° F for six different 
athers and for untanned hide are given in figure 
The curves are all of the sigmoid type, with a 
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TIME IN OAaYS 


Figure 4.—Rate at which equilibrium is reached at two 
different relative humidities by hide substance, chrome 
leather and belting leather when exposed in an atmosphere 
of water vapor at 100° F. 


O=hide substance; X =chrome leather; @= belting leather. 


fairly rapid initial rise, then a leveling off, and 
finally at a relative humidity of 75 percent, again 
a rapid rise. The very rapid rise shown by sole 
leather above 75-percent relative humidity is ap- 
parently caused by the magnesium sulfate and 
glucose it contains. The sample of sole leather 
which was degreased and washed shows no such 
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Figure 5.—Increase in the amount of woter vapor adsorbed 
untanned hide 


Adsorption of Water Vapor by Hide and Leather 


with increase in relative humidity for various leathers and 
at 100° F. 
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increase in the amount of moisture adsorbed. 
Untanned hide and chrome leather adsorb almost 
identical amounts of moisture at all relative 
humidities. 


3. Desorption Isotherms 


Adsorption and desorption curves for hide 
powder and belting leather are given in figure 6. 
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Figure 6.—Adsorpticn and desorption of water vapor by 
hide powder and belting leather at various relative humidities 
at 100° F. 


X =hide powder; O= belting leather. 


The hysteresis for untanned hide is no greater 
than that for leather, even though the amount of 
moisture adsorbed by the former is greater. The 
steep portion of the curve for collagen at high 
relative humidities is more nearly reversible than 
the flat part of the curve at the lower relative 
humidities. This may be explained by assuming 
that at high values of relative humidity, the 
adsorption consists mainly in the filling of the 
larger pores, which is reversible, because con- 
densation of water on water molecules rather 
than water molecules on leather surface is in- 
volved. At the lower relative humidities adsorp- 
tion consists in condensation where leather sur- 
face is involved, and this is not reversible because 
the surface area has increased because of swelling 
under the conditions of high relative humidity. 
All samples give off practically the same 
amounts of moisture on desorption as was gained 
on adsorption; the curves coincide at 0 relative 
humidity as shown for untanned bide and belting 
leather in figure 6. Sole leather (table 4) shows a 
loss of nearly 3 percent more weight on desorp- 
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tion than was gained on adsorption which mip 
be due to loss of salts or glucose. The weigh 
gained on adsorption and lost on desorption ; 
all the samples is given in table 4. 


TaBie 4.—Amounts of moisture adsorted by sample , 
adsorption and losi on desorption 


Adsorption Desorpticy 
moist ure Moistur 
gained lost 


Percent 
Collagen (hide powder) 55.9 
Sole leather 83.2 
Chrome leather 55.0 
Chrome-retanned leather, degreased 33.1 
Sole leather, degreased, washed 2.7 
Vegetable-tanned calf skin ‘ 25.0 
Belting leather 27.9 


Percent 


V. Calculation of Derived Quantities 


Application of Brunauer, Emmett, and Tella 
Equations 


Brunauer, Emmett, and Teller [10] have deve 
oped a theory to account for the adsorption 
gases in multimolecular layers. This theory wis 
developed by a generalization of Langmuir’s treat- 
ment of the unimolecular layer and has recently 
been used to calculate surface areas by a number 
of authors [7, 11). 

Two equations are developed {10] that may be 
applied to the data of this paper. 


Pp _t=1 2 
V(po—p) —(C Va Po 
V,,cx 1— (n+1)a"+nz"*" 
nti’? 


V= 

I—z 1+(e—1)r—er 

V=volume of water vapor adsorbed pe 
gram at 0° and 760 mm 


V,.=Volume of water vapor at 0° and 7i) 
mm, which, when adsorbed, coves 
the surface of 1 g of the adsorbent 
with a unimolecular layer 


p= pressure 


z=relative pressure of the water vap0 
that is equal to p/po, where py is the 
saturation pressure of the wale 
vapor at the temperature of the 
adsorption 
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a measure of the increase of the heat of 
adsorption £, of the first layer on 
the adsorbent over the heat of 
liquefaction (£,), according to the 
equation E,—E,=2.303RT log e. 


. 


From equation 1, which holds only at values of 
p, less than 0.5, a straight line may be obtained 

plotting experimental values of p/V(po—p) as 
dinate against p/po as abscissa. The intercept 
ll be 1/V,,¢ and the slope (c—1)/V,¢. From the 
perimental values of intercept and slope thus 
tained, c and V, may be calculated. Typical 
lots of this type are shown in figure 7, and the 
slues of ce and V,, for hide substance and all the 
athers are given in table 5. 


inte 5.—Values of surface properties for hide powder 
and various leathers calculated from the Brunauer, Em- 
mett, and Teller equations 


Sample 


ml/g 

110.0 | 
86.7 | 
110.1 | 
97.2 | 
89.0 
76.3 | 
83.3 | 


Collagen (hide powder) 

Sole leather 

Chrome leather . 
Chrome-retanned leather (degreased) 


8.7 263 
241 
206 


225 


Sole leather (degreased and washed) 
Vegetable-tanned calfskin 
Vegetable-tanned belting leather. 


8.2 
8.8) 


V.. is the volume of water vapor at 0°C and 760 mm that is required to 
ms unimolecular layer on the surface of the sample 
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IGURE 7 Linear plot of p/v (Po— P) against P/Pe for the 
adsorption of water vapor on vegetable-tanned calfekin ond 
hide powder at 100° F. 


0.3 0.4 


X = vegetable calfskin; O = hide powder. 
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VOLUME OF WATER VAPOR ADSORBED (ML/G AT STP) 
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P/Po 
Figure 8.—Observed and calculated isotherms for the ad- 
sorption of water vapor on untanned hide at 100° F. 
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The values for ¢ and V,, may be substituted into 
equation 2, and V calculated for different values 
of n, where n equals the number of adsorbed layers. 
These values may then be compared with the de- 
termined results for the volume of water vapor 
adsorbed to ascertain the number of adsorbed 
layers which are equivalent to them. Such re- 
sults are given in figure 8 for untanned hide, and 
in figure 9 for leather. The amounts of moisture 
adsorbed by leather at high relative humidities 
appear to be equivalent to five or six layers, where- 
as with untanned hide a large number of layers 
appears to be formed at these pressures, which 
probably indicates that the capillaries in the latter 
are much larger than in the former. 


2. Calculation of Surface Areas 


With the values for V,, (milliliters of waves 
vapor adsorbed per gram of sample) it is possible 
to calculate the surface areas of the materials 
if the volume occupied by each molecule of water 
is known. Taking the value determined by 
Smith and Green [12], namely, 10 square angstrom 
units as the area occupied by the water molecule 
and 2.70510", the number of molecules per 
milliliter at 0° and 760 mm, the surface area per 
gram of material is calculated as follows: 


(VX 2.705 X 10" 10 10-™)m?=2.705 Vp. 


The values for the surface area of all of the samples 


in square meters per gram are given in table 5. 
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VOLUME OF WATER VAPOR (ML/G AT STP) 


Ficure 9. 


Collagen and chrome-tanned leathers have the 
largest surface areas. The leathers tanned with 
the vegetable tanning materials have smaller 
surface areas, as would be expected, since they 
are composed of less hide substance and the capil- 
laries are reduced to smaller diameters, in some 
cases probably completely filled by tanning mate- 
rials. The surface areas of hide and the leathers 
are lower than that found in general for charcoals. 
This might be due to the granular nature of ground 
charcoal. The value of 298 m*/g of collagen is 
considered to be in fairly good agreement with 
the value of 337 m*/g obtained by Bull [7]. 


3. Heats of Adsorption 


By substituting the values for ¢ given in table 
4 in the equation L,—E,=2.303R T log e, it is 
possible to calculate the difference between the 
heat of adsorption and the heat of liquefaction of 
water on collagen and leather. The values for 
FE, — EF, are given in table 5. The positive value 
for £,—E, in contrast to the negative value 
obtained when water is adsorbed on charcoal [13] 
gives evidence of the effect of polar groups on 
the attraction of water. The fact that the magni- 
tude of the value obtained for collagen, 1,307 
cal/mole, is comparable with that obtained for 
the leathers, 1,284 to 1,434 cal/mole, may indicate 
that there is little alteration of groups that attract 
water when tanning occurs. The small difference 
between the heats of adsorption and heats of 


126 


P/P. 


Observed and calculated isotherms for the adsorption of water vapor on vegetable-tanned calfskin at 10/ | 


liquefaction, calculated by applying the B. £.1 
theory, suggests strong Van der Waals adsorptin 
Chemisorption is generally accompanied by high 
heats of adsorption (20,000 to 50,000 cal/mok 
[13]. Further work is planned in which calcu 
tions of the differential heats of adsorption will | 
made. 


VI. Discussion 
1. Reduction to Unit Surface 


If the volumes of water vapor adsorbed at each 
relative humidity are divided by the calculated 
values of V,, given in table 5, a factor is obtained 
that is nearly identical for all the samples at eae 
value of p/p) up to 0.75. These results are show 
infigure 10. Above a p/po value of 0.75, a defini 
variation occurs among the samples for the valu 
of this factor. Sole leather assumes the valu 
represented by the upper curve, collagen «0 
chrome-tanned leather are represented by th 
middle curve, and four samples of vegetable 
tanned leathers have the values represented \y 
the lower curve. 

These results indicate that up to 75-percet! 
relative humidity the adsorption is a functw 
of surface area alone. Above this value oil 
factors begin to influence the results. The mo 
important of these is probably differences in the 
sizes of the pores. The rapid increase for se 
leather is undoubtedly due to the deliquescet! 
nature of the glucose and magnesium sulfate 
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icure 10,—Adsorption isotherms of water vapor on various 

cathers and untanned hide at 100° F 

surface ina plot of V/Vm against p/ po. 


reduced to unit 


sed in the finishing process. The difference in 
he behavior of collagen and chrome-tanned 
ather, on the one hand, and the four vegetable 
annages, on the other, can be explained by as- 
uming that the former have larger pores on the 
verage and consequently adsorb more moisture 
t high relative humidities. 


2. Diameter of the Collagen Fiber 


As the surface area and the density of collagen 
re known, it is possible to calculate the diameter 
yf the fiber if it is assumed that it is a uniform 
viinder and that the calculated adsorption area 

identical with the total surface of collagen 
bbers. The areas of the ends of cylinder are 
begligible. 

As a cylinder of radius r and length Z has a 
urfaee S=2rrL and a volume V=-arL, the 
ameter 2r=4V/S. If the density of collagen 

taken as 1.4 g/em® (specific volume 0.7 em*/g), 

‘in the above equation is 0.7 em*/g and S (table 
Therefore, 


is 298 x 10* em?/g. 


<0.7 
2r me joe 94 X 10-8 em. 

Bear [14] found from X-ray data a distance of 
bout 11 A between the polypeptide chains in 
ollagen, which is equivalent to about one-eighth 
he dimension calculated here. Bear also found 

separation of polypeptide chains when he ob- 
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served the X-ray pattern of wet collagen. Ap- 
parently water goes between the chains and 
spreads them apart. These results indicate that 
moisture does not go between all of the chains at 
the low relative humidities at which the B. E. T. 
calculations were made. It is probable that this 
occurs at the higher relative humidities. 


3. Unimolecular Layer of Oil in Leather 


From the surface area of leather it is interesting 
to make a calculation as to the amount of oil that 
would be required to cover all the surface with a 
unimolecular film. Leather oils are composed of 
a mixture of triolein and tristearin and a number of 
other constituents that are decomposition products 
of these two compounds. 

From Langmuir’s data [15] on the molecular 
sizes of these compounds, it is possible to make a 
calculation as to the amount that would be re- 
quired for a unimolecular layer by using the data 
for washed and degreased sole leather. If the 
orientation of the chains is assumed to be per- 
pendicular to surface, with tristearin, 53.7 percent 
would be required for a single layer; with triolein, 
27.7 percent would be required. If it is assumed 
that the molecules lie flat on the leather surface, 
with tristearin, 17.7 percent would be necessary ; 
with triolein 24.4 percent would be required. 
These results are obviously not in harmony with 
practical experience. The reason is that surface 
available to water is not available to oils, as the 
oil molecules are 6 to 10 times larger in diameter. 


4. Relation of Adsorption to Capillarity 


It has been shown that the adsorption of water 
vapor by leather and hide shows correlation with 
the multimolecular adsorption theory of Brunauer, 
Emmett, and Teller for relative humidities below 
50 percent. The adsorption may also be con- 
sidered especially at the higher relative humidities 
to be the result of capillary condensation. The 
surface of leather may be considered as consisting 
of numerous capillaries of different sizes. This 
concept arises from its fibrous structure. Be- 
tween each fibril water may be adsorbed. The 
openings of various sizes between fibers give, in 
effect, the same characteristics as capillaries of 
various sizes. In the first stages of adsorption, 
the small capillaries will be filled, whereas the last 
stages will consist in filling the larger pores. 
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Calculations may be made in accordance with 
the Kelvin equation, which applies to the conden- 
sation of liquids in cylindrical capillaries. 


In (p/po) 2ye/rRT, 
in which 
Y=surface tension 
e—molar volume of the liquid 
r==radius of the capillary. 


A plot showing the amounts of liquid water ad- 
sorbed per gram of sample at values of p/po of 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9, in relation to 
the values for the radius calculated by substituting 
these values for p/pp in the Kelvin equation is 
shown in figure 11. 
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Ficure 11.—Adsorption of water vapor at 100° C on un™ 
tanned hide and two samples of leather in a plot of the 
volume of liquid water, V;, per gram of adsorbent against 
the corresponding radius of a cylindrical capillary as 
determined from the Kelvin equation. 


(The points for each sample correspond to p/pe=0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 
consecutively, reading from left to right.) O=hide powder; @==sole leather 
degreased and washed; X=vegetable-tanned calfskin. 


This figure gives a curve that represents a dis- 
tribution of capillary sizes in collagen and leather. 
According to this interpretation, nearly half the 
water adsorbed by leather below a p/p value of 
0.9 could take place in capillaries less than 10 A 
in radius, whereas about one-third of the water 
adsorbed by collagen could be taken up in capil- 
laries of this size. It can be seen from these re- 
sults that much of the surface available to water 





would not be available to oils having n Oleculgy 
cross sections of 6 to 10 A. Orientation coy 
occur in only about half of the capillaries availa), 
to water vapor, and therefore the amount » 
grease adsorbed would be greatly reduced » 
account of plugging of the pores. This interpy. 
tation is more nearly in line with what is found 
actual practice. 

An application with respect to particle sp 
might be made, in the light of the above theon 
to the practice of tanning, dyeing and fatliquory 
of leather. The use of materials having the prop, 
particle-size distribution to correspond with (| 
pore-size distribution in the leather might lead » 
better results in these processes. If small py. 
ticles are used along with larger ones, the great 
proportion of the smaller particles eventually y) 
be distributed into the small capillaries and ;) 
large particles will be in the large pores. Then. 
fore, if an oil having the correct size distributiy 
is used in fat-liquoring, a larger proportion of th 
surface should become lubricated. Likewise, \ 
tanning, a more uniform distribution of th 
tanning material would take place. 


Mole 
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oncentration of the Isotopes of Mercury by Free 
Evaporation in a 10-Cell Counter-Current Reflux Still’ 


By A. Keith Brewer*® and Samuel L. Madorsky 


Concentration of the isotopes of mercury was carried out in a counter-current molecular 


still consisting of 10 stages. 


Distillation temperatures varied from 55° to 100° C. 


A single- 


stage molecular still was operated, also with mercury, at temperatures varying from 108,1° 


to 128° C, 


ting with single-stage stills. 


I. Introduction 


Molecular distillation has been used by a num- 

r of investigators for concentrating the isotopes 
{ mereury. The technic employed in each in- 
ance was essentially the same. Brénsted and 
levesy [1] ° used a Dewar flask type still in which 
1 mercury evaporating in the heated outer flask 
as frozen out on the chilled inner flask. Maulli- 
en and Harkins [2] condensed the vapor on a 
ooled sloping glass roof, from which it ran into a 
ide reservoir. Harkins and Madorsky [3] used 
large steel annular trough with a A-shaped con- 
lensing roof. 

Each investigator used a single still with one 
vaporating and one condensing surface. In 
der to obtain an appreciable concentration it 
as necessary to repeat the process many times. 
This required an elaborate system for collecting, 
ecombining, and redistilling fractions. 

Wollner, Matchett, and Levine [4] describe a 
pultiple distillation apparatus which they used in 
urifying acetylated marihuana “red oil”. No 
jetails are given as to temperature, rate of dis- 
4 This paper covers experimental work carried out during the period June 

0, to April 3, 1941. Two statements describing the method and most 
the experimental work and dated June 24, 1940, and Nov. 25, 1940, respec- 

y, Were submitted to the Director of the National Bureau of Standards 

‘i recently it has been withheld from publication because of wartime 
urity restrictions. It will be included in vol. 6, Div. ITI, of the Manhattan 
roject Technical Series. 


Now at the U. 8. Navy Department. 
Numbers in brackets indicate the literature references at the end of this 
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The results agree with those obtained previously by other investigators opera- 


tillation, pressure inside the still, or as to whether 
a volatile fraction collected in the dry-ice-cooled 
trap. In view of the fact that organic liquids are 
subject to decomposition during high vacuum, 
high temperature, distillation, and that no pro- 
vision was made to cool the condensing part of the 
still, except contact with air, it is doubtful whether 
they operated under conditions of nonequilibrium 
distillation. 

The present paper describes an apparatus and 
technic in which a series of molecular stills have 
been arranged to operate under counter-current 
reflux. The new still is so designed that recom- 
bination of fractions is accomplished automat- 
ically by means of gravity feeds. The system acts 
as a reflux column in which each cell is a separate 
plate. 

The mechanism of operation, under total reflux, 
is shown graphically in figure 2, A. The vapor 
rising from the surface a; of the liquid in cell A, 
condenses on the sloping surface /; and runs down 
to drop into the adjacent cell A, higher up in the 
series. At the same time, the liquid level in cell 
A, is kept constant by means of a spillover, S, 
which allows the liquid to flow back into cell As. 
The flowback, in consequence, is exactly equal to 
the forward throughput. The over-all operation 
is such that the light fraction increases in concen- 
tration toward the upper end while the heavy 
fraction increases in concentration toward the 
lower end of the still. 
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II. Apparatus 


The evaporation apparatus is shown diagram- 
matically in figure 1. The evaporating reservoirs 
were machined from an iron plate, B. Each 
reservoir was 2 cm wide, 7.5 cm long, and 0.8 em 
deep at the shallow side. Spillovers S, figure 2, B, 
were cut at alternate ends for adjacent cells to 
provide thorough mixing of the liquid. The 
capacity of each cell was 9 ml, and that of tae 
entire system was 90 ml. The reservoirs were 
heated by a Nichrome coil imbedded in an alun- 
dum plate that was bolted directly to the reservoir 
block. 

The condensing roofs, C, figure 1, were made of 
iron wedges screwed to a supporting plate, D. 
The angle of slope could be varied by changing the 
wedges. The roofs were cooled by circulating 
alcohol at —18°C through an iron coil immersed 
in a pool of mercury, /, in contact with plate D. 
The temperature of the mercury was kept near 
0°C by a thermoregulator that operated the 
alcohol-circulating pump. 

Loss of mercury from the evaporators was 
prevented by glass plates provided with notches 
to allow proper evacuation of each cell. 

A high vacuum was obtained by enclosing the 
entire still in a glass cylinder with machined metal 
plates bolted across the ends; rubber gaskets 
coated with stopcock grease made the system 
vacuum-tight. The operating pressure, as read 
by a McLeod gage, was below 10~° mm Hg. 

Considerable difficulty was experienced in ob- 
taining roofs that did not allow most of the con- 
densed mercury to drop back into the same cell. 
In the first apparatus the machined surface of the 
iron wedges was used directly as the condensing 
surface; the angle was 25 degrees. A large frac- 
tion of the mercury dropped back into the same 
cell instead of running down into the next pool. 
Changing the slope of the wedges to 40 degrees 
made only a slight improvement. The condensing 
surfaces were then polished and blued at 300° C. 
The results were fairly satisfactory. A special 
study of condensing surfaces carried out with a 
single-cell apparatus showed that glass was by far 
the most effective material. Glass plates | mm 
thick were then cemented to the iron roofs. This 
change was effective in giving a complete forward 
transport of the condensing mercury. An angle 
of 33 degrees was satisfactory. 
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SECTION THRU A-A 
Figure 1.—Ten-cell cascade evaporation apparatus 
Purification and density determination of th 
mereury 


The mercury was first purified by aeration | 
24 hours in 10-percent nitric acid. It was the 
distilled three times, rejecting 10-percent fraction 
at the beginning and end of each distillation. 

The change in the isotope ratio was determine 
by the density method described by Mulliken a0 
Harkins [2] and Harkins and Madorsky 
Several modifications were introduced to incres 
the accuracy of the measurements. The picnot 
eter had a capacity of about 3 ml and held abou 
40 g of Hg. It was filled with the Hg sample * 
about 10° C, then placed in a thermostat at 3!¢ 
+0.003° C. The overflow of Hg was brushed 04 
with the sharp edge of a piece of paper exactly * 
the time when the reading on a thermometer 
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1oURE 2.—Diagramatic view of cascade distillation and 





apparatus for single-stage operation. 





liagram showing mechanism of cascade distillation; 7, section through 
s showing spillovers; C, Mulliken-type single-stage molecular still. 








ta given point on the downward temperature 
anation of the bath. At first a razor blade was 
sed to wipe off the mercury overflow, but it was 
und that the picnometer lost weight owing to 
raping of the glass. 

Weighings were made on a new chainomatic 
wlance, using a counterpois picnometer having 
he dimensions as close to the original as possible. 
lagnifying glasses were used to read the zero on 
he pointer scale and the divisions on the chain 
ial. An accuracy of about 2 parts per million 
as obtained. As the total weight of the Hg 
as 40 ¢, 2 ppm involved an accuracy of 0.08 mg. 
In order to check the density measurements 
gainst isotope abundance ratios, samples of 
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mercury of known density labeled No. 1 ordinary, 
No. 2 ordinary, and A, B, C, and D were sent to 
A. O. C. Nier, at the University of Minnesota, for 


mass spectrometric analysis. The ratios for 
“Hg to **Hg were measured. These ratios with 
comparative density measurements are given 


in table 1. 













Comparison between mass spectrograph analysis 
and density determinations of mercury 


TARLE 1. 


Change of 
density of 
mercury 


Ratio 


Sample 
He Hye 


ppm 


University of Minnesota, ordinary 1. 456 
No. 1, ordinary 1. 459 
No, 2, ordinary 1. 457 
A +10.9 1.450 
| B —11.2 1. 466 
Cc +30. 6 1. 438 
D 1 
, 1 


III. Temperature Measurements 


The temperature of the mercury in the cells 
was first determined by a thermometer and ther- 
mocouple inserted in plate B (fig. 1). It was 
evident, however, that the power input was much 
too small to account for the evaporation of Hg 
at a rate corresponding to the plate temperature. 
A Pt-Pt-Rh thermocouple inserted in the mercury 
gave incorrect readings because of stray contact 
potentials. In consequence, it was found neces- 
sary to estimate the temperature of the evap- 
orating Hg surface from the power input into the 
heater, after correcting for radiation and conduc- 
tion losses. 


IV. Operational Procedure 


After the cell had been filled with mercury, the 
system was closed and evacuated to a sticking 
vacuum. In all the runs the operation was con- 
tinuous. The mercury cooling bath was main- 
tained at about 0° C. 

Runs II, III, and IV were made with the 40- 
degree polished and blued roofs; the average 
distance from the evaporating surface was 19 mm. 
The duration was 2, 5, and 3 days, respectively. 
Run V+ VI was carried out for the first 6 days 
with blued iron roofs; the experiment was then 
stopped, and the glass-covered roofs were installed. 
The run was then continued for another 6 days. 
The slope angle was 33 degrees from the horizontal, 
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and the average distance from the evaporating 
surface was 15 mm. Runs VII, VIII, and IX 
were made with the same glass-covered roofs and 
lasted 3, 5, and 2 days, respectively. In run VII 
the mercury in the two lowest cells became mixed 
in sample taking. The conditions in run X were 
the same as in run LX, except that the temperature 
was higher. Soon after this run was started the 
glass plate fell off the roof in cell 5, and most 
of the condensate in this cell dropped back without 
reaching cell 6. In consequence, the lower four 
cells and the upper six cells in run X operated as 
separate systems. 


V. Throughput 


The throughput in the present still is determined 
by the quantity of mercury carried forward from 
cell to cell. As the stills were operated under 
total reflux, this is necessarily the same as the 
quantity of mercury spilling back. The forward 
movement of mercury is determined by three 
factors: (1) the rate of escape from the liquid 
surface, (2) the fraction condensed on the roofs, 
and (3) the fraction of condensate dropping back 
into the same cell. 


The rates of evaporation of Hg in grams 
square centimeter per hour were calculated yj 
the aid of Langmuir’s equation [5] 


Oe 
m= Ps) 5. RT” 


where 
m=(g of Hg evaporated per cm?) per secon 
P=pressure in dynes per em?® 
M=molecular weight of Hg 
R=gas constant=8.3 10" ergs per 1 deg | 
T=absolute temperature. 


The fraction of the mereury condensing on th 
roofs cannot be told with certainity. Assumip 
that the direction of escape of molecules from ths 
mercury surface is given by the cosine law, 4 
proximately 75 percent of the molecules leavin 
should strike the sloping surface. As this surfy 
is cooled to 0° C condensation thereon should | 
very nearly complete. 


VI. Results 


The pertinent data for nine consecutive mw 
are given in table 2. 


TaBLe 2.— Concentration of isotopes of mercury in the 10-cell apparatus 





Rate of 
evapora- 
tion 


Vapor 
pressure 


Tempera- 


_— 
Tim cure 


Experiment No 


mm He § (g/em*)/hr 

0.017 2. 95 
O48 75 
OAS .75 
O47 26 
067 
. 067 
067 


| Days 
I 


IV 


— 
eown won 


Total 
change 
in density 


Number 
of plates 


Angle to 


Kind of roof heriees 


Average,h 


Degrees 

Iron rT) 
do - 0 

do a 
Iron and glass , 33 
Glass 33 
. do 3 
do 33 
..do 33 


ppm 
21.8 
45.9 
vt 
129.4 
78. 2 
161.6 
85.7 


vocovwvnwece 


78 





The change in density of the mercury in the 
individual cells for each run is given in figures 3 
and 4. 

Runs II and IV are of the type to be expected in 
a reflux column in the initial stages of operation. 
The concentration predominates at the ends with 
very little effect in the central portion. Run III 
is characterized by a nonuniformity of transport 
between cells. All three runs were made with 
iron roofs in which a large fraction of the condensed 
mercury fell back into the cells from which it 
evaporated. As some of the roofs behaved more 
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efficiently than others, the erratic shape of line 
III is to be expected. 

Run V+ VI shows a‘uniform rate of change 1! 
density from cells 5 to 10. This is indicative of 
equilibrium having been established in that per- 
tion of the column. In the lower cells from 4 
down, the difference in density decreases slight! 
toward the bottom. This dropping off in the 
efficiency in the lower cells was expected as drp- 
lets of mercury were regularly observed to form 
in the upper part of the column and run along the 
walls to fall principally into cell 1. The constat! 
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that might be present would of necessity be 
lighter than mercury and at the same time have a 
lower vapor pressure, it would soon concentrate in 
the lowest cell and lower the density of the Hq in 
that cell. The possibility of impurities in the 
mercury, however, has been carefully considered 
and shown to be absent. The straight line ob- 
tained in run VIII is also definite proof for the 
absence of impurities. 

Runs VII, VIII, 1X, and X were made with 
glass plates on the condensing roofs. Visual ob- 
servation showed most of the mercury droplets 
falling into the next cell higher up. The quantity 








mn th | | T | T . . 
abs of mercury running down the walls and dropping 
suming : 
ap «COP 1. , into the lower cells was small. 
OM Tt a , r 

she . 2 oS ee a. Runs VII and IX are of the character to be 
W, ay CELL NUMBER FROM LOWEST TO HIGHEST ——~ ra hy : 
, expected when equilibrium has not been attained 
CAVIEEE cone 3.—Results obtained in runs II, 111, IV, and V_ between the various cells. The contour of the 
surf lus VI where polished ste sw TT indi 
_— plus VI where polished steel roofs were used. curve for run VII indicates that the throughput 


ual | {1,2 days at 85°; TEL, 8 days at 70°, IV, 3 days at 70°; V plus VI, 12days ig not uniform between the various cells as equi- 
librium is much nearer attained in the top three 
cells than in the remainder. 

Run VIII is characteristic of a counter-current 
system operating under total reflux when equi- 
librium has been established between all cells. 

The results obtained at 100° C are illustrated in 
run X._ As the glass plate on cell 5 fell off during 
operation, the data are not complete for the entire 
system. They do show, however, that cells from 
6 to 10 were in equilibrium. The separation per 
plate, as shown by figure 4, is 17.3 ppm. 

Figure 5 shows that the separation values are 
in general agreement with those to be expected 
from single-stage stills. As the exact temperature 
of the evaporating surface was not known in any 
instance, a more precise agreement cannot be 





expected. 








| 
arcs ta VII. Free Evaporation of Mercury at High 
Bs nsf Temperatures 





CELL NUMBER FROM LOWEST TO HIGHEST——= 
Results obtained in runs VII, VIII, 1X, and X, 


where glass-covered roofs were used. 


rT | 8 


The 10-cell cascade system was not designed for 
a study of free evaporation at temperatures above 
about 100° C. As such knowledge is essential in 
working out more efficient still designs, tests at 
temperatures above 100° C were made in a single- 
hildition of light mereury from the top cells to the — stage glass molecular still similar to the one used 
ottom cell gives rise to the observed decrease in by Mulliken and Harkins [2]. The still is il- 
rate lustrated in figure 2, C. 
It should be pointed out that a small difference The capacity of the evaporator was 200 g of 
bi density between the two lowest cells might indi- Hg (fig. 2, C), the evaporating surface was 10 
‘ule Impurity in the mercury. As any impurity cm’, and the condensing surface had an angle of 


lays at 75°; VIII, 5 days at 75°; IX, 2 days at 75°; X, 2 days at 100°. 
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30 degrees. The roof was cooled by an alcohol-— VIII. Comparison with previous work 
dry-ice mixture, and the evaporator was heated 
with an oxy-gas blast-lamp. The temperature of 
the mercury evaporating surface was calculated 
from the quantity of mercury evaporated per unit 
time, using Langmuir’s equation. 

Four runs at different temperatures were made. 
In the final run the fu'l foree of the oxy-gas flame 
was used. The cut, C, in each run was approxi- oral Sculeteniecaiaal 


mately 2. The pertinent data for these experi- © HARKINS AND MADORS«y 
@ MULLIKEN AND HARKINS 

© 10-CELL CASCADE 

© SINGLE STAGE 


_~— 


It is interesting to compare the present resy)is 
obtained in the 10-cell cascade apparatus with thp 
single-cell work at the University of Chicagy 
Harkins and Madorsky [3] obtained about 1459 
ppm decrease in density for the light fraction whe 
a cut of 2 was made. This value corresponds 


ments are given in table 3. 





Taste 3.—Concentration of isotopes of mercury in a single- 
stage apparatus 








Change | ,, 
7 | ; i. Change 
Experiment = Vapor | Rate of in density in dom> 
- rime eture pres- evapo- of light sit cal 
sure | ration | fraction | *' ; per 
| per cut | Plate 


lemper- 


| (g/cm *) 
min , mm He hr ppm ppm 
11 0.417 62.8 10. 50 15.15 
8 570 87.3 10.40 | 15.01 
6 770 113.7 8. 05 11. 62 
4 075 161.0 4.80 6. 


+ 


CHANGE OF DENSITY IN PPM PER PLATE 














| 
The change in density, AD, of the light fractions | | 
70 


was obtained by dividing the total difference in 60 80 8690 = 00 


density between the light and heavy fractions TEMPERATURE °C 
7 D The ec » » ity » » y 4 . - 
by —_ Phe change of dle nsity pe r plate ? AD , Was Figure 5. Comparison between the present work and 


computed from the equation given by Mulliken by other investigators 

and Harkins [2]: 

16.50/0.69315=23.80 ppm per plate. In the work 
of Mulliken and Harkins [2], a cut of 2 gave: 
decrease in density for the light fractions of |’ 
ppm. This corresponds to 16/0.69315=23.\8 


,ne 
AD=AD' a> 


where AD and AD’ are expressed in parts per * : 
million. Then ppm per plate. Table 4 and figure 5 show a con- 


parison between the work done at the Universi 
AD’ = AD/0.69315. 2) of Chicago and the present work. 


TABLE 4.—Concentration of isotopes of mercury: comparison with previous work 


Change in 

Vapor . density of 
pressure light frae- 
tion per cut 


Change 
in density 
per plate 


Number Rate of 
Reference of cells evapora- | Calculated 
in still tion 


(a/em®)/hr Tc mm Hg ppm ppm 
Harkins and Madorsky 2.37 52. 0.015 SJ 16. 50 23 
Mulliken and Harkins +00 59.5 . 025 q 16. 00 
Present work, experiment VIII 10. 27 . 07 
Present work, experiment XI 62.8 417 
Present work, experiment XII | 87.3 570 
Present work, experiment XIII | 113. 770 
Present work, experiment XIV 161.0 1.075 


10. % 
10. 40 
8.05 
4.08 


0 tS to 
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IX. Conclusions 


The results presented in this paper demonstrate 
at a series of individual molecular stills can be 
ranged in counter-current reflux in such a man- 
br that the throughput between cells is entirely 
otrolled by gravity feed. The advantages of 


» present reflux still over the single-stage still 


» many: (1) The quantity of material required 
vive a desired separation is materially smaller, 
the holdup is vastly lower, (3) the time required 
ra given separation is very much shorter, (4) 

labor involved is less; this is illustrated by the 
et that to realize the separation obtained in one 
wration in the present 10-compartment still, 55 
dividual recombinations and distillations would 
ive to be made with single-stage stills, (5) the 
wration can be made continuous; the material 
» be concentrated can be fed into the system at 
e end or at the midpoint, if both light and 
eavy fractions are of interest, and the concen- 


trate can be withdrawn continuously at the ends, 
and (6) The entire separation can be run without 
exposing the material to the atmosphere. This 
is particularly desirable for liquids that are sus- 
ceptible to decomposition upon exposure to air. 


The authors express their appreciation to J. W. 
Westhaver for his many basic contributions to this 
research. 
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